其他分享
首页 > 其他分享> > 【考试总结】2022-05-01

【考试总结】2022-05-01

作者:互联网

a

观察能构成合法区间的条件就是奇数成对,同时每对奇数之间没有 \(0\)

使用一类线段树维护历史信息的手段来解决本题,为了规避麻烦的标记下放,可以将每个 “偶数个奇数” 的出现区间左右端点刻画出来并使用 \(r-l+1\) 来计算

同时在每个节点要支持 \(01\) 翻转的部分维护 “对未翻转时奇数/偶数” 的权值增加量方便下方

Code Display
const int N=5e5+10;
int n,Q,a[N],his[N],ans[N];
vector<pair<int,int> > qu[N];
bool ban[N],od[N];
#define ls p<<1
#define rs p<<1|1
#define lson p<<1,l,mid
#define rson p<<1|1,mid+1,r
bool rev[N<<2];
int num[N<<2][2],delt[N<<2][2],sum[N<<2];
inline void push_up(int p){
    sum[p]=sum[ls]+sum[rs];
    num[p][0]=num[ls][0]+num[rs][0];
    num[p][1]=num[ls][1]+num[rs][1];
    return ;
}
inline void push_rev(int p){rev[p]^=1; swap(num[p][0],num[p][1]);}
inline void push_add(int p,int v,int fl){
    delt[p][rev[p]^fl]+=v;
    sum[p]+=num[p][fl]*v;
    return ;
}
inline void push_down(int p){
    for(int i=0;i<2;++i) if(delt[p][i]){
        push_add(ls,delt[p][i],i);
        push_add(rs,delt[p][i],i);
        delt[p][i]=0;
    }
    if(rev[p]){
        push_rev(ls); push_rev(rs);
        rev[p]=0;
    }
    return ;
}
inline void upd(int R,int v,int fl,int p=1,int l=1,int r=n){
    if(r<=R) return push_add(p,v,fl),void();
    int mid=(l+r)>>1; push_down(p);
    if(R>mid) upd(R,v,fl,rson);
    upd(R,v,fl,lson);
    return push_up(p);
}
inline int query(int st,int ed,int p=1,int l=1,int r=n){
    if(st<=l&&r<=ed) return sum[p];
    int mid=(l+r)>>1,res=0; push_down(p);
    if(ed>mid) res+=query(st,ed,rson);
    if(st<=mid) res+=query(st,ed,lson);
    return res;
}
inline void dfs(int p,int l,int r){
    if(!num[p][1]) return ;
    if(l==r){
        delt[p][0]=delt[p][1]=0;
        num[p][0]=num[p][1]=0;
        return ;
    }
    int mid=(l+r)>>1; push_down(p);
    dfs(lson); dfs(rson);
    return push_up(p);
}
inline void Reverse(int R,int p=1,int l=1,int r=n){
    if(r<=R) return push_rev(p);
    int mid=(l+r)>>1; push_down(p);
    if(R>mid) Reverse(R,rson);
    Reverse(R,lson);
    return push_up(p);
}
inline void build(int p,int l,int r){
    if(l==r){
        num[p][0]=1;
        sum[p]=-l;
        return ;
    }
    int mid=(l+r)>>1;
    build(lson); build(rson);
    return push_up(p);
}
inline void output(int p,int l,int r){
    if(l==r) return cout<<sum[p]<<" ",void();
    int mid=(l+r)>>1; push_down(p);
    output(lson); output(rson);
    return ;
}
#undef ls
#undef rs
#undef lson
#undef rson
signed main(){
    freopen("a.in","r",stdin); freopen("a.out","w",stdout);
    n=read(); Q=read();
    for(int i=1;i<=n;++i) a[i]=read();
    for(int i=1;i<=Q;++i){
        int l=read(),r=read();
        qu[r].emplace_back(l,i);
    }
    build(1,1,n);
    for(int i=1;i<=n;++i){
        if(!a[i]&&i) dfs(1,1,n);
        if(a[i]&1){
            upd(i,i,0);
            upd(i,-i,1);
            Reverse(i);
        }
        for(auto t:qu[i]){
            upd(i,i+1,0);
            ans[t.sec]+=query(t.fir,i);
            upd(i,-1-i,0);
        }
    }
    rep(i,1,Q) print(ans[i]);
    return 0;
}

b

对于 \(Q\) 个属于 \([1,m]\) 的 随机区间,将左端点排序后,对应的右端点最长上升子序列长度是 \(\Theta(\sqrt Q)\) 级别的

那么将所有询问分成 \(\dfrac{Q}{\sqrt Q}\) 组,每组中询问形成嵌套关系,使用势能线段树暴力扩展即可

但是据说会被卡常,所以代码写了个能过的暴力

