【BZOJ 1877】 [SDOI2009]晨跑
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。
每个点拆点限制流量只能为1,以防止重复经过,然后套最小费用流,中间统计增广几次
1 #include<cstdio> 2 #include<cstring> 3 using namespace std; 4 const int N=25100,inf=1000000; 5 struct ee{int to,next,f,w;}e[N*2]; 6 int S,T,cnt=1,n,k,ans,timer,m,u,v,w; 7 int head[N],dis[N],pre[N],q[N]; 8 bool inq[N]; 9 void ins(int u,int v,int f,int w){ 10 e[++cnt].to=v,e[cnt].next=head[u],e[cnt].f=f,e[cnt].w=w,head[u]=cnt; 11 e[++cnt].to=u,e[cnt].next=head[v],e[cnt].f=0,e[cnt].w=-w,head[v]=cnt; 12 } 13 14 bool spfa(){ 15 for (int i=1;i<=T;i++) dis[i]=inf; 16 int h=0,t=1; 17 q[t]=S;dis[S]=0;inq[S]=1; 18 while (h!=t){ 19 int now=q[++h];if(h==2501) h=0; 20 for (int i=head[now];i;i=e[i].next){ 21 int v=e[i].to; 22 if (dis[v]>dis[now]+e[i].w&&e[i].f){ 23 dis[v]=dis[now]+e[i].w; 24 pre[v]=i; 25 if (!inq[v]){ 26 q[++t]=v;if (t==2501) t=0; 27 inq[v]=1; 28 } 29 } 30 } 31 inq[now]=0; 32 } 33 if (dis[T]==inf) return 0; 34 return 1; 35 } 36 37 void updata(){ 38 int tmp=T; 39 while (tmp!=S){ 40 int l=pre[tmp],v=e[l].to; 41 e[l].f-=1;e[l^1].f+=1; 42 tmp=e[l^1].to; 43 } 44 ans+=dis[T]; 45 } 46 47 int main(){ 48 scanf("%d%d",&n,&m); 49 S=0;T=n*2+1; 50 ins(S,1,inf,0);ins(2*n,T,inf,0); 51 for (int i=2;i<n;i++) ins(i,i+n,1,0); 52 ins(1,1+n,inf,0);ins(n,n+n,inf,0); 53 for (int i=1;i<=m;i++){ 54 scanf("%d%d%d",&u,&v,&w); 55 ins(u+n,v,1,w); 56 } 57 while (spfa()) { 58 timer++; 59 updata(); 60 } 61 printf("%d %d",timer,ans); 62 }