记录编号 326358 评测结果 AAAAAATAAT
题目名称 学姐的巧克力盒 最终得分 80
用户昵称 GravatarAntiLeaf 是否通过 未通过
代码语言 C++ 运行时间 8.742 s
提交时间 2016-10-21 07:31:59 内存使用 30.47 MiB
显示代码纯文本
  1. #include<cstdio>
  2. #include<cstring>
  3. #include<cmath>
  4. #include<algorithm>
  5. #define LL long long
  6. using namespace std;
  7. const int maxn=1100010;//1048576
  8. LL euler_phi(LL);
  9. LL qmul(int,int,const LL*,LL);
  10. LL qpow(LL,LL,LL);
  11. int n,m,M=1,l,r,d;
  12. LL k,a[maxn<<1],b[maxn<<1],p1,p2,phi;
  13. int main(){
  14. #define MINE
  15. #ifdef MINE
  16. freopen("chocolatebox.in","r",stdin);
  17. freopen("chocolatebox.out","w",stdout);
  18. #endif
  19. scanf("%d%d%lld%lld%lld",&n,&m,&k,&p1,&p2);
  20. if(!p1)p1=1ll;
  21. if(!p2)p2=1ll;
  22. while(M<n+2)M<<=1;
  23. phi=euler_phi(p2);
  24. for(int i=1;i<=n;i++){
  25. scanf("%lld",&a[i+M]);
  26. b[i+M]=a[i+M]%phi;
  27. a[i+M]%=p1;
  28. }
  29. for(int i=n+M+1;i<(M<<1);i++)a[i]=b[i]=1ll;
  30. a[M]=b[M]=1ll;
  31. for(int i=M-1;i;i--){
  32. a[i]=a[i<<1]*a[i<<1|1]%p1;
  33. b[i]=b[i<<1]*b[i<<1|1]%phi;
  34. }
  35. while(m--){
  36. scanf("%d%d%d",&d,&l,&r);
  37. if(d==1)printf("%lld\n",qmul(l,r,a,p1));
  38. else printf("%lld\n",qpow(k,qmul(l,r,b,phi),p2));
  39. }
  40. //printf("\n");for(int i=1;i<(M<<1);i++)printf("%lld ",a[i]);
  41. #ifndef MINE
  42. printf("\n-------------------------DONE-------------------------\n");
  43. for(;;);
  44. #endif
  45. return 0;
  46. }
  47. LL euler_phi(LL n){
  48. LL ans=n;
  49. int m=(int)(sqrt(n)+0.5);
  50. for(int i=2;i<=m;i++)if(n%i==0){
  51. ans=ans/i*(i-1);
  52. do n/=i;while(n%i==0);
  53. }
  54. if(n>1)ans=ans/n*(n-1);
  55. return ans;
  56. }
  57. LL qmul(int l,int r,const LL *a,LL p){
  58. LL ans=1ll;
  59. for(l+=M-1,r+=M+1;l^r^1;l>>=1,r>>=1){
  60. if(~l&1)ans=ans*a[l^1]%p;
  61. if(r&1)ans=ans*a[r^1]%p;
  62. }
  63. return ans;
  64. }
  65. LL qpow(LL a,LL b,LL p){
  66. LL ans=1ll;
  67. for(;b;b>>=1,a=a*a%p)if(b&1ll)ans=ans*a%p;
  68. return ans%p;
  69. }