华南师范大学·2012年ACM程序设计竞赛·决赛 1011

构建概率矩阵，然后矩阵快速幂，直接输出答案。

#include<cstdio>
#include<vector>
#include<cstring>
using namespace std;
int n,m,k;
vector<int> g[201];
struct matrix
{
    float m[201][201];
    matrix(){memset(m,0,sizeof(m));}
    matrix operator *(const matrix &a)const
    {
        matrix tmp;
        for(int i=0;i<n;i++)
            for(int j=0;j<n;j++)
              for(int k=0;k<n;k++)
                tmp.m[i][j]+=m[i][k]*a.m[k][j];
        return tmp;
    }
};
matrix posibility;
matrix pow(matrix a,int n)
{
    if(n==1) return a;
    if(n%2)
    {
        matrix tmp=pow(a,n/2);
        return tmp*tmp*a;
    }
    else
    {
        matrix tmp=pow(a,n/2);
        return tmp*tmp;
    }
}
int main()
{
    int t;
    scanf("%d",&t);
    for(int CASE=1;CASE<=t;CASE++)
    {
        scanf("%d%d%d",&n,&m,&k);
        for(int i=0;i<n;i++)
            g[i].clear();
        while(m--)
        {
            int u,v;
            scanf("%d%d",&u,&v);
            u--;v--;
            g[u].push_back(v);
            g[v].push_back(u);
        }
        memset(posibility.m,0,sizeof(posibility.m));
        for(int i=0;i<n;i++)
        {
           double cnt=g[i].size();
           for(int j=0;j<cnt;j++)
           {
               int v=g[i][j];
               posibility.m[i][v]=1/cnt;
           }
        }
        matrix ans=pow(posibility,k);
        printf("Case #%d: %.4f
",CASE,ans.m[0][1]);
    }
    return 0;
}