*寒假水87——How many ways?? 【矩阵快速幂】（不懂）-优快云博客

本文介绍了一种利用矩阵快速幂算法解决特定图论问题的方法，即计算从起点A到终点B恰好经过k个点的不同路径数量。通过具体示例展示了如何构建矩阵并使用快速幂运算来高效求解该问题。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

春天到了, HDU校园里开满了花, 姹紫嫣红, 非常美丽. 葱头是个爱花的人, 看着校花校草竞相开放, 漫步校园, 心情也变得舒畅. 为了多看看这迷人的校园, 葱头决定, 每次上课都走不同的路线去教室, 但是由于时间问题, 每次只能经过k个地方, 比方说, 这次葱头决定经过2个地方, 那他可以先去问鼎广场看看喷泉, 再去教室, 也可以先到体育场跑几圈, 再到教室. 他非常想知道, 从A 点恰好经过k个点到达B点的方案数, 当然这个数有可能非常大, 所以你只要输出它模上1000的余数就可以了. 你能帮帮他么?? 你可决定了葱头一天能看多少校花哦

Input输入数据有多组, 每组的第一行是2个整数 n, m(0 < n <= 20, m <= 100) 表示校园内共有n个点, 为了方便起见, 点从0到n-1编号,接着有m行, 每行有两个整数 s, t (0<=s,t<n) 表示从s点能到t点, 注意图是有向的.接着的一行是两个整数T,表示有T组询问(1<=T<=100),
接下来的T行, 每行有三个整数 A, B, k, 表示问你从A 点到 B点恰好经过k个点的方案数 (k < 20), 可以走重复边。如果不存在这样的走法, 则输出0
当n, m都为0的时候输入结束
Output计算每次询问的方案数, 由于走法很多, 输出其对1000取模的结果Sample Input

Sample Output

#include<iostream>
#include<stdio.h>
#include<string.h>

using namespace std;

int i,j,k,n,m,mod=1000,q;

struct node
{
    int a[25][25];
    friend node operator *(node a,node b)
    {
        node c;
        for(i=0;i<n;i++)
        {
            for(j=0;j<n;j++)
            {
                c.a[i][j] = 0;
                for(k=0;k<n;k++)
                    c.a[i][j] += a.a[i][k]*b.a[k][j],c.a[i][j]%=mod;
            }
        }
        return c;
    }
}A;

node qkm(node A,int m)
{
	int i,j;
    node ans,base = A;
    for(i=0;i<n;i++)
    {
        for(j=0;j<n;j++)
        {
            if(i==j)
                ans.a[i][j]=1;
            else
                ans.a[i][j]=0;
        }
    }
    while(m)
    {
        if(m%2)
            ans=ans*base;
        base=base*base;
        m/=2;
    }
    return ans;
}

int main()
{
	int a,b,k,i;
    while(scanf("%d%d",&n,&m),n+m)
    {
        memset(A.a,0,sizeof A.a);
        for(i=0;i<m;i++)
        {
            int s,t;
            scanf("%d%d",&s,&t);
            A.a[s][t]=1;
        }
        scanf("%d",&q);
        while(q--)
        {
            
            scanf("%d%d%d",&a,&b,&k);
            node C;
            C = qkm(A,k);
            printf("%d\n",C.a[a][b]);
        }
    }
    return 0;
}

题解：关于矩阵快速幂点击打开链接