思路：

首先 先普及一個性質： 循環矩陣*循環矩陣=循環矩陣
由于此題是距離小于d的都加上一個數。
那么 構造矩陣的時候 我們發現 誒呦 這是個循環矩陣
看看數據范圍 n^2log(k)可以過。
那就把這個矩陣改一改。
因為這是個循環矩陣， 所以呢 只用保存一行就可以了。
每回做乘法的時候只做第一行的乘法。
for(i) for(j) temp[i]+=a[j]*b[(i+j)%n];
就這么著 搞搞就能過了。

（好像可以用FFT？ 表示并不會）

// by SiriusRen #include <cstdio> #include <cstring> #define LL long long using namespace std; long long f[666],co[666],c[666]; int n,m,d,k; void mul(LL *a,LL *b){memset(c,0,sizeof(c));for(int i=0;i<n;i++)for(int j=0;j<n;j++)c[i]+=a[j]*b[(i+j)%n];for(int i=0;i<n;i++)b[i]=c[i]%m; } int main(){scanf("%d%d%d%d",&n,&m,&d,&k);for(int i=0;i<n;i++)scanf("%lld",&f[i]);for(int i=0;i<=d;i++)co[i]=co[n-i]=1;while(k){if(k&1)mul(co,f);mul(co,co),k>>=1;}for(int i=0;i<n;i++)printf("%lld ",f[i]); }

// by SiriusRen #include <cstdio> #define LL long long #define f(Q,R) for(int Q=0;Q<R;Q++) using namespace std; LL f[666],co[666],c[666]; int n,m,d,k; void mul(LL *a,LL *b){f(i,n)c[i]=0;f(i,n)f(j,n)c[(i+j)%n]+=a[i]*b[j];f(i,n)b[i]=c[i]%m; } int main(){scanf("%d%d%d%d",&n,&m,&d,&k);f(i,n)scanf("%lld",&f[i]);f(i,d+1)co[i]=co[n-i]=1;while(k){if(k&1)mul(co,f);mul(co,co),k>>=1;}f(i,n)printf("%lld ",f[i]); }

Code length能進前三的存在哈哈哈

轉載于:https://www.cnblogs.com/SiriusRen/p/6532384.html

總結

以上是生活随笔為你收集整理的POJ 3150 循环矩阵的应用的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。