BZOJ 3037: 创世纪 环套树dp

3037: 创世纪

Time Limit: 5 Sec  Memory Limit: 128 MB
Submit: 524  Solved: 233
[Submit][Status][Discuss]

Description

applepi手里有一本书《创世纪》，里面记录了这样一个故事……
上帝手中有着N 种被称作“世界元素”的东西，现在他要把它们中的一部分投放到一个新的空间中去以建造世界。每种世界元素都可以限制另外一种世界元素，所以说上帝希望所有被投放的世界元素都有至少一个没有被投放的世界元素能够限制它，这样上帝就可以保持对世界的控制。
由于那个著名的有关于上帝能不能制造一块连自己都不能举起的大石头的二律背反命题，我们知道上帝不是万能的，而且不但不是万能的，他甚至有事情需要找你帮忙——上帝希望知道他最多可以投放多少种世界元素，但是他只会O(2^N) 级别的算法。虽然上帝拥有无限多的时间，但是他也是个急性子。你需要帮助上帝解决这个问题。

Input

第一行是一个整数N，表示世界元素的数目。
第二行有 N 个整数A1, A2, …, AN。Ai 表示第i 个世界元素能够限制的世界元素的编号。

Output

一个整数，表示最多可以投放的世界元素的数目。

Sample Input

6
2 3 1 3 6 5

Sample Output

3

HINT

样例说明
选择2、3、5 三个世界元素即可。分别有1、4、6 来限制它们。

数据范围与约定
对于30% 的数据，N≤10。
对于60% 的数据, N≤10^5。
对于 100% 的数据，N≤10^6，1≤Ai≤N，Ai≠i。

BJ真是水水啊 环套树dp都写不对。。。

https://www.cnblogs.com/CQzhangyu/p/6534765.html

#include<cmath>
#include<ctime>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<string>
#include<bitset>
#include<queue>
#include<set>
#include<map>
using namespace std;

typedef long long ll;

inline int read()
{
	int x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return x*f;
}
void print(int x)
{if(x<0)putchar('-'),x=-x;if(x>=10)print(x/10);putchar(x%10+'0');}

const int N=1001000;

int last[N],ecnt;
struct EDGE{int to,nt;}e[N];
inline void add(int u,int v)
{e[++ecnt]=(EDGE){v,last[u]};last[u]=ecnt;}

int Fa[N];
inline int find(int x)
{return Fa[x]==x ? x : Fa[x]=find(Fa[x]);}

int f[N][2];

int fbd;

void dp(int u)
{
	f[u][0]=0;
	int t=1<<30;
	for(int i=last[u],v;i;i=e[i].nt)
	{
		v=e[i].to;
		if(v!=fbd) dp(v);
		f[u][0]+=max(f[v][0],f[v][1]);
		t=min(t,max(f[v][1],f[v][0])-f[v][0]);
	}
	f[u][1]=f[u][0]+1-t;
}

int S[N],T[N],tot;

int main()
{
	int n=read();
	register int i,u,x,y,tmp;
	for(i=1;i<=n;++i) Fa[i]=i;
	for(i=1;i<=n;++i)
	{
		u=read();
		x=find(u);y=find(i);
		x!=y ? (add(u,i),Fa[x]=y) : (S[++tot]=u,T[tot]=i);
	}
	int ans(0);
	for(i=1;i<=tot;++i)
	{
		dp(S[i]);
		fbd=S[i];
		dp(T[i]);
		tmp=f[T[i]][1];
		f[S[i]][1]=f[S[i]][0]+1;
		dp(T[i]);
		ans+=max(tmp,f[T[i]][0]);
	}
	cout<<ans<<endl;
	return 0;
}