bzoj 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

2
5

HINT

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

我猜是因为bzoj没有check。。所以把最小割的方案输出给去掉了。

这题可以先跑一遍spfa。然后就可以求出最短路图了。然后再求图中最小割即可。

当时数组开小了。WA和RE若干。【复制模板是种病】

#include<queue>
#include<cstdio>
#include<string>
#include<cstring>
using namespace std;
inline int min(int x,int y)
{
     if(x<y)
          return x;
     return y;
}
int head[100001],exhead[100001];
struct map
{
	 int f;
	 int s,t;
     int next;
}a[400001],exa[400001];
int edge,exedge;
int p;
int q[400001],d[400001];
inline void add(int s,int t,int f)
{
     a[edge].next=head[s];
     head[s]=edge;
     a[edge].s=s;
     a[edge].t=t;
     a[edge].f=f;
}
inline void exadd(int s,int t,int f)
{
     exa[exedge].next=exhead[s];
     exhead[s]=exedge;
     exa[exedge].s=s;
     exa[exedge].t=t;
     exa[exedge].f=f;
}
inline bool bfs()
{
     int l=0,r=0;
     memset(q,0,sizeof(q));
     r++;
     q[r]=1;
     memset(d,-1,sizeof(d));
     d[1]=0;
     int i,k;
     while(l<r)
     {
     	  l++;
          int k=q[l];
          for(i=head[k];i!=0;i=a[i].next)
          {
               if(a[i].f>0&&d[a[i].t]==-1)
               {
                    d[a[i].t]=d[k]+1;
                    r++;
                    q[r]=a[i].t;
               }
          }
     }
     if(d[p]>=0)
          return true;
     return false;
}
inline int dfs(int k,int s)
{
     if(k==p)
          return s;
     int t=s;
     int i;
     for(i=head[k];i!=0;i=a[i].next)
     {
          if(d[a[i].t]==d[k]+1&&a[i].f>0)
          {
               int xx=dfs(a[i].t,min(s,a[i].f));
               a[i].f-=xx;
               if(i%2==0)
                    a[i-1].f+=xx;
               else
                    a[i+1].f+=xx;
               s-=xx;
          }
     }
     return t-s;
}
inline int maxflow()
{
     int s=0;
     while(bfs())
          s+=dfs(1,2100000000);
     return s;
}
int dis[1001];
bool v[1001];
inline void spfa()
{
	 memset(dis,127/3,sizeof(dis));
	 memset(v,false,sizeof(v));
     queue <int>Q;
     while(!Q.empty())
          Q.pop();
     Q.push(1);
     dis[1]=0;
     v[1]=true;
     while(!Q.empty())
     {
          int d=Q.front();
          v[d]=false;
          Q.pop();
          int i;
          for(i=exhead[d];i!=0;i=exa[i].next)
          {
               int t=exa[i].t;
               if(dis[d]+exa[i].f<dis[t])
               {
                    dis[t]=dis[d]+exa[i].f;
                    if(!v[t])
                    {
                         Q.push(t);
                         v[t]=true;
                    }
               }
          }
     }
}
int sum;
inline void bfs(int d)
{
     int i;
     for(i=head[d];i!=0;i=a[i].next)
     {
          int t=a[i].t;
          if(!v[t]&&a[i].f>0)
          {
               sum++;
          	   v[t]=true;
          }
     }
}
int dk[400001];
int main()
{
     int T;
     //freopen("route.in","r",stdin);
     //freopen("route.out","w",stdout);
    // scanf("%d",&T);
    // while(T>0)
    // {
     //     T--;
          memset(head,0,sizeof(head));
          memset(a,0,sizeof(a));
          memset(exhead,0,sizeof(exhead));
          memset(exa,0,sizeof(exa));
          edge=0;
          exedge=0;
		  int n,m;
          scanf("%d%d",&n,&m);
          int i;
          int s,t,x,c;
          for(i=1;i<=m;i++)
          {
               scanf("%d%d%d%d",&s,&t,&x,&c);
               exedge++;
               exadd(s,t,x);
               dk[exedge]=c;
               exedge++;
               exadd(t,s,x);
               dk[exedge]=c;
          }
          spfa();
          printf("%d\n",dis[n]);
          for(i=1;i<=exedge;i++)
          {
               if(dis[exa[i].s]+exa[i].f==dis[exa[i].t])
               {
                    edge++;
                    x=dk[i];
                    add(exa[i].s,exa[i].t,x);
                    edge++;
                    add(exa[i].t,exa[i].s,0);
               }
          }
          p=n;
          int ans=maxflow();
          printf("%d\n",ans);
    // }
     return 0;
}