题目描述
有一天,一个名叫顺旺基的程序员从石头里诞生了。又有一天,他学会了冒泡排序和独立集。在一个图里,独立集就是一个点集,满足任意两个点之间没有边。于是他就想把这两个东西结合在一起。众所周知,独立集是需要一个图的。那么顺旺基同学创造了一个算法,从冒泡排序中产生一个无向图。
这个算法不标准的伪代码如下:
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最大独立集的大小。但是事情不止于此。顺旺基同学有时候心情会不爽,这个时候他就会要求你再回答多一个问题:最大独立集可能不是唯一的,但有些点是一定要选的,问哪些点一定会在最大独立集里。今天恰好他不爽,被他问到的同学就求助于你了。
输入
两行。第一行为N,第二行为1到N的一个全排列。
输出
两行。第一行输出最大独立集的大小,第二行从小到大输出一定在最大独立集的点的编号(输入时的序号)。
样例输入
3
3 1 2
样例输出
2
2 3
数据范围
30%的数据满足 N<=16
60%的数据满足 N<=1,000
100%的数据满足 N<=100,000
这题难想是难想,但是……
标程更难理解好不?!
所以做法直接扔标程里。
所以可能和路由器的代码风格不一致。
#include <cstdio>
#include <iostream>
#include <algorithm>
#include <map>
#include <cstring>
using namespace std;
const int MAXN = 100000;
const int INF = 1e9+7;
int B[MAXN+10], LenB;
int a[MAXN+10], b[MAXN+10];
int A[MAXN+10], LenA;
int L[MAXN+10], R[MAXN+10];
map<pair<int, int>, int> check;
int main(){
int n;
scanf("%d", &n);
for(int i=1;i<=n;i++){
scanf("%d", &a[i]);
b[i] = -a[i];//变求下降为求上升
}
//下面求出最长上升子序列
for(int i=1;i<=n;i++){
int pos = lower_bound(A+1, A+1+LenA, a[i]) - A;//找到a[i]放在A的位置
L[i] = pos;
LenA += int(A[pos] == 0);//计算A的长度
A[pos] = a[i];
}
for(int i=1;i<=n;i++) B[i] = -INF;
//下面求出最长下降子序列
for(int i=n;i>=1;i--){
int pos = lower_bound(B+1, B+1+LenB, b[i]) - B;
R[i] = pos;
LenB += int(B[pos] == -INF);
B[pos] = b[i];
}
printf("%d\n", LenA);//最长上升子序列长
for(int i=1;i<=n;i++)
if(L[i] + R[i] - 1 == LenA)
check[make_pair(L[i], R[i])]++;//创建二元组
LenB=0;
for(int i=1;i<=n;i++){
//printf("%d\n",check[make_pair(L[i], R[i])]);
if(L[i] + R[i] - 1 == LenA && check[make_pair(L[i], R[i])] == 1){
//有一个奇怪的定理
//当一个数的向左最长下降长度与向右最长上升长度
//对于所有数的左最长下降长度与向右最长上升长度中是独一无二的
//那么这个点就符合条件
//简单的证明:
//如果这个点并非唯一
//那么这个点就可以被另一个和它相同的点所替代
B[++LenB] = i;
}
}
for(int i=1;i<=LenB;i++)
printf("%d ", B[i]);
return 0;
}