BZOJ1712: [Usaco2007 China]Summing Sums 加密

n<=50000个数做t<=1?????????反正不爆int次操作，每次把一个数组中每一个数变成当前，当前，当前数组的其他数的和，就是然后求最后每个Ai

每个数做的操作是一样的而且每个数加上的其他数的系数都是一样的，实验可得最后每个数都会变成：

至于x，y用矩阵算，base矩阵如下：

#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#include<algorithm>
//#include<iostream>
using namespace std;

#define LL long long
typedef LL mat[4][4];
mat start,a,t;
const int mod=98765431;
void init(mat &a)
{
    memset(a,0,sizeof(a));
    a[1][1]=a[2][2]=1;
}
void copy(mat &a,mat b) {memcpy(a,b,sizeof(mat));}
void mul(mat a,mat b,mat &ans)
{
    mat t;memset(t,0,sizeof(t));
    for (int i=1;i<=2;i++)
        for (int j=1;j<=2;j++)
            for (int k=1;k<=2;k++)
                t[i][j]=(t[i][j]+a[i][k]*b[k][j]%mod)%mod;
    copy(ans,t);
}
void pow(mat a,LL b,mat &ans)
{
    mat t,tmp;init(t);copy(tmp,a);
    while (b)
    {
        if (b&1) mul(t,tmp,t);
        mul(tmp,tmp,tmp);
        b>>=1;
    }
    copy(ans,t);
}
int n,K;
#define maxn 50011
LL num[maxn],sum;
int main()
{
    scanf("%d%d",&n,&K);sum=0;
    for (int i=1;i<=n;i++) scanf("%lld",&num[i]),sum+=num[i],sum-=sum>mod?mod:0;
    start[1][1]=1;start[2][1]=0;
    a[1][1]=0;a[1][2]=n-1;a[2][1]=1;a[2][2]=n-2;
    pow(a,K,t);mul(t,start,t);
    for (int i=1;i<=n;i++) printf("%lld
",((t[1][1]*num[i]%mod+t[2][1]*(sum-num[i])%mod)%mod+mod)%mod);
    return 0;
}