codeforces575A Fibonotci

题目大意：f[k]=f[k-1]*s[(n-1)%n]+f[(k-2)]*s[(k-2)%n];会修改某一位置的s值，但循环不变，求f[k];

矩阵快速幂裸题，由于有修改，所以需要线段树优化

 #include<bits/stdc++.h>
 using namespace std;
 #define ll long long
 #define maxn 50005
 #define f(c,d) ((1<<(c))*(d))
 struct node{
     ll x[][];
 }mt,mi,t,tt,x[][maxn];
 int n,m,p;
 ll k,s[maxn],ti;
 pair<ll,ll>v[maxn];
 node operator * (node a,node b){
     node c=mt;
     for(int i=;i<;++i)
         for(int j=;j<;++j)
             for(int k=;k<;++k)
                 c.x[i][j]=(c.x[i][j]+a.x[i][k]*b.x[k][j]%p)%p;
     return c;
 }
 node operator ^ (node a,ll y){
     if(y==)return mi;
     node c=a;y--;
     while(y){
         if(y&)c=c*a;
         y>>=;a=a*a;
     }return c;
 }
 node ask(int a,int b,int c,int d){
     if(f(c,d)==a&&f(c,d+)==b)return x[c][d];
     c--;d=d<<|;
     if(b<=f(c,d))return ask(a,b,c,d-);
     if(a>=f(c,d))return ask(a,b,c,d);
     return ask(a,f(c,d),c,d-)*ask(f(c,d),b,c,d);
 }
 node ask(int a,int b){
     if(a==b)return mi;
     return ask(a,b,,);
 }
 void to(ll a){
     if(a/n-ti/n){
         t=t*ask(ti%n,n);
         t=t*(x[][]^(a/n-ti/n-));
         t=t*ask(,a%n);
     }else t=t*ask(ti%n,a%n);
     ti=a;
 }
 int main(){
     scanf("%I64d%d",&k,&p);
     mi.x[][]=mi.x[][]=%p;
     scanf("%d",&n);
     for(int i=;i<n;++i)scanf("%I64d",&s[i]);
     for(int i=;i<n;++i){
         x[][i].x[][]=s[i]%p;
         x[][i].x[][]=;
         x[][i].x[][]=s[(i+)%n]%p;
     }
     for(int i=;i<;++i)
         for(int j=;f(i,j)<n;++j){
             x[i][j]=x[i-][j<<];
             if(f(i-,(j<<|))<n)x[i][j]=x[i-][j<<]*x[i-][j<<|];
         }
     scanf("%d",&m);
     for(int i=;i<m;++i)scanf("%I64d%I64d",&v[i].first,&v[i].second);
     sort(v,v+m);
     ti=;
     t=mi;
     for(int i=;i<m;++i){
         if(v[i].first>k)break;
         if(ti!=v[i].first){
             to(v[i].first-);
             tt.x[][]=s[ti%n]%p;
             tt.x[][]=;
             tt.x[][]=v[i].second%p;
             t=t*tt;
             ti++;
         }
         if(ti==k)break;
         ti++;
         tt.x[][]=v[i].second%p;
         tt.x[][]=;
         tt.x[][]=(v[i+].first==ti?v[i+].second:s[ti%n])%p;
         t=t*tt;
     }
     to(k);
     printf("%I64d\n",t.x[][]);
     return ;
 }

575A