bzoj 1266 1266: [AHOI2006]上学路线route

Description

可可和卡卡家住合肥市的东郊，每天上学他们都要转车多次才能到达市区西端的学校。直到有一天他们两人参加了学校的信息学奥林匹克竞赛小组才发现每天上学的乘车路线不一定是最优的。可可：“很可能我们在上学的路途上浪费了大量的时间，让我们写一个程序来计算上学需要的最少时间吧！” 合肥市一共设有N个公交车站，不妨将它们编号为1…N的自然数，并认为可可和卡卡家住在1号汽车站附近，而他们学校在N号汽车站。市内有M条直达汽车路线，执行第i条路线的公交车往返于站点pi和qi之间，从起点到终点需要花费的时间为ti。(1<=i<=M, 1<=pi, qi<=N) 两个人坐在电脑前，根据上面的信息很快就编程算出了最优的乘车方案。然而可可忽然有了一个鬼点子，他想趁卡卡不备，在卡卡的输入数据中删去一些路线，从而让卡卡的程序得出的答案大于实际的最短时间。而对于每一条路线i事实上都有一个代价ci：删去路线的ci越大卡卡就越容易发现这个玩笑，可可想知道什么样的删除方案可以达到他的目的而让被删除的公交车路线ci之和最小。 [任务] 编写一个程序： ? 从输入文件中读取合肥市公交路线的信息； ? 计算出实际上可可和卡卡上学需要花费的最少时间； ? 帮助可可设计一个方案，删除输入信息中的一些公交路线，使得删除后从家到学校需要的最少时间变大，而被删除路线的ci和最小；向输出文件输出答案。

Input

输入文件中第一行有两个正整数N和M，分别表示合肥市公交车站和公交汽车路线的个数。以下M行，每行（第i行，总第(i+1)行）用四个正整数描述第i条路线：pi, qi, ti, ci；具体含义见上文描述。

Output

输出文件最多有两行。第一行中仅有一个整数，表示从可可和卡卡家到学校需要的最短时间。第二行输出一个整数C,表示Ci之和

Sample Input

6 7
1 2 1 3
2 6 1 5
1 3 1 1
3 4 1 1
4 6 1 1
5 6 1 2
1 5 1 4

Sample Output

HINT

2<=N<=500, 1<=M<=124 750, 1<=ti, ci<=10 000
合肥市的公交网络十分发达，你可以认为任意两个车站间都可以通过直达或转车互相到达，当然如果在你提供的删除方案中，家和学校无法互相到达，那么则认为上学需要的最短为正无穷大：这显然是一个合法的方案。

Source

思路：跑spfa+网络流最小割。跑正反两边spfa，然后拎出一条边，判断边权+左端点到1号节点的最短路+右端点到n号节点的最短路（离1号节点近的点为左端点，反之为右端点）是否等于1到n的最短路，如果相等左右端点连边，权值为这条边代价。然后跑dinic。

错因：网络流建图错误，不应该是跑spfa时建图，因为spfa有松弛操作，可能边数建多了。

#include<queue> #include<vector> #include<cstdio> #include<cstring> #include<iostream> #include<algorithm> #define N 501 #define M 124751 using namespace std; queue<int>que; struct nond{ int x,y,z,zz; }edge[M]; int n,m;; int tot,tot1=1; int scr,decc,ans,maxn; int lev[M],cur[M]; int dis[N],dis1[N],vis[N]; int to1[M*2],net1[M*2],cap1[M*2],head1[M*2]; int to[M*2],net[M*2],cap[M*2],val[M*2],head[M*2]; void add1(int u,int v,int w){ to1[++tot1]=v;net1[tot1]=head1[u];cap1[tot1]=w;head1[u]=tot1; to1[++tot1]=u;net1[tot1]=head1[v];cap1[tot1]=0;head1[v]=tot1; } void add(int u,int v,int w,int z){ to[++tot]=v;net[tot]=head[u];cap[tot]=w;val[tot]=z;head[u]=tot; to[++tot]=u;net[tot]=head[v];cap[tot]=w;val[tot]=z;head[v]=tot; } void spfa(int s){ memset(vis,0,sizeof(vis)); memset(dis,0x3f,sizeof(dis)); while(!que.empty()) que.pop(); que.push(s); vis[s]=1;dis[s]=0; while(!que.empty()){ int now=que.front(); que.pop(); vis[now]=0; for(int i=head[now];i;i=net[i]) if(dis[to[i]]>dis[now]+cap[i]){ dis[to[i]]=dis[now]+cap[i]; if(!vis[to[i]]){ vis[to[i]]=1; que.push(to[i]); } } } } void spfa1(int s){ memset(vis,0,sizeof(vis)); memset(dis1,0x3f,sizeof(dis)); while(!que.empty()) que.pop(); que.push(s); vis[s]=1;dis1[s]=0; while(!que.empty()){ int now=que.front(); que.pop(); vis[now]=0; for(int i=head[now];i;i=net[i]) if(dis1[to[i]]>dis1[now]+cap[i]){ dis1[to[i]]=dis1[now]+cap[i]; if(!vis[to[i]]){ vis[to[i]]=1; que.push(to[i]); } } } } bool bfs(){ for(int i=scr;i<=decc;i++){ lev[i]=-1; cur[i]=head1[i]; } while(!que.empty()) que.pop(); que.push(scr); lev[scr]=0; while(!que.empty()){ int now=que.front(); que.pop(); for(int i=head1[now];i;i=net1[i]) if(lev[to1[i]]==-1&&cap1[i]>0){ lev[to1[i]]=lev[now]+1; que.push(to1[i]); if(to1[i]==decc) return true; } } return false; } int dinic(int now,int flow){ if(now==decc) return flow; int rest=0,detal; for(int & i=cur[now];i;i=net1[i]) if(lev[to1[i]]==lev[now]+1&&cap1[i]){ detal=dinic(to1[i],min(flow-rest,cap1[i])); if(detal){ cap1[i]-=detal; cap1[i^1]+=detal; rest+=detal; if(rest==flow) break; } } if(rest!=flow) lev[now]=-1; return rest; } int main(){ scanf("%d%d",&n,&m); scr=1;decc=n; for(int i=1;i<=m;i++){ int u,v,w,z; scanf("%d%d%d%d",&u,&v,&w,&z); edge[i].x=u; edge[i].y=v; edge[i].z=w; edge[i].zz=z; add(u,v,w,z); } spfa(1); cout<<dis[n]<<endl; spfa1(n); for(int i=1;i<=m;i++){ int u=edge[i].x,v=edge[i].y; if(dis[u]>dis[v]) swap(u,v); if(dis[u]+dis1[v]+edge[i].z==dis[n]) add1(u,v,edge[i].zz); } while(bfs()) ans+=dinic(scr,0x7f7f7f7f); cout<<ans; }