bzoj 4870: [Shoi2017]组合数问题

题意：

求题面上那个乱七八糟的柿子。

题解：

其实是道大傻逼题
千万不要像我一样被那个玩意吓坏了。
地球版题意：在n*k个物品中选若干个物品，使其数量$mod\ p=r$。
然后dp+矩乘。
我真是大傻逼
code：

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#define LL long long
using namespace std;
struct node{
    LL n,m,a[55][55];
}em,f;
LL n,mod,k,r;
node operator * (node a,node b)
{
    node ans;ans=em;
    ans.n=a.n;ans.m=b.m;
    for(LL i=0;i<a.n;i++)
        for(LL j=0;j<b.m;j++)
            for(LL k=0;k<a.m;k++)
                (ans.a[i][j]+=a.a[i][k]*b.a[k][j])%=mod;
    return ans;
}
node operator ^ (node a,LL b)
{
    node ans;ans=em;
    ans.n=a.n;ans.m=a.m;
    for(LL i=0;i<a.n;i++) ans.a[i][i]=1;
    while(b)
    {
        if(b&1) ans=ans*a;
        a=a*a;b>>=1;
    }
    return ans;
}
int main()
{
    scanf("%lld %lld %lld %lld",&n,&mod,&k,&r);
    f.n=1;f.m=k;f.a[0][0]=1;
    node a;a.n=a.m=k;
    for(LL i=0;i<k;i++)
        a.a[i][i]++,a.a[(i-1+k)%k][i]++;
    f=f*(a^(n*k));
    printf("%lld",f.a[0][r]);
}