题目链接:https://vjudge.net/problem/UVALive-3704
具体思路:用一维数组表示二维数组。
AC代码:
#include
#include
#include
#include
#include
#include
#include
#include
#include
using namespace std;
# define inf 0x3f3f3f3f
# define maxn 500+100
# define ll long long
ll n,m,d,k;
ll ans[maxn];
ll fir[maxn];
struct Matrix
{
ll a[maxn];
} t1,t2;
Matrix mul(Matrix w1, Matrix w2)
{
Matrix temp;
for(ll j=1; j<=n; j++)
{
temp.a[j]=0;
for(ll k=1; k<=n; k++)
{
temp.a[j]=(temp.a[j]+w1.a[k]*w2.a[(j-k+n)%n+1])%m;
}
}
return temp;
}
Matrix quickpow() // t1 - > n*n , t2 - > 1*n
{
Matrix tt=t2;
k--;
while(k)
{
if(k&1)
{
tt=mul(tt,t1);
}
t1=mul(t1,t1);
k=k>>1;
}
return tt;
}
int main()
{
while(~scanf("%lld%lld%lld%lld",&n,&m,&d,&k))
{
for(ll i=1; i<=n; i++)
{
scanf("%lld",&fir[i]);
}
for(ll i=1; i<=n; i++)
{
t1.a[i]=0;
}
for(ll i=1; i<=d+1; i++)
{
t1.a[i]=1;
}
for(ll i=1; i<=d; i++)
{
t1.a[n-i+1]=1;
}
t2=t1;
Matrix w=quickpow();
memset(ans,0,sizeof(ans));
for(ll i=1; i<=n; i++)
{
for(ll j=1; j<=n; j++)
{
ans[i]=(ans[i]+w.a[(i-j+n)%n+1]*fir[j])%m;
}
}
for(ll i=1; i<=n; i++)
{
printf("%lld",ans[i]);
if(i!=n)printf(" ");
}
printf("\n");
}
return 0;
}