[JZOJ3485]独立集

Posted skylee的OI博客

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了[JZOJ3485]独立集相关的知识,希望对你有一定的参考价值。

题目大意:
  给你一个序列,对序列中所有逆序对之间连一条边,问图中最大独立集为多大,有哪些点一定在最大独立集中。

思路:
  在纸上画一下发现最大独立集一定是元序列的一个LIS,最大独立集必经点就是所有LIS的公共部分。
  考虑把所有的LIS记录下来,然后构建一个DAG,DAG的割点即为LIS的公共部分。
  由于我们需要先知道所有的LIS的具体方案,只能写O(n^2)的DP记录转移,能拿60分。
  不过事实上我们也不需要知道具体的状态。
  我们可以记录一下所有LIS中的每个点在LIS上的哪个位置,看一下这个位置是否只有它一个点的。
  如果不,那么肯定可以被别的点替代。
  如果是,那么它肯定是LIS的公共部分。
  LIS可以O(n log n)求,最后判断是O(n)的,总的时间复杂度是O(n log n)。

#include<cstdio>
#include<cctype>
#include<cstring>
#include<algorithm>
inline int getint() {
    register char ch;
    while(!isdigit(ch=getchar()));
    register int x=ch^0;
    while(isdigit(ch=getchar())) x=(((x<<2)+x)<<1)+(ch^0);
    return x; 
}
const int N=100001;
int w[N],id[N],f[N],g[N],h[N];
int n;
class FenwickTree {
    private:
        int val[N];
        int lowbit(const int &x) const {
            return x&-x;
        } 
    public:
        void reset() {
            memset(val,0,sizeof val);
        }
        void modify(int p,const int &x) {
            while(p<=n) {
                val[p]=std::max(val[p],x);
                p+=lowbit(p);
            }
        }
        int query(int p) const {
            int ret=0;
            while(p) {
                ret=std::max(ret,val[p]);
                p-=lowbit(p);
            }
            return ret;
        }
};
FenwickTree t;
int main() {
    n=getint();
    for(register int i=1;i<=n;i++) {
        const int v=getint();
        w[i]=v;
        id[v]=i;
    }
    for(register int i=1;i<=n;i++) {
        f[w[i]]=t.query(w[i])+1;
        t.modify(w[i],f[w[i]]);
    }
    t.reset();
    for(register int i=n;i;i--) {
        g[w[i]]=t.query(n-w[i]+1)+1;
        t.modify(n-w[i]+1,g[w[i]]);
    }
    int ans=0;
    for(register int i=1;i<=n;i++) {
        ans=std::max(ans,f[w[i]]+g[w[i]]-1);
    }
    printf("%d\n",ans);
    for(register int i=1;i<=n;i++) {
        if(f[w[i]]+g[w[i]]-1==ans) {
            if(!~h[f[w[i]]]) continue;
            if(!h[f[w[i]]]) {
                h[f[w[i]]]=1;
            } else {
                h[f[w[i]]]=-1;
            }
        }
    }
    for(register int i=1;i<=n;i++) {
        if(f[w[i]]+g[w[i]]-1==ans&&h[f[w[i]]]==1) {
            printf("%d ",i);
        }
    }
    return 0;
}

 

以上是关于[JZOJ3485]独立集的主要内容,如果未能解决你的问题,请参考以下文章

JZOJ 3515 软件公司

[jzoj]3456.NOIP2013模拟联考3恭介的法则(rule)

SP3485

解决未能加载文件或程序集“Newtonsoft.Json ...."或它的某一个依赖项。找到的程序集清单定义与程序集引用不匹配。 (异常来自 HRESULT:0x80131040)(代码片段

sp3485推荐电路(转)

手写数字识别——基于全连接层和MNIST数据集