记录编号 335518 评测结果 WWTTTTTTTTTTTTTTTTETTETTT
题目名称 [keller战记·外传][HZOI 2015]keller的土行孙 最终得分 0
用户昵称 GravatarTabing010102 是否通过 未通过
代码语言 C++ 运行时间 45.605 s
提交时间 2016-11-02 14:14:33 内存使用 2.68 MiB
显示代码纯文本
  1. //这是搞不过去的
  2. #include <cstdio>
  3. #include <cstdlib>
  4. using namespace std;
  5. typedef long long LL;
  6. const int MOD = 1000000007;
  7. const int maxn = 50000+10;
  8. FILE *fin, *fout;
  9. int n, m, A[maxn];
  10. int gcd(int a, int b) { return b==0?a:gcd(b, a%b); }
  11. LL lcm(LL a, LL b) { return (a/gcd(a, b)*b)%MOD; }
  12. struct SEG {
  13. int gcdv[maxn*4];
  14. LL lcmv[maxn*4];
  15. SEG(int n) {
  16. this->build(1, 1, n);
  17. }
  18. void maintain(int o, int L, int R) {
  19. if(R > L) {
  20. int lc=o*2, rc=o*2+1;
  21. gcdv[o] = gcd(gcdv[lc], gcdv[rc]);
  22. LL tmp = lcm(lcmv[lc], lcmv[rc]);
  23. // tmp %= MOD;
  24. lcmv[o] = tmp;
  25. }
  26. }
  27. void build(int o, int L, int R) {
  28. if(L == R) {
  29. lcmv[o] = gcdv[o] = A[L];
  30. } else {
  31. int lc=o*2, rc=o*2+1;
  32. int M = L + (R-L)/2;
  33. build(lc, L, M);
  34. build(rc, M+1, R);
  35. }
  36. maintain(o, L, R);
  37. // printf("build(o=%d, L=%d, R=%d) gcdv=%d, lcmv=%I64d\n", o, L, R, gcdv[o], lcmv[o]);
  38. }
  39. int oL, oR;
  40. int qgcd(int o, int L, int R) {
  41. if(L == R) {
  42. return gcdv[o];
  43. } else {
  44. int lc=o*2, rc=o*2+1;
  45. int M = L + (R-L)/2;
  46. int tans = -1;
  47. if(M >= oL) tans = qgcd(lc, L, M);
  48. if(M < oR) {
  49. int tmp = qgcd(rc, M+1, R);
  50. if(tans == -1) tans = tmp;
  51. else tans = gcd(tans, tmp);
  52. }
  53. return tans;
  54. }
  55. }
  56. LL qlcm(int o, int L, int R) {
  57. if(L == R) {
  58. return lcmv[o];
  59. } else {
  60. int lc=o*2, rc=o*2+1;
  61. int M = L + (R-L)/2;
  62. LL tans = -1;
  63. if(M >= oL) tans = qlcm(lc, L, M);
  64. if(M < oR) {
  65. LL tmp = qlcm(rc, M+1, R);
  66. if(tans == -1) tans = tmp;
  67. else tans = lcm(tans, tmp);
  68. }
  69. return tans;
  70. }
  71. }
  72. void execute_command() {
  73. int l, r;
  74. fscanf(fin, "%d%d", &l, &r);
  75. oL = l; oR = r;
  76. int gv = qgcd(1, 1, n);
  77. LL lv = qlcm(1, 1, n);
  78. LL tmp = gv * lv;
  79. // printf("gv=%d, lv=%I64d, tmp=%I64d, ", gv, lv, tmp);
  80. tmp %= MOD;
  81. fprintf(fout, "%d\n", (int)tmp);
  82. // printf("ans=%d \n", (int)tmp);
  83. }
  84. };
  85. int main() {
  86. fin = fopen("Keller_T233.in", "r");
  87. fout = fopen("Keller_T233.out", "w");
  88. // fout = stdout;
  89. fscanf(fin, "%d%d", &n, &m);
  90. if(n <= 0) exit(0);
  91. for(int i = 1; i <= n; i++) fscanf(fin, "%d", &A[i]);
  92. SEG *D = new SEG(n);
  93. for(int i = 1; i <= m; i++) D->execute_command();
  94. return 0;
  95. }