洛谷 P3746 [六省联考2017]组合数问题矩阵乘法

最新推荐文章于 2022-03-20 20:08:16 发布

Amber_lylovely

最新推荐文章于 2022-03-20 20:08:16 发布

阅读量246

点赞数

分类专栏：矩阵乘法

本文链接：https://blog.csdn.net/liangzihao1/article/details/89446956

版权

矩阵乘法专栏收录该内容

8 篇文章 0 订阅

订阅专栏

该博客探讨了洛谷P3746题目，涉及组合数问题的解决策略，重点在于利用矩阵乘法进行状态转移。博主详细分析了如何通过递推公式f[i][j]=f[i−1][j]+f[i−1][(j+k−1) mod k]来计算在m mod k=j条件下的(mi)之和，并指出由于k较小，可以采用矩阵快速幂优化算法提高效率。文章包含问题链接和解决方案代码。

摘要由CSDN通过智能技术生成

题目：
https://www.luogu.org/problemnew/show/P3746

分析：
设 $f [i] [j]$ 表示在 $m\ mod\ k=j$ 的 $\binom{i}{m}$ 的和。
转移显然， $f[i][j]=f[i-1][j]+f[i-1][(j+k-1)\ mod\ k]$
因为 $k$ 较小，可以构造出矩阵进行递推，其中 $f [0] [0] = 1$ 。

代码：

// luogu-judger-enable-o2
#include <iostream>
#include <cstdio>
#include <cmath>
#define LL long long

const int N=57;

using namespace std;

LL n,mod;
int k,r;

struct matrix{
    int n,m;
    LL a[N][N];
}A,B,C;

matrix operator *(matrix a,matrix b)
{
    matrix c;
    c.n=a.n,c.m=b.m;
    for (int i=1;i<=c.n;i++)
    {
        for (int j=1;j<=c.m;j++) c.a[i][j]=0;
    }
    for (int k=1;k<=a.m;k++)
    {
        for (int i=1;i<=c.n;i++)
        {
            for (int j=1;j<=c.m;j++)
            {
                c.a[i][j]=(c.a[i][j]+a.a[i][k]*b.a[k][j]%mod)%mod;
            }
        }
    }
    return c;
}

void ksm(LL p)
{
    if (p==0) return;
    ksm(p/2);
    B=B*B;
    if (p&1) B=B*A;
}

int main()
{
    scanf("%lld%lld%d%d",&n,&mod,&k,&r);
    A.n=A.m=B.n=B.m=k;
    for (int i=1;i<=k;i++)
    {
        A.a[i][i]++;
        if (i==1) A.a[i][k]++;
             else A.a[i][i-1]++;
        B.a[i][i]=1;
    }	
    ksm(n*(LL)k);
    C.n=1,C.m=k;
    C.a[1][1]=1;
    C=C*B;	
    printf("%lld\n",C.a[1][r+1]);
}