4052: [Cerc2013]Magical GCD

4052: [Cerc2013]Magical GCD

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 148  Solved: 70
[Submit][Status][Discuss]

Description

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

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

 

Input

 

Output

 

Sample Input

1 
5 
30 60 20 20 20

Sample Output

80

HINT

 

Source

题解：我会说这题暴力就水过了?= =实际上此题应该采用的方法也差不多就是暴力= =

很明显，当我们从后往前一路求最大公约数时，最多不会超过logN种（HansBug:因为要么就不变，要变就得至少少一半^_^），所以只需要记录下哪些地方可能引起公因数变化即可——

具体方法：求出两两相邻的两数的最大公约数，然后将所有的前一个公因数不是后一个的倍数的位置记录下来，后面专门留意这些位置即可（HansBug：很明显，要是前者是后者的倍数的话，那么既然都是后面的因数，则必然是前面的因数，不可能引起变化，所以记录下这些就够了，之所以不直接把原数列这么干是因为事实证明两两求下之后的数列突然变得异常优美，想想为什么^_^）

还有个萌萌哒优化：当当前的最大公因数即使到数列最前面都难以超越当前最大值的话，就可以跳掉了

（如上是对比，上面的无优化的，下面的是优化的，事实证明快了不是一点= =）

/**************************************************************
    Problem: 4052
    User: HansBug
    Language: Pascal
    Result: Accepted
    Time:1644 ms
    Memory:2180 kb
****************************************************************/
 
var
   i,j,k,l,m,n,t:longint;
   ans,v:int64;
   a,b:array[0..100005] of int64;
   c:array[0..100005] of longint;
function gcd(x,y:int64):int64;inline;
         var z:int64;
         begin
              while y<>0 do
                    begin
                         z:=x mod y;
                         x:=y;
                         y:=z;
                    end;
              exit(x);
         end;
begin
     readln(t);
     while t>0 do
           begin
                readln(n);ans:=0;m:=1;c[1]:=1;
                for i:=1 to n do
                    begin
                         read(a[i]);
                         if a[i]>ans then ans:=a[i];
                    end;
                readln;
                for i:=1 to n-1 do b[i]:=gcd(a[i],a[i+1]);
                for i:=1 to n-2 do
                    if (b[i] mod b[i+1])<>0 then
                       begin
                            inc(m);c[m]:=i+1;
                       end;
                c[m+1]:=maxlongint;c[0]:=-1;
                l:=0;
                for i:=1 to n-1 do
                    begin
                         while c[l]<=i do inc(l);
                         dec(l);v:=b[i];
                         for j:=l downto 1 do
                             begin
                                  if int64(int64(i-c[j]+2)*v)>ans then ans:=(int64(i-c[j]+2)*v);
                                  v:=gcd(v,b[c[j]-1]);
                                  if int64(int64(i+1)*v)<=ans then break;
                                  if v=1 then
                                     begin
                                          if (i+1)>ans then ans:=i+1;
                                          break;
                                     end;
                             end;
                    end;
                writeln(ans);
                dec(t);
           end;
     readln;
end.