bzoj 5177: [Jsoi2013]贪心的导游

题意：

求区间$mod\ p_i$最大值。

题解：

搞笑的题……
因为$p$不大，$a_i$也不大，所以可以枚举$p$的倍数，二分+主席树找区间这个范围的最大值。
这样是$log^2$，好像复杂度炸的一批，然而成功在10秒内过了，很搞笑，应该是数据水。
直接在主席树上二分就少一个$log$，然而并没有快多少。
code：

#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<cstring>
using namespace std;
int n,m;
struct trnode{
    int lc,rc,c;
}tr[12000000];int tot=0,root[1000010];
int a[1000010];
int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
void update(int &x,int fr,int l,int r,int k)
{
    x=++tot;tr[x]=tr[fr];
    if(l==r) {tr[x].c++;return;}
    int mid=(l+r)/2;
    if(k<=mid) update(tr[x].lc,tr[fr].lc,l,mid,k);
    else update(tr[x].rc,tr[fr].rc,mid+1,r,k);
    tr[x].c=tr[tr[x].lc].c+tr[tr[x].rc].c;
}
bool flag;int ans=0;
int findans(int x1,int x2,int l,int r,int L)
{
    while(l<r)
    {
        int mid=(l+r)/2,rc=tr[tr[x2].rc].c-tr[tr[x1].rc].c;
        if(rc>0) l=mid+1,x1=tr[x1].rc,x2=tr[x2].rc;
        else r=mid,x1=tr[x1].lc,x2=tr[x2].lc;
    }
    ans=max(ans,l-L);
}
const int inf=1000;
void solve(int x1,int x2,int l,int r,int fl,int fr,int L)
{
    if(tr[x2].c-tr[x1].c==0) return;
    if(l==fl&&r==fr) {flag=true;findans(x1,x2,l,r,L);return;}
    int mid=(l+r)/2;
    if(fl>mid) solve(tr[x1].rc,tr[x2].rc,mid+1,r,fl,fr,L);
    else if(fr<=mid) solve(tr[x1].lc,tr[x2].lc,l,mid,fl,fr,L);
    else
    {
        solve(tr[x1].rc,tr[x2].rc,mid+1,r,mid+1,fr,L);
        if(flag) return;
        solve(tr[x1].lc,tr[x2].lc,l,mid,fl,mid,L);
    }
}
int main()
{
    n=read();m=read();
    for(int i=1;i<=n;i++)
    {
        a[i]=read();
        update(root[i],root[i-1],0,inf,a[i]);
    }
    while(m--)
    {
        int l=read()+1,r=read()+1,p=read();
        if(l>r) swap(l,r);
        int P=0;ans=0;
        while(P<=inf)
        {
            flag=false;
            solve(root[l-1],root[r],0,inf,P,min(inf,P+p-1),P),P+=p;
            if(ans==p-1) break;
        }
        printf("%d\n",ans);
    }
}