CF359D：Pair of Numbers——题解

https://vjudge.net/problem/CodeForces-359D

http://codeforces.com/problemset/problem/359/D

题目大意：

给一串数，问一个区间内所有的数是否能被其其中一个数所全部整除，求出满足条件的区间的长度最大值，并输出这样的区间的个数与它们的左端点。

换句话将，求区间GCD=区间MIN的最大长度区间。

明显st表解决。

对于最大区间长度，二分判断即可。

（因为在poj做过类似的题所以思路能很快……就是题看不懂有点难，所以特地给出翻译）

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<vector>
using namespace std;
const int N=1e5*+;
inline int read(){
    int X=,w=; char ch=;
    while(ch<'' || ch>'') {w|=ch=='-';ch=getchar();}
    while(ch>='' && ch<='') X=(X<<)+(X<<)+(ch^),ch=getchar();
    return w?-X:X;
}
inline int qpow(int a){
    if(a==)return ;
    return (<<a);
}
int gcdd[N][];
int minn[N][];
int log[N];
int n;
int gcd(int a,int b){
    if(!b)return a;
    return gcd(b,a%b);
}
void st(){
    for(int j=;j<=log[n];j++){
    for(int i=;i<=n;i++){
        int p=i+qpow(j)-;
        if(p>n)continue;
        gcdd[i][j]=gcd(gcdd[i][j-],gcdd[i+qpow(j-)][j-]);
        minn[i][j]=min(minn[i][j-],minn[i+qpow(j-)][j-]);
    }
    }
    return;
}
bool pan(int k){
    for(int i=;i<=n;i++){
    int j=i+k;
    if(j>n)break;
    int l=log[k];
    int p=qpow(l);
    int GCD=gcd(gcdd[i][l],gcdd[j-p+][l]);
    int MIN=min(minn[i][l],minn[j-p+][l]);
    if(GCD==MIN)return ;
    }
    return ;
}
int len;
void erfen(int l,int r){
    if(l>r)return;
    int mid=(l+r)>>;
    if(pan(mid)){
    len=mid;
    erfen(mid+,r);
    }else{
    erfen(l,mid-);
    }
    return;
}
vector<int>ans;
int main(){
    n=read();
    log[]=-;
    for(int i=;i<=n;i++){
    gcdd[i][]=minn[i][]=read();
    log[i]=log[i>>]+;
    }
    log[]=;
    st();
    erfen(,n-);
    for(int i=;i<=n;i++){
    int j=i+len;
    if(j>n)break;
    int l=log[len];
    int p=qpow(l);
    int GCD=gcd(gcdd[i][l],gcdd[j-p+][l]);
    int MIN=min(minn[i][l],minn[j-p+][l]);
    if(GCD==MIN)ans.push_back(i);
    }
    printf("%d %d\n",(int)ans.size(),len);
    for(int i=;i<ans.size();i++){
    printf("%d ",ans[i]);
    }
    return ;
}