syhd142  
日历
<2010年9月>
2930311234
567891011
12131415161718
19202122232425
262728293012
3456789
统计
  • 随笔 - 23
  • 文章 - 122
  • 评论 - 31
  • 引用 - 0

导航

常用链接

留言簿(2)

随笔档案(23)

文章分类(270)

文章档案(122)

我的豆瓣

搜索

  •  

最新评论

阅读排行榜

评论排行榜

 
传递闭包,矩阵相乘。
#include <stdio.h>
#include 
<string.h>

#define N 105

struct Matrix
{
    
int mat[N][N];
};

void print(Matrix g, int n)
{
    
for(int i = 1; i <= n; i++)
    {
        
for(int j = 1; j <= n; j++)
        {
            printf(
"%3d", g.mat[i][j]);
        }
        printf(
"\n");
    }
}

Matrix MatrixMul(Matrix t, Matrix g, 
int n)
{
    Matrix ans;
    
for(int i = 1; i <= n; i++)
    {
        
for(int j = 1; j <= n; j++)
        {
            
int sum = 0;
            
for(int k = 1; k <= n; k++)
            {
                sum 
|= t.mat[i][k] & g.mat[k][j];
            }
            ans.mat[i][j] 
= sum;
        }
    }
    
return ans;
}

Matrix MatrixPow(Matrix g, 
int n, int d)
{
    Matrix t;
    
for(int i = 1; i <= n; i++)
    {
        
for(int j = 1; j <= n; j++)
        {
            
if(i == j) t.mat[i][j] = 1;
            
else t.mat[i][j] = 0;
        }
    }
//    print(t, n);
    while(d)
    {
        
if(d & 1) t = MatrixMul(t, g, n);
        g 
= MatrixMul(g, g, n);
        d 
>>= 1;
    }
    
return t;
}

int main()
{
    
//freopen("in", "r", stdin);
    int n, m, s, e, d;
    Matrix g;
    
while(scanf("%d %d"&n, &m), n + m)
    {
        memset(g.mat, 
0sizeof(g.mat));
        
for(int i = 0; i < m; i++)
        {
            scanf(
"%d %d"&s, &e);
            g.mat[s][e] 
= g.mat[e][s] = 1;
        }
        scanf(
"%d %d %d"&s, &e, &d);
        g 
= MatrixPow(g, n, d);
        
if(g.mat[s][e]) puts("Yes, Teobaldo can travel.");
        
else puts("No, Teobaldo can not travel.");
    }
    
return 0;
}
posted on 2010-09-12 10:52 Fucker 阅读(303) 评论(0)  编辑 收藏 引用 所属分类: ACM/ICPC图论

只有注册用户登录后才能发表评论。
网站导航: 博客园   IT新闻   BlogJava   博问   Chat2DB   管理


 
Copyright © Fucker Powered by: 博客园 模板提供:沪江博客