【NOIP模拟】数字对（RMQ，二分）

题意：小H是个善于思考的学生，现在她又在思考一个有关序列的问题。 她的面前浮现出一个长度为n的序列{ai}，她想找出一段区间[L, R](1 <= L <= R <= n)。

这个特殊区间满足，存在一个k(L <= k <= R)，并且对于任意的i(L <= i <= R)，ai都能被ak整除。这样的一个特殊区间 [L, R]价值为R – L。

小H想知道序列中所有特殊区间的最大价值是多少，而有多少个这样的区间呢？这些区间又分别是哪些呢？你能帮助她吧。

100%: 1 <= n <= 500000 , 1 <= ai < 2 ^ 31.

思路：

 #include<map>
 #include<set>
 #include<cmath>
 #include<cstdio>
 #include<vector>
 #include<cstring>
 #include<cstdlib>
 #include<iostream>
 #include<algorithm>
 #define MAXN 510000
 using namespace std;

 int f[MAXN][],g[MAXN][];
 int a[MAXN],b[MAXN],c[MAXN];
 int n,num;

 int gcd(int x,int y)
 {
     int a,b,r,t;
     a=x; b=y;
     if(a<b){t=a;a=b;b=t;}
     while()
     {
         r=a%b;
         a=b;
         b=r;
         if(!r) break;

     }
     return a;
 } 

 int query1(int x,int y)
 {
   int len,l;
   len=y-x+; l=log(len)/log();
   return gcd(f[x][l],f[y-(<<l)+][l]);
 }

 int query2(int x,int y)
 {
   int len,l;
   len=y-x+; l=log(len)/log();
   return min(g[x][l],g[y-(<<l)+][l]);
 }
 int isok(int k)
 {
     num=;
     for(int i=;i+k-<=n;i++)
     {
         int s1=query1(i,i+k-);
         int s2=query2(i,i+k-);
         if(s1==s2) b[++num]=i;
     //    printf("%d %d %d %d\n",i,i+k-1,s1,s2);
     }
     if(num) return ;
      else return ;
 }
 int main()
 {
     freopen("1.in","r",stdin);
     freopen("1.out","w",stdout);
     scanf("%d",&n);
     for(int i=;i<=n;i++)
     {
         scanf("%d",&a[i]);
         f[i][]=g[i][]=a[i];
     }
     int len=log(n)/log();
     for(int i=;i<=len;i++)
      for(int j=;j+(<<i)-<=n;j++)
      {
          f[j][i]=gcd(f[j][i-],f[j+(<<(i-))][i-]);
          g[j][i]=min(g[j][i-],g[j+(<<(i-))][i-]);
      }
     int l,r,mid,last,ans;
     l=;r=n;last=ans=;c[]=;
     while(l<=r)
     {
         int mid=(l+r)/;
         if(isok(mid))
         {
             last=mid; ans=num;
             for(int i=;i<=num;i++) c[i]=b[i];
             l=mid+;
         }
          else r=mid-;
     }
     printf("%d %d\n",ans,last-);
     for(int i=;i<=ans-;i++) printf("%d ",c[i]);
     printf("%d",c[ans]);
 }