[BZOJ4052][Cerc2013]Magical GCD

[BZOJ4052][Cerc2013]Magical GCD

试题描述

给出一个长度在 100 000 以内的正整数序列，大小不超过 10^12。 

求一个连续子序列，使得在所有的连续子序列中，它们的GCD值乘以它们的长度最大。

输入

本题为多组数据。

第一行一个整数 T，表示数据组数。

每组数据第一行为一个整数 n，表示序列长度；第二行为 n 个整数表示序列。

输出

对于每组数据，输出 max{ gcd * length }

输入示例

1
5
30 60 20 20 20

输出示例

80

数据规模及约定

见“试题描述”

题解

可以发现，随着子序列的增长，gcd 要么不变，若变化则变化后的 gcd 不会超过原来的 1/2，所以 gcd 的变化次数不会超过 log₂10¹² = 12 · log₂10 次。

有了上面的结论，我们就可以枚举左端点，然后向右二分 gcd 即将变化的位置（维护区间 gcd 可以用 ST 表），在每个这样的位置更新答案，然后移动子序列右端点，重复这个过程。设数字大小为 A，那么更新答案次数为 log₂A；每次二分需要 log₂n；哦对了虽然 RMQ 询问复杂度 O(1)，但是每次需要取 gcd，最坏情况 log₂A；别忘了我们还枚举了左端点有一个 n，所以最终复杂度不会超过 O(n · log₂²A · log₂n)。

+ View Code?

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81

#include <iostream> #include <cstdio> #include <algorithm> #include <cmath> #include <stack> #include <vector> #include <queue> #include <cstring> #include <string> #include <map> #include <set> using namespace std; #define LL long long   const int BufferSize = 1 << 16; char buffer[BufferSize], *Head, *Tail; inline char Getchar() {     if(Head == Tail) {         int l = fread(buffer, 1, BufferSize, stdin);         Tail = (Head = buffer) + l;     }     return *Head++; } LL read() {     LL x = 0, f = 1; char c = Getchar();     while(!isdigit(c)){ if(c == '-') f = -1; c = Getchar(); }     while(isdigit(c)){ x = x * 10 + c - '0'; c = Getchar(); }     return x * f; }   #define maxn 100010 #define maxlog 17 int n; LL A[maxn];   LL gcd(LL a, LL b) { return b ? gcd(b, a % b) : a; }   LL Gcd[maxlog][maxn]; int Log[maxn]; void rmq_init() {     Log[1] = 0;     for(int i = 2; i <= n; i++) Log[i] = Log[i>>1] + 1;     for(int i = 1; i <= n; i++) Gcd[0][i] = A[i];     for(int j = 1; (1 << j) <= n; j++)         for(int i = 1; i + (1 << j) - 1 <= n; i++)             Gcd[j][i] = gcd(Gcd[j-1][i], Gcd[j-1][i+(1<<j-1)]);     return ; } LL query(int l, int r) {     int len = r - l + 1, t = Log[len];     return gcd(Gcd[t][l], Gcd[t][r-(1<<t)+1]); }   int main() {     int T = read();     while(T--) {         n = read();         for(int i = 1; i <= n; i++) A[i] = read();                   rmq_init();         LL ans = 0;         for(int L = 1; L <= n; L++) {             int R = L; LL tmp = query(L, R);             while(R < n) {                 int l = R, r = n + 1;                 while(r - l > 1) {                     int mid = l + r >> 1;                     if(query(L, mid) == tmp) l = mid; else r = mid;                 }                 R = l;                 ans = max(ans, tmp * (R - L + 1));                 R++; tmp = query(L, R);             }             if(R <= n) ans = max(ans, tmp * (R - L + 1)); // mark!!         }                   printf("%lld\n", ans);     }           return 0; }

大视野数据好强！我就因为少加 "mark!!" 那一行导致 WA 了一上午。。。