Code Display
const int N=15010;
int n,m,Q,sum[N<<2];
struct operation{
    int l,r,v,id;
}opt[N<<1];
inline void upd(int p,int l,int r,int st,int ed){
    if(sum[p]==r-l+1) return ;
    if(st<=l&&r<=ed) return sum[p]=r-l+1,void();
    int mid=(l+r)>>1;
    if(st<=mid) upd(p<<1,l,mid,st,ed);
    if(ed>mid) upd(p<<1|1,mid+1,r,st,ed);
    sum[p]=sum[p<<1]+sum[p<<1|1];
    return ;
}
inline void build(int p,int l,int r){
    sum[p]=0; if(l==r) return ;
    int mid=(l+r)>>1;
    build(p<<1,l,mid); build(p<<1|1,mid+1,r);
}
signed main(){
    freopen("b.in","r",stdin); freopen("b.out","w",stdout);
    n=read(); m=read(); Q=read();
    for(int i=1;i<=n;++i) opt[i]={i,i,read(),0};
    int cnt=n;
    for(int i=1;i<=m;++i){
        int l=read(),r=read(),v=read();
        opt[++cnt]={l,r,v,i};
    }
    sort(opt+1,opt+cnt+1,[&](const operation A,const operation B){return A.v<B.v;});
    while(Q--){
        int l1=read(),r1=read(),l2=read(),r2=read();
        int lst=0;
        int ans=0;
        int siz=r2-l2+1;
        build(1,1,siz);
        for(int i=1;i<=cnt;++i){
            if(opt[i].id&&(opt[i].id<l1||opt[i].id>r1)) continue;
            if(opt[i].r<l2||opt[i].l>r2) continue;
            upd(1,1,siz,max(opt[i].l,l2)-l2+1,min(opt[i].r,r2)-l2+1);
            ans+=(sum[1]-lst)*opt[i].v;
            if(sum[1]==siz) break;   
            lst=sum[1];
        }
        print(ans);
    }
    return 0;
}

c

发现本质上要完成的工作是对 \(n\times K\) 个点值做背包,但是每个元素的值域只有 \(\pm1\) ,也就是在 \(n\) 个数字里面找到 \(m\) 个并将它们乘起来再对所有选择方案的权值求和

设有 \(x\) 个元素贡献 \(1\) 而 \(n-x\) 个元素贡献 \(-1\),对于所有 \(x\) 求出来 \([z^m](1-z)^{n-x}(1+z)^{x}\)

而这通过分配幂次变成 \(\displaystyle \sum_{i=0}^m(-1)^{m-i}\binom{n-x}{m-i}\binom{x}i\)

再拆成阶乘逆元可以通过一次卷积进行计算

在对 \(a_i\) 形成的桶做 \(\rm FWT\) 后得到的点值解方程可以得到上面所写的 \(x\),那么直接代换成背包后的结果再做逆变换得到答案

Code Display
const int N=2e5+10;
int n,m,K;
inline void FWT(vector<int> &f,int lim){
    for(int p=2;p<=lim;p<<=1){
        int len=p>>1;
        for(int k=0;k<lim;k+=p) for(int l=k;l<k+len;++l){
            int x=f[l],y=f[l+len];
            f[l+len]=x-y;
            f[l]=x+y;
        }
    }
    return ;
}
inline void iFWT(vector<int> &f,int lim){
    for(int p=2;p<=lim;p<<=1){
        int len=p>>1;
        for(int k=0;k<lim;k+=p) for(int l=k;l<k+len;++l){
            int x=f[l],y=f[l+len];
            f[l+len]=del(x,y);
            f[l]=add(x,y);
        }
    }
    for(int tmp=ksm(lim,mod-2),i=0;i<lim;++i) ckmul(f[i],tmp);
    return ;
}
int W[1<<20],r[1<<20];
inline void NTT(vector<int> &f,int lim,int opt){
    f.resize(lim);
    for(int i=0;i<lim;++i){
        r[i]=r[i>>1]>>1|((i&1)?(lim>>1):0);
        if(i<r[i]) swap(f[i],f[r[i]]);
    }
    for(int p=2;p<=lim;p<<=1){
        int len=p>>1; W[0]=1; 
        W[1]=ksm(3,(mod-1)/p);
        if(opt==-1) W[1]=ksm(W[1],mod-2);
        for(int j=2;j<len;++j) W[j]=mul(W[j-1],W[1]);
        for(int k=0;k<lim;k+=p){
            for(int l=k;l<k+len;++l){
                int tt=mul(f[l+len],W[l-k]);
                f[l+len]=del(f[l],tt);
                ckadd(f[l],tt);
            }
        }        
    }
    if(opt==-1) for(int tmp=ksm(lim,mod-2),i=0;i<lim;++i) ckmul(f[i],tmp);
}
int ifac[N],fac[N];
signed main(){
    freopen("c.in","r",stdin); freopen("c.out","w",stdout);
    n=read(); m=read(); K=read();
    vector<int> f,g,a(K);
    for(int i=1;i<=n;++i) a[read()]++;
    if(m==1){
        rep(i,0,K-1) print(a[i]);
        putchar('\n');
        exit(0);
    }
    fac[0]=1;
    for(int i=1;i<=n;++i) fac[i]=mul(fac[i-1],i);
    ifac[n]=ksm(fac[n],mod-2);
    for(int i=n;i>=1;--i) ifac[i-1]=mul(ifac[i],i);
    
    FWT(a,K);
    int lim=1;
    while(lim<=n) lim<<=1;
    f.resize(m+1);
    g.resize(n-m+1);
    for(int i=0;i<=m;++i){
        f[i]=mul(ifac[i],ifac[m-i]);
        if((m-i)&1) f[i]=del(0,f[i]);
    }
    for(int i=0;i<=n-m;++i) g[i]=mul(ifac[i],ifac[n-m-i]);
    NTT(f,lim,1); NTT(g,lim,1);
    for(int i=0;i<lim;++i) ckmul(f[i],g[i]);
    NTT(f,lim,-1);
    for(int i=0;i<=n;++i) ckmul(f[i],mul(fac[i],fac[n-i]));
    for(int i=0;i<K;++i) a[i]=f[(n-a[i])/2];
    iFWT(a,K);
    for(int i=0;i<K;++i) print(a[i]); putchar('\n');
    return 0;
}

标签:opt,01,return,05,int,lim,rson,2022,push
来源: https://www.cnblogs.com/yspm/p/16216788.html