ZOJ 2760 How Many Shortest Path（floyd + 最大流）

思路：

首先用floyd跑出最短路，然后对于每条边判断，dis[j]==dis[i]+<i,j>
满足条件即加到网络中，（i，j，1），求解最大流，即为答案。

坑点：

注意他给的样例有诈。。
当i等于j的时候需要自己手动置0，而不是指望他给你0。。
靠。因为这个狂wa不止。。

#include <iostream>#include <cstdio>#include <queue>#include <string.h>#define inf 0x3f3f3f3ftypedef long long int lli;using namespace std;const int maxn = 5200;struct edge{    lli to,next,v;}ed[600000+20000];lli ans;lli head[maxn],dis[maxn],cnte,cur[maxn];void ini(){    memset(head,-1,sizeof(head));    ans = 0;cnte = 0;}inline void ae(int x,int y,int v){    ed[cnte].to = y;    ed[cnte].v = v;    ed[cnte].next=head[x];    head[x]=cnte++;    ed[cnte].to = x;    ed[cnte].v = 0;    ed[cnte].next=head[y];    head[y]=cnte++;}int bfs(int s,int t){    memset(dis,-1,sizeof(dis));    queue<int> q;q.push(s);    dis[s] = 0;    int now;    while(!q.empty()){        now = q.front();q.pop();        for(int i = head[now];i!=-1;i=ed[i].next){            int v = ed[i].to;            if(dis[v] == -1 && ed[i].v > 0){                dis[v] = dis[now]+1;q.push(v);                if(v == t) return 1;            }        }    }    if(dis[t] < 0) return 0;    return 1;}lli dfs(int s,int t,lli low){    if(s==t || !low) return low;    lli a,ret = 0;    for(lli &i = cur[s];i!=-1;i=ed[i].next){        int v= ed[i].to;        if(dis[v]==dis[s]+1 && ed[i].v>0 && (a=dfs(v,t,min(low,ed[i].v)))){            ed[i].v -= a; ed[i^1].v += a;            low -= a; ret += a;            if(!low) break;        }    }    return ret;}void dinic(int s,int t){    while(bfs(s,t)){        memcpy(cur,head,sizeof(head));        ans += dfs(s,t,inf);    }}int ma[110][110];int ma2[110][110];void floyd(int n){    for(int k = 1;k <= n;k++){        for(int i = 1;i <= n;i++){            for(int j = 1;j <= n;j++){                ma2[i][j] = min(ma2[i][j],ma2[i][k]+ma2[k][j]);            }        }    }}int ddis[maxn];int main(){    int cas,n;    while(~scanf("%d",&n)){        int temp,s,t;        ini();        memset(ma,0x3f,sizeof(ma));        for(int i = 1;i <= n;i++){            for(int j = 1;j <= n;j++){                scanf("%d",&temp);                if(temp != -1)                    ma[i][j] = temp;                if(i == j)                    ma[i][j] = 0;            }        }        scanf("%d%d",&s,&t);        if(s==t){            puts("inf");continue;        }        memcpy(ma2,ma,sizeof(ma));        floyd(n);        for(int i = 1;i <= n;i++){            ddis[i] = ma2[s+1][i];        }        for(int i = 1;i <= n;i++){            if(ddis[i] == inf) continue;            for(int j = 1;j <= n;j++){                if(ma[i][j] != inf &&i!=j&& ddis[j] == ddis[i]+ma[i][j]){                    ae(i,j,1);                }            }        }        dinic(s+1,t+1);        printf("%lld\n",ans);    }    return 0;}