【BZOJ】【1877】【SDOI2009】晨跑
网络流/费用流
费用流入门题……根本就是模板题好吗!
拆点搞定度数限制,也就是每个点最多经过一次……源点汇点除外。
1 /************************************************************** 2 Problem: 1877 3 User: Tunix 4 Language: C++ 5 Result: Accepted 6 Time:900 ms 7 Memory:5972 kb 8 ****************************************************************/ 9 10 //BZOJ 1877 11 #include<vector> 12 #include<cstdio> 13 #include<cstring> 14 #include<cstdlib> 15 #include<iostream> 16 #include<algorithm> 17 #define rep(i,n) for(int i=0;i<n;++i) 18 #define F(i,j,n) for(int i=j;i<=n;++i) 19 #define D(i,j,n) for(int i=j;i>=n;--i) 20 #define pb push_back 21 using namespace std; 22 inline int getint(){ 23 int v=0,sign=1; char ch=getchar(); 24 while(ch<'0'||ch>'9'){ if (ch=='-') sign=-1; ch=getchar();} 25 while(ch>='0'&&ch<='9'){ v=v*10+ch-'0'; ch=getchar();} 26 return v*sign; 27 } 28 const int N=1010,M=200000,INF=~0u>>2; 29 typedef long long LL; 30 /******************tamplate*********************/ 31 int n,m,ans,flow; 32 struct edge{int from,to,v,c;}; 33 struct Net{ 34 edge E[M]; 35 int head[N],next[M],cnt; 36 void ins(int x,int y,int z,int c){ 37 E[++cnt]=(edge){x,y,z,c}; 38 next[cnt]=head[x]; head[x]=cnt; 39 } 40 void add(int x,int y,int z,int c){ 41 ins(x,y,z,c); ins(y,x,0,-c); 42 } 43 int from[N],Q[M],d[N],S,T; 44 bool inq[N]; 45 bool spfa(){ 46 int l=0,r=-1; 47 F(i,1,T) d[i]=INF; 48 d[S]=0; Q[++r]=S; inq[S]=1; 49 while(l<=r){ 50 int x=Q[l++]; 51 inq[x]=0; 52 for(int i=head[x];i;i=next[i]) 53 if(E[i].v>0 && d[x]+E[i].c<d[E[i].to]){ 54 d[E[i].to]=d[x]+E[i].c; 55 from[E[i].to]=i; 56 if (!inq[E[i].to]){ 57 Q[++r]=E[i].to; 58 inq[E[i].to]=1; 59 } 60 } 61 } 62 return d[T]!=INF; 63 } 64 void mcf(){ 65 int x=INF; 66 for(int i=from[T];i;i=from[E[i].from]) 67 x=min(x,E[i].v); 68 for(int i=from[T];i;i=from[E[i].from]){ 69 E[i].v-=x; 70 E[i^1].v+=x; 71 } 72 flow+=x; 73 ans+=x*d[T]; 74 } 75 void init(){ 76 n=getint(); m=getint(); cnt=1; 77 S=1; T=n+n; 78 F(i,1,n) add(i,i+n,1,0); 79 add(1,1+n,INF,0); 80 add(n,n+n,INF,0); 81 int x,y,z; 82 F(i,1,m){ 83 x=getint(); y=getint(); z=getint(); 84 add(x+n,y,1,z); 85 } 86 while(spfa()) mcf(); 87 printf("%d %d\n",flow,ans); 88 } 89 }G1; 90 91 int main(){ 92 #ifndef ONLINE_JUDGE 93 freopen("1877.in","r",stdin); 94 freopen("1877.out","w",stdout); 95 #endif 96 G1.init(); 97 return 0; 98 } 99
1877: [SDOI2009]晨跑
Time Limit: 4 Sec Memory Limit: 64 MBSubmit: 1371 Solved: 712
[Submit][Status][Discuss]
Description
Elaxia最近迷恋上了空手道,他为自己设定了一套健身计划,比如俯卧撑、仰卧起坐等
等,不过到目前为止,他坚持下来的只有晨跑。
现在给出一张学校附近的地图,这张地图中包含N个十字路口和M条街道,Elaxia只能从
一个十字路口跑向另外一个十字路口,街道之间只在十字路口处相交。Elaxia每天从寝室出发
跑到学校,保证寝室编号为1,学校编号为N。
Elaxia的晨跑计划是按周期(包含若干天)进行的,由于他不喜欢走重复的路线,所以
在一个周期内,每天的晨跑路线都不会相交(在十字路口处),寝室和学校不算十字路
口。Elaxia耐力不太好,他希望在一个周期内跑的路程尽量短,但是又希望训练周期包含的天
数尽量长。
除了练空手道,Elaxia其他时间都花在了学习和找MM上面,所有他想请你帮忙为他设计
一套满足他要求的晨跑计划。
Input
第一行:两个数N,M。表示十字路口数和街道数。
接下来M行,每行3个数a,b,c,表示路口a和路口b之间有条长度为c的街道(单向)。
Output
两个数,第一个数为最长周期的天数,第二个数为满足最长天数的条件下最短的路程长
度。
Sample Input
7 10
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 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
HINT
对于30%的数据,N ≤ 20,M ≤ 120。
对于100%的数据,N ≤ 200,M ≤ 20000。