其他分享
首页 > 其他分享> > bzoj 4869: [Shoi2017]相逢是问候

bzoj 4869: [Shoi2017]相逢是问候

作者:互联网

Description

Informatikverbindetdichundmich. 信息将你我连结。B君希望以维护一个长度为n的数组,这个数组的下标为从1到n的正整数。一共有m个操作,可以 分为两种:0 l r表示将第l个到第r个数(al,al+1,...,ar)中的每一个数ai替换为c^ai,即c的ai次方,其中c是 输入的一个常数,也就是执行赋值ai=c^ai1 l r求第l个到第r个数的和,也就是输出:sigma(ai),l<=i<=rai因为 这个结果可能会很大,所以你只需要输出结果mod p的值即可。  

Input

第一行有三个整数n,m,p,c,所有整数含义见问题描述。 接下来一行n个整数,表示a数组的初始值。 接下来m行,每行三个整数,其中第一个整数表示了操作的类型。 如果是0的话,表示这是一个修改操作,操作的参数为l,r。 如果是1的话,表示这是一个询问操作,操作的参数为l,r。 1 ≤ n ≤ 50000, 1 ≤ m ≤ 50000, 1 ≤ p ≤ 100000000, 0 < c <p, 0 ≤ ai < p  

Output

对于每个询问操作,输出一行,包括一个整数表示答案mod p的值。  

Sample Input

4 4 7 2
1 2 3 4
0 1 4
1 2 4
0 1 4
1 1 3

Sample Output

0
3

HINT

 

 鸣谢多名网友提供正确数据,已重测!

  思路: 这个题目也是扩展欧拉定理的应用,只不过和线段树结合起来了。  根据扩展欧拉定理的性质,我们知道题目的操作经过几次以后就会让每个ai变成定值(不超过log次),所以我们做区间修改的时候只要暴力修改就可以了。  也就是我们最多只要修改nlogn次,每次修改的时候都要用扩展欧拉定理重新求一下值,由于要涉及到快速幂,所以每次重新求值就需要log*log 这样总的时间复制度就是n*log^3。 比较容易被卡。  我们可以对幂次做预处理,这样时间复杂度可以优化到n*log^2,这样就很快了。 
  1 #include<bits/stdc++.h>
  2 using namespace std;   
  3 #define LL long long 
  4 #define lc (x<<1) 
  5 #define rc (x<<1|1)  
  6 #define mid (l+r)/2  
  7 int const N=50000+10;  
  8 int c[40][10001],cc[40][10001],p,n,m,C,a[N],s[N<<2],sum[N<<2],phi[40],num;   
  9 int ph(int x){
 10     int ret=x;  
 11     for(int i=2;i*i<=x;i++)  
 12         if(x%i==0){
 13             ret=ret/i*(i-1);  
 14             while (x%i==0) x/=i;  
 15         }
 16     if(x>1) ret=ret/x*(x-1);  
 17     return ret;  
 18 }
 19 void build(int x,int l,int r){
 20     if(l==r){
 21         sum[x]=a[l]%phi[0];  return;  
 22     }
 23     build(lc,l,mid);  
 24     build(rc,mid+1,r);  
 25     sum[x]=(sum[lc]+sum[rc])%phi[0];    
 26 }
 27 
 28 int calc(unsigned long long  x,int p){
 29     int check=0;   
 30     int tx=x,tp=p;  
 31 //    if(x>=phi[p+1]) x=x%phi[p+1]+phi[p+1];  
 32     while (p>-1){
 33         if(check){
 34             x=x%phi[p+1]+phi[p+1]; 
 35             int num1=x/10000,num2=x%10000;  
 36             x=1LL*c[p][num2]*cc[p][num1]%phi[p]; 
 37         }else {
 38             unsigned long long tmp=1; 
 39             if(C!=1)  
 40             for(int i=1;i<=x;i++){
 41                 tmp=tmp*C;  
 42                 if(tmp>=phi[p]) {
 43                     check=1; 
 44                     int num1=x/10000,num2=x%10000;  
 45                     tmp=1LL*c[p][num2]*cc[p][num1]%phi[p];  
 46                     break;  
 47                 }   
 48             } 
 49             x=tmp;  
 50         }  
 51         p--; 
 52     }
 53     //if(tp>=1) cout<<tx<<" **** "<<tp<<" "<<x<<endl;   
 54     return x;   
 55 }
 56 void insert(int x,int l,int r,int ll,int rr){
 57     if(s[x]>num) return;    
 58     if(l==r){
 59         s[x]++;   
 60         sum[x]=calc(a[l],s[x]-1); 
 61         return;  
 62     }
 63     if(ll<=mid) insert(lc,l,mid,ll,rr);  
 64     if(rr>mid)  insert(rc,mid+1,r,ll,rr);  
 65     s[x]=min(s[lc],s[rc]);  
 66     sum[x]=(sum[lc]+sum[rc])%phi[0];     
 67 }
 68 int query(int x,int l,int r,int ll,int rr){
 69     if(ll<=l && r<=rr) return sum[x];  
 70     int ret=0;  
 71     if(ll<=mid) ret+=query(lc,l,mid,ll,rr);  
 72     if(rr>mid)  ret+=query(rc,mid+1,r,ll,rr);  
 73     return ret%phi[0];   
 74 }
 75 int main(){
 76     //freopen("9.in","r",stdin); 
 77     scanf("%d%d%d%d",&n,&m,&p,&C);
 78     phi[0]=p;   
 79     while (1){
 80         c[num][0]=1;  
 81         unsigned long long tc=C;  
 82         for(int i=1;i<=10000;i++) c[num][i]=1LL*c[num][i-1]*C%p; 
 83         cc[num][0]=1;    
 84         for(int i=1;i<=10000;i++) cc[num][i]=1LL*cc[num][i-1]*c[num][10000]%p; 
 85         if(p==1) break;  
 86         phi[++num]=p=ph(p);    
 87     }  
 88     phi[++num]=1;   
 89     //cout<<num<<" "<<phi[0]<<" "<<phi[1]<<endl; 
 90     for(int i=1;i<=n;i++) scanf("%d",&a[i]);  
 91     build(1,1,n);  
 92     while (m--){
 93         int opt,l,r;  
 94         scanf("%d%d%d",&opt,&l,&r);  
 95         if(opt==0)insert(1,1,n,l,r);  
 96         else printf("%d\n",query(1,1,n,l,r));  
 97         //if(opt==0) break; 
 98 
 99     }
100     return 0; 
101 }
View Code

 

 

标签:phi,int,sum,long,ai,4869,rc,Shoi2017,bzoj
来源: https://www.cnblogs.com/ZJXXCN/p/11505015.html