BZOJ1877 [SDOI2009] 晨跑
题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=1877
Description
Elaxia最近迷恋上了空手道，他为自己设定了一套健身计划，比如俯卧撑、仰卧起坐等 等，不过到目前为止，他坚持下来的只有晨跑。 现在给出一张学校附近的地图，这张地图中包含N个十字路口和M条街道，Elaxia只能从 一个十字路口跑向另外一个十字路口，街道之间只在十字路口处相交。Elaxia每天从寝室出发 跑到学校，保证寝室编号为1，学校编号为N。 Elaxia的晨跑计划是按周期（包含若干天）进行的，由于他不喜欢走重复的路线，所以 在一个周期内，每天的晨跑路线都不会相交（在十字路口处），寝室和学校不算十字路 口。Elaxia耐力不太好，他希望在一个周期内跑的路程尽量短，但是又希望训练周期包含的天 数尽量长。 除了练空手道，Elaxia其他时间都花在了学习和找MM上面，所有他想请你帮忙为他设计 一套满足他要求的晨跑计划。
Input
第一行：两个数N,M。表示十字路口数和街道数。 接下来M行，每行3个数a,b,c，表示路口a和路口b之间有条长度为c的街道（单向）。
Output
两个数，第一个数为最长周期的天数，第二个数为满足最长天数的条件下最短的路程长度。
水到不行的裸费用流
顺带一提：最近一直写费用流，STOIt2都用费用流写果断炸，以为真的可以用费用流写，结果今天改了很久才发现根本不可能……
我智商不够，早日滚粗报平安
1 #include <iostream> 2 #include <cstdio> 3 #include <algorithm> 4 #include <cstring> 5 #include <queue> 6 #define rep(i,l,r) for(int i=l; i<=r; i++) 7 #define clr(x,y) memset(x,y,sizeof(x)) 8 #define travel(x) for(Edge *p=last[x]; p; p=p->pre) 9 using namespace std; 10 const int INF = 0x3f3f3f3f; 11 const int maxn = 210; 12 inline int read(){ 13 int ans = 0, f = 1; 14 char c = getchar(); 15 for(; !isdigit(c); c = getchar()) 16 if (c == '-') f = -1; 17 for(; isdigit(c); c = getchar()) 18 ans = ans * 10 + c - '0'; 19 return ans * f; 20 } 21 struct Edge{ 22 Edge *pre,*rev; int to,cap,cost; 23 }edge[50010],*last[maxn<<1],*cur[maxn<<1],*pt; 24 int n,m,x,y,z,S,T,flow,cost,d[maxn<<1]; 25 bool isin[maxn<<1],vis[maxn<<1]; 26 queue <int> q; 27 inline void add(int x,int y,int z,int w){ 28 pt->pre = last[x]; pt->to = y; pt->cap = z; pt->cost = w; last[x] = pt++; 29 pt->pre = last[y]; pt->to = x; pt->cap = 0; pt->cost = -w; last[y] = pt++; 30 last[x]->rev = last[y]; last[y]->rev = last[x]; 31 } 32 bool spfa(){ 33 clr(isin,0); isin[S] = 1; q.push(S); 34 clr(d,INF); d[S] = 0; 35 while (!q.empty()){ 36 int now = q.front(); q.pop(); isin[now] = 0; 37 travel(now){ 38 if (d[p->to] > d[now] + p->cost && p->cap > 0){ 39 d[p->to] = d[now] + p->cost; 40 if (!isin[p->to]) isin[p->to] = 1, q.push(p->to); 41 } 42 } 43 } 44 return d[T] != INF; 45 } 46 int dfs(int x,int flow){ 47 if (x == T || (!flow)) return flow; vis[x] = 1; int w = 0; 48 for(Edge *p = cur[x]; p && w < flow; p = p->pre){ 49 if (p->cap > 0 && (!vis[p->to]) && d[p->to] == d[x] + p->cost){ 50 int delta = dfs(p->to,min(p->cap,flow-w)); 51 p->cap -= delta; p->rev->cap += delta; w += delta; 52 cost += p->cost * delta; 53 if (p->cap) cur[x] = p; 54 } 55 } 56 return w; 57 } 58 void mincost(){ 59 while (spfa()){ 60 clr(vis,0); 61 rep(i,S,T) cur[i] = last[i]; 62 flow += dfs(S,INF); 63 } 64 } 65 int main(){ 66 n = read(); m = read(); clr(last,0); pt = edge; 67 rep(i,1,m){ 68 x = read(); y = read(); z = read(); 69 add(x+n,y,1,z); 70 } 71 rep(i,2,n-1) add(i,i+n,1,0); 72 add(1,n+1,INF,0); 73 add(n,n<<1,INF,0); 74 S = 1; T = n << 1; 75 mincost(); 76 printf("%d %d\n",flow,cost); 77 return 0; 78 }
