【NOIP2013模拟联考10】独立集(bubble)题解

Description

有一天,一个名叫顺旺基的程序员从石头里诞生了。又有一天,他学会了冒泡排序和独立集。在一个图里,独立集就是一个点集,满足任意两个点之间没有边。于是他就想把这两个东西结合在一起。众所周知,独立集是需要一个图的。那么顺旺基同学创造了一个算法,从冒泡排序中产生一个无向图。

这个算法不标准的伪代码如下:
procedure bubblesortgraph(n, a[]) :

                 /*输入:点数n,1到n的全排列a。

                   输出:一个点数为n的无向图G。*/

              创建一个有n个点,0条边的无向图G。

               repeat

                    swapped = false

                    for i 从 1 到 n-1 :

                          if a[i] > a[i + 1] :

                                在G中连接点a[i]和点a[i + 1]

                                交换a[i]和a[i + 1]

                                swapped = true

               until not swapped

               输出图G。

                //结束。

那么我们要算出这个无向图G最大独立集的大小。但是事情不止于此。顺旺基同学有时候心情会不爽,这个时候他就会要求你再回答多一个问题:最大独立集可能不是唯一的,但有些点是一定要选的,问哪些点一定会在最大独立集里。今天恰好他不爽,被他问到的同学就求助于你了。

Input

两行。第一行为N,第二行为1到N的一个全排列。

Output

两行。第一行输出最大独立集的大小,第二行从小到大输出一定在最大独立集的点的编号(输入时的序号)。

Sample Input

3

3 1 2

Sample Output

2

2 3

Data Constraint

30%的数据满足 N<=16

60%的数据满足 N<=1,000

100%的数据满足 N<=100,000

思路:

最长不下降子序列及输出序列(铺垫,必学)

对于两个点i<j若a[i]<=a[j]则i,j可以在同一个独立集
所以最大的独立集应该为a的最长不下降子序列长度

引用一下OJ题解,于是配上铺垫我们只有一个疑问,如何输出一定在最大独立集的点的编号 ???

如何求一定在序列中的点呢?
把a反过来做一次最长下降子序列,两个序列取交集即可

代码:

#include<bits/stdc++.h>
#define rg register
#define max(a,b) a>b?a:b
#define Fu(i,a,b) for(int i=(a);i<=(b);i++)
#define Fd(i,a,b) for(int i=(a);i>=(b);i--)
#define fre(x) freopen(#x".in","r",stdin),freopen(#x".out","w",stdout)
using namespace std;
int read(){
    
     int f=1,x=0; char s=getchar(); while(s<'0'||s>'9') {
    
     if (s=='-') f=-1; s=getchar(); } while(s>='0'&&s<='9') {
    
     x=x*10+s-'0'; s=getchar(); } return x*f; }
int n,a[100005],f[100005],f1[1000005],u[1000005],cnt,d[100005];
int main(){
    
    
	n=read();
	Fu(i,1,n) a[i]=read();
    for(int i=1;i<=n;i++){
    
    
        int t=lower_bound(d,d+cnt+1,a[i]) - d;
        if(t>cnt) cnt=t;
        f[i]=t;
        d[t]=a[i];
    }
    cnt=0,d[0]=-(0x7fffffff);
    for(int i=n;i>=1;i--){
    
    
        int t=lower_bound(d,d+cnt+1,-a[i]) - d;
        if(t>cnt) cnt=t;
        f1[i]=t;
        d[t]=-a[i];
    }
    printf("%d\n",cnt);
    for(int i=1;i<=n;i++)
        if(f[i]+f1[i] == cnt+1)
            ++u[f[i]];
    for(int i=1;i<=n;i++)
        if(f[i]+f1[i] == cnt+1 && u[f[i]]==1)
            printf("%d ",i);
	return 0;
}

谢谢观看!!!

猜你喜欢

转载自blog.csdn.net/zhy_Learn/article/details/110245436