bzoj 1712 [Usaco2007 China]Summing Sums 加密

1712: [Usaco2007 China]Summing Sums 加密
Time Limit: 5 Sec Memory Limit: 64 MB
Submit: 480 Solved: 188
[Submit][Status][Discuss]
Description

那N只可爱的奶牛刚刚学习了有关密码的许多算法，终于，她们创造出了属于奶牛的加密方法．由于她们并不是经验十足，她们的加密方法非常简单：第i只奶牛掌握着密码的第i个数字，起始的时候是Ci(0≤Ci<90000000).加密的时候，第i只奶牛会计算其他所有奶牛的数字和，并将这个数字和除以98765431取余．在所有奶牛计算完毕之后，每一只奶牛会用自己算得的数字代替原有的数字．也就是说，

这样，她们就完成了一次加密． 在十一月，奶牛们把这个加密法则告诉了驼鹿卡门，卡门惊呆了．之后，在一个浓雾弥漫的平安夜，卡门与奶牛们：“你们的算法十分原始，很容易就被人破解．所以你们要重复这个加密过程T(1≤T≤1414213562)次，才能达到加密效果．” 这回轮到奶牛们惊呆了．很显然，奶牛们特别讨厌做同样的无聊的事情很多次．经过了漫长的争论，卡门和奶牛们终于找到的解决办法：你被刚来加密这些数字．

Input

第1行输入N和T，之后N行每行一个整数表示初始的Ci．

Output

共N行，每行一个整数，表示T次加密之后的Ci．

Sample Input

3 4

1

0

4

INPUT DETAILS:

Three cows, with starting numbers 1, 0, and 4; four repetitions of the

encryption algorithm.

Sample Output

26

25

29

OUTPUT DETAILS:

The following is a table of the cows’ numbers for each turn:

      Cows' numbers

Turn Cow1 Cow2 Cow3

0 1 0 4

1 4 5 1

2 6 5 9

3 14 15 11

4 26 25 29

HINT

N<=50000

Source

Gold

---

【分析】

听说可以用矩阵乘法做…但是其实化简一下式子可以套上等比数列求和公式求解。

首先设$sum=\sum c[i]$，可以发现
第一次操作后
$sum'=(n-1)sum，c'[i]=sum-c[i]$

第二次操作后
$sum'=(n-1)^2sum，c'[i]=(n-1)sum-sum+c[i]$

第三次操作后
$sum'=(n-1)^3sum，c'[i]=(n-1)^2sum-(n-1)sum+sum-c[i]$

我们惊讶的发现竟然还有这种操作…
于是T次操作后

$c'[i]=( (n-1)^{T-1}-(n-1)^{T-2}+(n-1)^{T-3}-...±(n-1)^0 )sum±c[i]$

分T为奇数或者偶数讨论一番，括号里的式子把正负分开处理就可以用等比数列求和了。

(模数是素数)

---

【代码】

//bzoj 
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define ll long long
#define M(a) memset(a,0,sizeof a)
#define fo(i,j,k) for(int i=j;i<=k;i++)
using namespace std;
const int mxn=50005;
const int p=98765431;
int n,T;
ll c[mxn];
inline ll ksm(ll x,int k)
{
    if(k==0) return 1;
    if(k==1) return x;
    ll tmp=ksm(x,k>>1);
    if(k&1) return tmp*tmp%p*x%p;
    else return tmp*tmp%p;
}
int main()
{
    ll x,sum=0,a1,q;
    scanf("%d%d",&n,&T);
    T--,x=n-1;
    fo(i,1,n)
    {
        scanf("%lld",&c[i]);
        c[i]%=p,sum=(sum+c[i])%p;
    }
    if(T&1)
    {
        a1=1,q=x*x%p;
        ll tmp=a1*(ksm(q,(T+1)/2)-1+p)%p;
        tmp=tmp*ksm((q-1+p)%p,p-2)%p;
        tmp=tmp*(x-1)%p,tmp=tmp*sum%p;
        fo(i,1,n) printf("%lld\n",(tmp+c[i])%p);
    }
    else
    {
        a1=x,q=x*x%p;
        ll tmp=a1*(ksm(q,(T+1)/2)-1+p)%p;
        tmp=tmp*ksm((q-1+p)%p,p-2)%p;
        tmp=tmp*(x-1)%p,tmp=(tmp+1)%p,tmp=tmp*sum%p;
        fo(i,1,n) printf("%lld\n",(tmp-c[i]+p)%p);
    }
    return 0;
}