[bzoj1266][网络流][spfa]上学路线

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
合肥市的公交网络十分发达，你可以认为任意两个车站间都可以通过直达或转车互相到达，当然如果在你提供的删除方案中，家和学校无法互相到达，那么则认为上学需要的最短为正无穷大：这显然是一个合法的方案。

题解

我我我我我**网络流建图
贡献了1版wa的我。。
首先观察题目，不难想第一问就是跑一个裸的spfa嘛。。
那么为什么会出来第二问？？想一想最短路径可能不止一条嘛。原来想往回走路删最短路最小边的希望破灭。那怎么办？
对1~n的所有最短路径建网络流图，跑一遍dinic，最大流就是第二问答案
为啥？最小割=最大流啊
至于我如何wa的。。
建图的ins少打了两行

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
struct node
{
    int x,y,next,c,other,d; 
}a[310000];
int last[310000],len;
void ins(int x,int y,int c,int d)
{
    len++;
    a[len].x=x;a[len].y=y;a[len].c=c;a[len].d=d;
    a[len].next=last[x];last[x]=len;
}
int d[310000];
int list[310000];
bool v[310000];
int n,m;
int addx[310000],addy[310000],addd[310000];
void insert(int k)
{
    int k1,k2;
    int x=addx[k],y=addy[k],c=addd[k];
    len++;k1=len;//!!!!记住继承啊啊啊啊啊 
    a[len].x=x;a[len].y=y;a[len].c=c;
    a[len].next=last[x];last[x]=len;

    len++;k2=len;
    a[len].x=y;a[len].y=x;a[len].c=0;
    a[len].next=last[y];last[y]=len;

    a[k1].other=k2;a[k2].other=k1;
}
int h[310000];
int head,tail;
int st,ed;
bool bt_h()
{
    memset(h,0,sizeof(h));
    h[st]=1;list[1]=1;
    head=1;tail=2;
    while(head!=tail)
    {
        int x=list[head];
        for(int k=last[x];k;k=a[k].next)
        {
            int y=a[k].y;
            if(h[y]==0 && a[k].c>0)
            {
                h[y]=h[x]+1;
                list[tail++]=y;
            }
        }
        head++;
    }
    if(h[ed]>0)return true;
    return false;
}
int findflow(int x,int f)
{
    if(x==ed)return f;
    int t=0,s=0;
    for(int k=last[x];k;k=a[k].next)
    {
        int y=a[k].y;
        if(a[k].c>0 && h[y]==h[x]+1 && s<f)
        {
            s+=(t=findflow(y,min(a[k].c,f-s)));
            a[k].c-=t; a[a[k].other].c+=t;
        }
    }
    if(s==0)h[x]=0;
    return s;
}
int main()
{
    scanf("%d%d",&n,&m);
    memset(last,0,sizeof(last));len=0;
    for(int i=1;i<=m;i++)
    {
        int x,y,c,d;
        scanf("%d%d%d%d",&x,&y,&c,&d);
        ins(x,y,c,d);ins(y,x,c,d);
    }
    memset(v,false,sizeof(v));
    for (int i=1;i<=n;i++) d[i]=999999999;
    head=1;tail=2;
    list[1]=1;v[1]=true;d[1]=0;
    while(head!=tail)
    {
        int x=list[head];
        for(int k=last[x];k;k=a[k].next)
        {
            int y=a[k].y;
            if(d[y]>d[x]+a[k].c)
            {
                d[y]=d[x]+a[k].c;
                if(v[y]==false)
                {
                    v[y]=true;
                    list[tail]=y;
                    tail++; if(tail==n+1) tail=1;
                }
            }
        }
        head++; if(head==n+1) head=1;
        v[x]=false;
    }
    printf("%d\n",d[n]);
    int nd=0;
    for(int i=1;i<=len;i++)if(d[a[i].x]+a[i].c==d[a[i].y])addx[++nd]=a[i].x,addy[nd]=a[i].y,addd[nd]=a[i].d;
    len=0;memset(last,0,sizeof(last));
    for(int i=1;i<=nd;i++)insert(i);
    st=1;ed=n;int ans=0;
    while(bt_h())ans+=findflow(st,999999999);
    printf("%d\n",ans);
    return 0;
}