【NOIP2014模拟8.17】Magical GCD

Description

对于一个由正整数组成的序列， Magical GCD 是指一个区间的长度乘以该区间内所有数字的最大公约数。给你一个序列，求出这个序列最大的 Magical GCD。

Input

单个测试点包含多组数据。

输入的第一行是一个整数T表示数据组数。

每组数据的第一行是一个整数N，描述序列长度。

接下来N个数字，描述这个序列元素A[i]。

Output

对于每组测试数据输出一行，包含一个整数，表示序列最大的 Magical GCD。

题解

暴力出正解。
注：int64！！

代码

var
  t,n,m:longint;
  ans:int64;
  a,b:array [0..100001] of int64;
  l:array[0..100001] of longint;
function gcd(o,p:int64):int64;
begin
  if p=0 then exit(o);
  exit(gcd(p,o mod p));
end;

function max(o,p:int64):int64;
begin
  if o>p then exit(o);
  exit(p);
end;

procedure init;
var
  i:longint;
begin
  readln(n);
  ans:=0; m:=1; l[1]:=1;
  for i:=1 to n do
    begin
      read(a[i]);
      ans:=max(ans,a[i]);
    end;
  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); l[m]:=i+1;
      end;
  l[0]:=-1; l[m+1]:=maxlongint;
end;

procedure main;
var
  i,j,len:longint;
  gcdt:int64;
begin
  len:=0;
  for i:=1 to n-1 do
    begin
      while l[len]<=i do inc(len);
      dec(len); gcdt:=b[i];
      for j:=len downto 1 do
        begin
          ans:=max(ans,(i-l[j]+2)*gcdt);
          gcdt:=gcd(gcdt,b[l[j]-1]);
          if gcdt=1 then
            begin
              ans:=max(ans,i+1);
              break;
            end;
        end;
    end;
end;

begin
  readln(t);
  while t>0 do
    begin
      init;
      main;
      writeln(ans);
      dec(t);
    end;
end.