1877. [SDOI2009]晨跑【費用流】
阿新 • • 發佈:2018-03-31
clu use post AC while tdi nod %d 最短
第一行:兩個數N,M。表示十字路口數和街道數。
接下來M行,每行3個數a,b,c,表示路口a和路口b之間有條長度為c的街道(單向)。
N ≤ 200,M ≤ 20000。
1 2 1
1 3 1
2 4 1
3 4 1
4 5 1
4 6 1
2 5 5
3 6 6
5 7 1
6 7 1
唯一和模板不同的就是我們要限制每個點只能到一次。
這一看就是常規拆點啊
所以
我們將點裂成兩個,然後在兩點間連一個容量為1費用為0的邊,用來限制每個點只能走一次
然後再在x+n和y中的連邊,把容量設為1,費用設為邊長
跑一邊從(n+1)到n的最小費用最大流(因為1節點可以重復走所以要n+1)
最大流和最小費用即為答案。
Description
Elaxia最近迷戀上了空手道,他為自己設定了一套健身計劃,比如俯臥撐、仰臥起坐等 等,不過到目前為止,他 堅持下來的只有晨跑。 現在給出一張學校附近的地圖,這張地圖中包含N個十字路口和M條街道,Elaxia只能從 一 個十字路口跑向另外一個十字路口,街道之間只在十字路口處相交。Elaxia每天從寢室出發 跑到學校,保證寢室 編號為1,學校編號為N。 Elaxia的晨跑計劃是按周期(包含若幹天)進行的,由於他不喜歡走重復的路線,所以 在一個周期內,每天的晨跑路線都不會相交(在十字路口處),寢室和學校不算十字路 口。Elaxia耐力不太好, 他希望在一個周期內跑的路程盡量短,但是又希望訓練周期包含的天 數盡量長。 除了練空手道,Elaxia其他時間 都花在了學習和找MM上面,所有他想請你幫忙為他設計 一套滿足他要求的晨跑計劃。Input
Output
兩個數,第一個數為最長周期的天數,第二個數為滿足最長天數的條件下最短的路程長 度。
Sample Input
7 101 2 1
1 3 1
2 4 1
3 4 1
4 5 1
4 6 1
2 5 5
3 6 6
5 7 1
6 7 1
Sample Output
2 11 一道近乎最小費用最大流的模板題唯一和模板不同的就是我們要限制每個點只能到一次。
這一看就是常規拆點啊
所以
我們將點裂成兩個,然後在兩點間連一個容量為1費用為0的邊,用來限制每個點只能走一次
然後再在x+n和y中的連邊,把容量設為1,費用設為邊長
跑一邊從(n+1)到n的最小費用最大流(因為1節點可以重復走所以要n+1)
最大流和最小費用即為答案。
#include<iostream> #include<cstring> #include<cstdio> #include<queue> using namespace std; struct node { int Flow; int Cost; int to; int next; }edge[100010]; queue<int>q; int INF,dis[601],pre[601]; int head[601],num_edge; bool used[601]; void add(int u,int v,int l,int c) { edge[++num_edge].to=v; edge[num_edge].Flow=l; edge[num_edge].Cost=c; edge[num_edge].next=head[u]; head[u]=num_edge; } bool SPFA(int s,int e) { memset(pre,-1,sizeof(pre)); memset(dis,0x7f,sizeof(dis)); q.push(s); dis[s]=0; used[s]=true; while (!q.empty()) { int x=q.front(); q.pop(); for (int i=head[x];i!=0;i=edge[i].next) if (edge[i].Flow>0 && dis[x]+edge[i].Cost<dis[edge[i].to]) { dis[edge[i].to]=dis[x]+edge[i].Cost; pre[edge[i].to]=i; if (!used[edge[i].to]) { used[edge[i].to]=true; q.push(edge[i].to); } } used[x]=false; } return (dis[e]!=INF); } void MCMF(int s,int e) { int d=INF,Ans=0,Fee=0; while (SPFA(s,e)) { for (int i=e;i!=s;i=i=edge[((pre[i]-1)^1)+1].to) d=min(edge[pre[i]].Flow,d); for (int i=e;i!=s;i=i=edge[((pre[i]-1)^1)+1].to) { edge[pre[i]].Flow-=d; edge[((pre[i]-1)^1)+1].Flow+=d; } Ans+=d; Fee+=dis[e]*d; } printf("%d %d",Ans,Fee); } int main() { memset(&INF,0x7f,sizeof(INF)); int n,m,u,v,l; scanf("%d%d",&n,&m); for (int i=1;i<=n;++i) { add(i,i+n,1,0); add(i+n,i,0,0); } for (int i=1;i<=m;++i) { scanf("%d%d%d",&u,&v,&l); add(u+n,v,1,l); add(v,u+n,0,-l); } MCMF(n+1,n); }
1877. [SDOI2009]晨跑【費用流】