给定有向图，每条边只能用一次，求给定两点间有几条最短路。

1. 先用flyd求出所有最短路，之后用最短路建图就网络流。
2. 重点是如何利用最短路建图！！
3. 
   
4. 
   
5. 
   
6. #include<cstdio>  
7. #include<cstring>  
8. #include<queue>  
9. #define find_min(a,b) a<b?a:b  
10. using namespace std;  
11.   
12. const int N = 210;  
13. const int MAX = 9999999999;  
14.   
15. struct Edge{  
16.     int s,e,v,next;  
17. }edge[5*N];  
18.   
19. int e_num,head[N],d[N],sp,tp;  
20. int n,dist[N][N],mat[N][N];  
21.   
22. void AddEdge(int a,int b,int c){  
23.     edge[e_num].s=a; edge[e_num].e=b; edge[e_num].v=c;  
24.     edge[e_num].next=head[a]; head[a]=e_num++;  
25.   
26.     edge[e_num].s=b; edge[e_num].e=a; edge[e_num].v=0;  
27.     edge[e_num].next=head[b]; head[b]=e_num++;  
28. }  
29.   
30. int bfs(){  
31.     queue <int> q;  
32.     memset(d,-1,sizeof(d));  
33.     d[sp]=0;  
34.     q.push(sp);  
35.     while(!q.empty()){  
36.         int cur=q.front();  
37.         q.pop();  
38.         for(int i=head[cur];i!=-1;i=edge[i].next){  
39.             int u=edge[i].e;  
40.             if(d[u]==-1 && edge[i].v>0){//没有标记,且可行流大于0  
41.                 d[u]=d[cur]+1;  
42.                 q.push(u);  
43.             }  
44.         }  
45.     }  
46.     return d[tp] != -1;//汇点是否成功标号,也就是说是否找到增广路  
47. }  
48.   
49. int dfs(int a,int b){//a为起点  
50.     int r=0;  
51.     if(a==tp)return b;  
52.     for(int i=head[a];i!=-1 && r<b;i=edge[i].next){  
53.         int u=edge[i].e;  
54.         if(edge[i].v>0 && d[u]==d[a]+1){  
55.             int x=find_min(edge[i].v,b-r);  
56.             x=dfs(u,x);  
57.             r+=x;  
58.             edge[i].v-=x;  
59.             edge[i^1].v+=x;  
60.         }  
61.     }  
62.     if(!r)d[a]=-2;  
63.     return r;  
64. }  
65.   
66. void dinic(int sp,int tp){  
67.     int total=0,t;  
68.     while(bfs()){  
69.         while(t=dfs(sp,MAX))  
70.         total+=t;  
71.     }  
72.     printf("%d\n",total);  
73. }  
74. int main()  
75. {  
76.     int i,j,k;  
77.     while(~scanf("%d",&n)){      
78.         for(i=0;i<n;i++){  
79.             for(j=0;j<n;j++){  
80.                 scanf("%d",&mat[i][j]);  
81.                 if(i==j)mat[i][j]=0;  
82.                 if(mat[i][j]==-1)mat[i][j]=MAX;  
83.                 dist[i][j]=mat[i][j];  
84.             }  
85.         }scanf("%d%d",&sp,&tp);  
86.           
87.         if(sp!=tp){  
88.             for(k=0;k<n;k++){  
89.                 for(i=0;i<n;i++)  
90.                     if(dist[i][k]<MAX)for(j=0;j<n;j++)  
91.                         if(dist[k][j]<MAX)dist[i][j]=find_min(dist[i][j],dist[i][k]+dist[k][j]);  
92.             }  
93.               
94.             e_num=0;  
95.             memset(head,-1,sizeof(head));  
96.             for(i=0;i<n;i++){  
97.                 if(dist[sp][i]<MAX)for(j=0;j<n;j++)  
98.                     if(dist[j][tp]<MAX && mat[i][j]<MAX && dist[sp][i]+mat[i][j]+dist[j][tp]==dist[sp][tp])AddEdge(i,j,1);  
99.             }  
100.   
101.             dinic(sp,tp);  
102.         }  
103.         else puts("inf");  
104.     }      
105.     return 0;      
106. }