洛谷3746

一开始没去想正解，对着部分分口胡了75分。
然后想正解，跟常人不同的是，我并没有把考虑k，接着就死活想不出，我真是太菜了。

#include<cstdio>#include<cstring>const int N=55;int n,p,k,r,i;long long m;struct matrix{    int a[N][N];    inline void clear(){        memset(a,0,sizeof a);    }    matrix operator*(const matrix&x){        matrix ans;ans.clear();        register int i,j,l;        for(i=0;i<k;++i)            for(l=0;l<k;++l)                for(j=0;j<k;++j)ans.a[i][j]=(ans.a[i][j]+1ll*a[i][l]*x.a[l][j])%p;        return ans;    }}a,b;int main(){    scanf("%d%d%d%d",&n,&p,&k,&r);    for(i=0;i<k;++i)a.a[i][i]=1,++b.a[i][i],++b.a[i][(i+1)%k];    m=1ll*n*k;    for(;m;m>>=1,b=b*b)if(m&1)a=a*b;    printf("%d\n",a.a[r][0]);    return 0;}