1 void floyd() 2 { 3 int i,j,k; 4 for(k=1;k<=n;k++) 5 for(i=1;i<=n;i++) 6 for(j=1;j<=n;j++) 7 { 8 if(dp[i][j]>dp[i][k]+dp[k][j]) 9 dp[i][j]=dp[i][k]+dp[k][j]; 10 } 11 }
显然我们可以从代码里看出,该算法的时间复杂度为O()。
1 void floyd() 2 { 3 int i,j,k; 4 for(k=1;k<=n;k++) 5 for(i=1;i<=n;i++) 6 for(j=1;j<=n;j++) 7 if(matrix[i][k]&&matrix[k][j]) 8 matrix[i][j]=1; 9 }
若有会Warshall算法求传递闭包的朋友便会发现Floyd算法和Warshall算法高度相似!
相关题目:ZOJ P41241 #include <bits/stdc++.h> 2 3 using namespace std; 4 5 int n,m,matrix[105][105],num[2][105],flag=0; 6 7 void floyd() 8 { 9 int i,j,k; 10 11 for(k=1;k<=n;k++) //floyd 12 for(i=1;i<=n;i++) 13 for(j=1;j<=n;j++) 14 if(matrix[i][k]&&matrix[k][j]) 15 matrix[i][j]=1; 16 17 for(i=1;i<=n;i++) //判断是否有自环 18 for(j=1;j<=n;j++) 19 if(matrix[i][j]&&matrix[j][i]) 20 { 21 flag=1; 22 return; 23 } 24 25 for(i=1;i<=n;i++) //维护更新num数组 26 for(j=1;j<=n;j++) 27 if(matrix[i][j]) 28 { 29 num[0][i]++; 30 num[1][j]++; 31 } 32 } 33 34 int main() 35 { 36 int t; 37 scanf("%d",&t); 38 while(t--) 39 { 40 int i,a,b; 41 42 memset(matrix,0,sizeof(matrix)); //初始化 43 memset(num,0,sizeof(num)); 44 flag=0; 45 46 scanf("%d%d",&n,&m); 47 for(i=0;i<m;i++) //预处理 48 { 49 scanf("%d%d",&a,&b); 50 matrix[a][b]=1; 51 } 52 floyd(); 53 if(flag==0) //输出,当一个数比它大的和比它小的都小于n/2时,可视为该项为中间项 54 { 55 for(i=1;i<=n;i++) 56 { 57 if(num[0][i]<=n/2&&num[1][i]<=n/2) printf("1"); 58 else printf("0"); 59 } 60 printf("\n"); 61 } 62 else 63 { 64 for(i=0;i<n;i++) printf("0"); 65 printf("\n"); 66 } 67 } 68 return 0; 69 }
(用Floyd算法解此题算是一个比较巧妙的解法)
四、相关题目
1.例题 Luogu P2910
代码:
1 #include <bits/stdc++.h> 2 3 using namespace std; 4 5 int matrix[105][105],n; 6 7 void floyd() 8 { 9 int i,j,k; 10 for(k=1;k<=n;k++) 11 for(i=1;i<=n;i++) 12 for(j=1;j<=n;j++) 13 { 14 matrix[i][j]=min(matrix[i][j],matrix[i][k]+matrix[k][j]); 15 } 16 } 17 18 19 int main() 20 { 21 int i,j,m,order[10005]={0},ans=0; 22 scanf("%d%d",&n,&m); 23 for(i=0;i<m;i++) scanf("%d",&order[i]); 24 for(i=1;i<=n;i++) 25 for(j=1;j<=n;j++) scanf("%d",&matrix[i][j]); 26 floyd(); 27 for(i=0;i<m-1;i++) 28 ans+=matrix[order[i]][order[i+1]]; 29 printf("%d",ans); 30 return 0; 31 }
Author : Houge Date : 2019.6.1
Update log :
原文:https://www.cnblogs.com/CSGOBESTGAMEEVER/p/10936333.html