hdu 5001 矩阵 快速幂

题目大意：给n个点，m条边，每个点选择每条边都是任意的，求走K步，不经过每一个点的概率。

这里用了一个struct代表一种类型（矩阵类型），这样写比较方便，以后可以这样用！！！

一开始我题目理解错了，写错了也就算了，可是当题目理解对了之后，还是错，真的是。。。虽然我知道我的代码能力差，但是总是犯一些很傻比的错误，我真的很忧伤，所以呢，以后写题要更加小心才可以，也要完完全全想清楚才行！！！！点击打开链接

#include <cstdio>#include <cstring>#include <iostream>#include <algorithm>using namespace std;int node[60], n, m, step;bool flag[60][60];struct matrix{    double a[60][60];};matrix rela, ans, temp;matrix mut(matrix A, matrix B){    int i, j;    matrix res;    for(i = 1; i <= n; i++)    {        for(j = 1; j <= n; j++)        {            res.a[i][j] = 0;            for(int t = 1; t <= n; t++)                res.a[i][j] += A.a[i][t] * B.a[t][j];        }    }    return res;}void solve( int step ){    int i, j;    for(i = 1; i <= n; i++)        for(j = 1; j <= n; j++)            ans.a[i][j] = (i == j);    while(step)    {        if(step & 1)           ans = mut(ans, rela);        rela = mut(rela, rela);        step /= 2;    }}int main(){    int T, cas, i, j;    int pp, qq;    scanf("%d", &T);    for(cas = 1; cas <= T; cas++)    {        scanf("%d%d%d", &n, &m, &step);        memset(node, 0, sizeof(node));        memset(flag, 0, sizeof(flag));        for(i = 1; i <= m; i++)        {            scanf("%d%d", &pp, &qq);            node[pp]++;   node[qq]++;            flag[pp][qq] = flag[qq][pp] = 1;        }        for(i = 1; i <= n; i++)        {            for(j = 1; j <= n; j++)            {                rela.a[j][j] = 0;                for(int t = 1; t <= n; t++)                {                    if(t!=j)                    {                        if(j==i || t==i)                            rela.a[j][t] = 0;                        else if(flag[j][t])                            rela.a[j][t] = 1.0/node[j];                        else                            rela.a[j][t] = 0;                    }                }            }            //一开始写得是solve（），所以导致step在第一次计算后就变了，所以后面的只不对、            solve( step );            double mm;            double sum = 0;            for ( int j = 1 ; j <= n ; j ++ ) {                for ( int k = 1 ; k <= n; k ++)                    sum += ans.a[j][k] ;            }            printf("%.8lf\n", sum/n);        }    }    return 0;}