bzoj1266最短路+最小割

本来写了spfa

wa了

看到网上有人写Floyd过了

表示不开心￣へ￣

改成Floyd试试。。。

还是wa

ヾ(｡｀Д´｡)原来是建图错了（样例怎么过的）

结果T了

于是把Floyd改回spfa

还是T了

。。。

￣へ￣

看来问题不在最短路，改回Floyd（mdzz）

。。。

好像dinic有点问题

( ⊙ o ⊙ )A了

#include <cstdio>  
#include <cstring>  
#include <iostream>  
#include <algorithm>  
#define M 510  
#define S 1  
#define T n  
#define INF 1000000000 
using namespace std;  
struct edge{  
    int x,y,f,c;  
}e[124800];  
int n,m,ans,i,j,k,x,y,z;  
int map[M][M];
int head[M],tot=1;
int dpt[M]; 
struct abcd{  
    int to,next,f;  
}table[124800<<2];
void add(int x,int y,int z)  
{  
    table[++tot].to=y;  
    table[tot].f=z;  
    table[tot].next=head[x];  
    head[x]=tot;  
    table[++tot].to=x;  
    table[tot].f=0;  
    table[tot].next=head[y];  
    head[y]=tot; 
}  
    bool bfs()  
    {  
        static int q[M];  
        int i,r=0,h=0;  
        memset(dpt,-1,sizeof dpt);  
        dpt[S]=1;q[++r]=S;  
        while(r!=h)  
        {  
            int x=q[++h];  
            for(i=head[x];i;i=table[i].next)  
                if(table[i].f&&!~dpt[table[i].to])  
                {  
                    dpt[table[i].to]=dpt[x]+1;  
                    q[++r]=table[i].to;  
                    if(table[i].to==T)  
                        return true;  
                }  
        }  
        return false;  
    }  
    int dfs(int x,int flow)  
    {  
        int i,left=flow;  
        if(x==T) return flow;  
        for(i=head[x];i&&left;i=table[i].next)  
            if(table[i].f&&dpt[table[i].to]==dpt[x]+1)  
            {  
                int temp=dfs(table[i].to,min(left,table[i].f) );  
                left-=temp;  
                table[i].f-=temp;  
                table[i^1].f+=temp;  
            }  
        if(left) dpt[x]=-1;  
     return flow-left;  
    }  
int main()  
{  
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
            map[i][j]=INF;
    for(i=1;i<=n;i++)  
        map[i][i]=0;  
    for(i=1;i<=m;i++)  
    {  
        scanf("%d%d%d%d",&x,&y,&z,&e[i].c);  
        e[i].x=x;e[i].y=y;e[i].f=z;  
        map[x][y]=map[y][x]=z;  
    }  
    for(k=1;k<=n;k++)  
        for(i=1;i<=n;i++)  
            for(j=1;j<=n;j++)  
                map[i][j]=min(map[i][j],map[i][k]+map[k][j]);  
    printf("%d
",map[1][n]);
    for(i=1;i<=m;i++)  
    {  
        x=e[i].x;y=e[i].y;z=e[i].f;  
        if( map[1][x]+map[y][n]+z==map[1][n] )  
            add(x,y,e[i].c);  
        if( map[1][y]+map[x][n]+z==map[1][n] )  
            add(y,x,e[i].c);  
    }  
    while( bfs() )  
        ans+=dfs(1,INF);  
    printf("%d",ans);
    return 0;  
}