Poetize4 创世纪

3037: 创世纪

Time Limit: 5 Sec Memory Limit: 128 MB
Submit: 123 Solved: 66
[Submit][Status]

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

HINT

样例说明

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

数据范围与约定

对于30% 的数据，N≤10。

对于60% 的数据, N≤10^5。

对于 100% 的数据，N≤10^6，1≤Ai≤N，Ai≠i。

Source

Poetize4

题解：

搬运题解：http://blog.csdn.net/popoqqq/article/details/39965603

题目大意：给定一张有向图，每个点有且仅有一条出边，要求若一个点x扔下去，至少存在一个保留的点y，y的出边指向x，求最多扔下去多少个点

首先原题的意思就是支配关系我们反向考虑求最少保留的点要求一个点若扔出去则必须存在一个保留的点指向它

于是这就是最小支配集不过由于是有向图所以一个点要么选择要么被子节点支配所以就只剩下2个状态了

设f[x]为以x为根的子树选择x的最小支配集 g[x]为不选择x的最小支配集

然后由于是基环树林所以我们选择一个环上的点拆掉它的出边设这个点为x 出边指向的点为y 讨论

1.若x选择则y一开始就是被支配状态 g[y]初值为0 求一遍最小支配集

2.若x不选正常求最小支配集即可

两种情况取最小值计入ans 最后输出n-ans即可

然后说一下我关于这张图的形态的理解:

因为原本每个点只有一个出度，入度不定，然后n条边，就是一个内向树（直觉上。。。）（用词不专业，不知道对不对。。。）

然后我们存的时候是反过来存的存每个点能够被哪几个点支配，这样就成了入度为1，出度不定了，就成了一个外向树

所以我们随便找到环上的一个点，然后就可以遍历到所有点。。。而我们找环的时候是顺着原来的边找的，因为这样在环的分叉上回到环。

而这个过程使得一些点没有被标记，所以在DP的时候还要加标记

一些注释写在代码里：

代码：

 #include<cstdio>

 #include<cstdlib>

 #include<cmath>

 #include<cstring>

 #include<algorithm>

 #include<iostream>

 #include<vector>

 #include<map>

 #include<set>

 #include<queue>

 #include<string>

 #define inf 1000000000

 #define maxn 1000000+100

 #define maxm 500+100

 #define eps 1e-10

 #define ll long long

 #define pa pair<int,int>

 #define for0(i,n) for(int i=0;i<=(n);i++)

 #define for1(i,n) for(int i=1;i<=(n);i++)

 #define for2(i,x,y) for(int i=(x);i<=(y);i++)

 #define for3(i,x,y) for(int i=(x);i>=(y);i--)

 #define mod 1000000007

 using namespace std;

 inline int read()

 {

     int x=,f=;char ch=getchar();

     while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}

     while(ch>=''&&ch<=''){x=*x+ch-'';ch=getchar();}

     return x*f;

 }

 int n,rt,ban,ans,a[maxn],fa[maxn],head[maxn],f[maxn],g[maxn];

 bool v[maxn];

 struct edge{int go,next;}e[maxn];

 inline void dfs(int x)

 {

     v[x]=;

     if(v[a[x]])rt=x;else dfs(a[x]);

 }

 inline void dp(int x)

 {

     f[x]=;g[x]=inf;v[x]=;

     if(x==ban)g[x]=;

     for(int i=head[x];i;i=e[i].next)

      if(i!=rt&&e[i].go!=fa[x])

       {

           int y=e[i].go;

           fa[y]=x;

         dp(y);

           g[x]+=min(f[y],g[y]);

           g[x]=min(g[x],f[x]+f[y]-);//这里巧妙的避免了使用临时变量来存储f[y]与g[y]的最小差值

           f[x]+=min(f[y],g[y]);

       }

 }

 int main()

 {

     freopen("input.txt","r",stdin);

     freopen("output.txt","w",stdout);

     n=read();int x;

     for1(i,n)e[i].go=i,a[i]=x=read(),e[i].next=head[x],head[x]=i;

     for1(i,n)

     if(!v[i])

      {

          dfs(i);//找出这棵基环树

         ban=a[rt];//选取rt作为x，a[rt]作为y，g[y]=0

         dp(rt);

         int tmp=f[rt];//限定rt必须选

         ban=;

         dp(rt);

         tmp=min(tmp,g[rt]);//不用选rt

         ans+=tmp;

      }

     printf("%d\n",n-ans);

     return ;

 }