hdu7107 GCD on Sequence

solution

从大到小枚举gcd
然后每次把gcd的倍数拿出来，枚举左端点区间。
比如有[1,10,14]，左端点区间就是[1,1],[2,10],[11,14]
然后发现，定义一个r[l]表示l最远能到达的没有确定过的边界。
然后这个一定是连续的，就是比r[l]小的他一定都是没用过的，比r[l]大的一定都是用过的区间。
线段树维护这个区间,这个r[]也是整体单调递增的。

code

#include <bits/stdc++.h>
#define FOR(i,a,b) for(int i=a;i<=b;++i)
#define ll long long
#define int long long
using namespace std;
const int _=1e6+7;
// const int mod=1e9+7;
int read() {
    int x=0,f=1;char s=getchar();
    for(;s>'9'||s<'0';s=getchar()) if(s=='-') f=-1;
    for(;s>='0'&&s<='9';s=getchar()) x=x*10+s-'0';
    return x*f;
}
int n,a[_],b[_],rk[_],ans[_];
namespace seg {
    #define ls rt<<1
    #define rs rt<<1|1
    struct node {int l,r,siz,ma,tot,lazy;}e[_<<2];
    void build(int l,int r,int rt) {
        e[rt].siz=r-l+1;
        e[rt].l=l,e[rt].r=r;
        e[rt].lazy=0;
        if(l==r) {
            e[rt].tot=n;
            e[rt].ma=n;
            return;
        }
        int mid=(l+r)>>1;
        build(l,mid,ls);
        build(mid+1,r,rs);
        e[rt].tot=(e[ls].tot+e[rs].tot);
        e[rt].ma=max(e[ls].ma,e[rs].ma);
    }
    void tag(int rt,int val) {
        e[rt].tot=e[rt].siz*val;
        e[rt].ma=val;
        e[rt].lazy=val;
    }
    void pushdown(int rt) {
        if(e[rt].lazy) {
            tag(ls,e[rt].lazy);
            tag(rs,e[rt].lazy);
            e[rt].lazy=0;
        }
    }
    void modify(int L,int R,int val,int rt) {
        // cout<<L<<" "<<R<<" "<<val<<"\n";
        if(L<=e[rt].l&&e[rt].r<=R) return tag(rt,val);
        // cout<<e[1].tot<<"!\n";
        // exit(0);
        pushdown(rt);
        int mid=(e[rt].l+e[rt].r)>>1;
        if(L<=mid) modify(L,R,val,ls);
        if(R>mid) modify(L,R,val,rs);
        e[rt].tot=(e[ls].tot+e[rs].tot);
        e[rt].ma=max(e[ls].ma,e[rs].ma);
    }
    int query(int val,int rt) {
        if(e[rt].ma<val) return 0x3f3f3f3f;
        if(e[rt].l==e[rt].r) return e[rt].l;
        pushdown(rt);
        if(e[ls].ma>=val) return query(val,ls);
        else return query(val,rs);
    }
    void dfs(int rt) {
        if(e[rt].l==e[rt].r) {
            cout<<e[rt].tot<<" ";
            return;
        }
        pushdown(rt);
        dfs(ls);
        dfs(rs);
        e[rt].tot=e[ls].tot+e[rs].tot;
        e[rt].ma=max(e[ls].tot,e[rs].tot);
    }
}
void solve() {
    n=read();
    FOR(i,1,n) a[i]=read();
    FOR(i,1,n) rk[a[i]]=i;
    seg::build(1,n,1);
    for(int i=n;i>=1;--i) {
        int len=0;
        for(int j=i;j<=n;j+=i) b[++len]=rk[j];
        sort(b+1,b+1+len);
        ans[i]=seg::e[1].tot;
        // FOR(j,1,len) cout<<b[j]<<"<<< ";
        FOR(j,1,len-1) {
            int k=seg::query(b[j+1]-1,1);
            if(k>b[j]) continue;
            // cout<<"["<<max(b[j-1]+1,k)<<","<<b[j]<<"] -> "<<b[j+1]-1<<"\n";
            seg::modify(max(b[j-1]+1,k),b[j],b[j+1]-1,1);
        }
        // seg::dfs(1);
        // cout<<"\n";
        ans[i]-=seg::e[1].tot;
    }
    FOR(i,1,n) cout<<ans[i]<<"\n";
}
signed main() {
    #ifdef ONLINE_JUDGE
    #else
        freopen("a.in","r",stdin);
        // freopen("a.out","w",stdout);
    #endif
    int T=read();
    while(T--) {
        solve();
    }   
    return 0;
}