UVA - 125 Numbering Paths floyd的变形

题目大意：给出N个城市，要求求出这每个城市之间有多少条道路可走

解题思路：floyd算法的思想是DP思想，用f[i][j]表示i到j的值，这个值可以是多种多样的，比如前面的UVA 10048的噪音题f[i][j]就表示从i到j通过的点的最大值的最小值，通过中间节点k来过度，这题的话f[i][j]就表示从i到j有多少条线路可走，转换方程是，当f[i][k] && f[k][j]时，也就是i到j时通过节点k可通，则f[i][j] += f[i][k] * f[k][j]，当f[k][k]不等于0时，就表示在k点处有一个环了，再通过判断哪几个点的连接时需要通过k时，即可判断路径数是否为无限

#include<cstdio>
#include<cstring>
#define maxn 35
int f[maxn][maxn];
int main() {
	int mark = 0,num,t1,t2;
	while(scanf("%d",&num) != EOF && num) {
		int N = 0;	
		memset(f,0,sizeof(f));
		for(int i = 0; i < num; i++) {
			scanf("%d%d",&t1,&t2);	
			f[t1][t2] = 1;
			if(t1 > N)
				N = t1;
			if(t2 > N)
				N = t2;
		}

		for(int k = 0; k <= N; k++)
			for(int i = 0; i <= N; i++)
				for(int j = 0 ; j <= N; j++)
					if(f[i][k] && f[k][j])
						f[i][j] += f[i][k] * f[k][j];
			
		for(int k = 0; k <= N; k++)
			if(f[k][k]) 
				for(int i = 0; i <= N; i++)
					for(int j = 0; j <= N; j++)
						if(f[i][k] && f[k][j])
							f[i][j] = -1;
		printf("matrix for city %d\n",mark++);
		for(int i = 0; i <= N; i++) {
			for(int j = 0; j <= N; j++) {
				if(j == 0)
					printf("%d",f[i][j]);
				else
					printf(" %d",f[i][j]);	
			}
			printf("\n");
		}	
	}	
	return 0;
}