BZOJ 3037 创世纪

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

---

【题目分析】
基环树DP，对树的情况正常dp，如果遇到了环，随意拆开一个，考虑断开的点去还是不去。求最大的话比较麻烦，所以求一个最小支配集（就是min），在用总共的点数减去就好了。

---

【代码】

#include <cstdio>
#include <cstring>
#include <queue>
using namespace std;
int en=0,p,vis[1000001],to[1000001],n,h[1000001],ne[1000001],v[1000001],ans=0,b;
int f[1000001],g[1000001];
inline void add(int a,int b)
{v[en]=b;ne[en]=h[a];h[a]=en++;}
inline void dfs(int k)
{
    vis[k]=1;
    if (vis[to[k]]) p=k;
    else dfs(to[k]);
}
void solve(int now,int fa)
{
    f[now]=1; g[now]=1e9; vis[now]=2;
    if (now==b)g[now]=0;
    int u=h[now];
    while (u>=0)
    {
        if (v[u]!=p)
        {
            solve(v[u],now);
            g[now]+=min(f[v[u]],g[v[u]]);
            g[now]=min(g[now],f[now]+f[v[u]]-1);
            f[now]+=min(f[v[u]],g[v[u]]);
        }
        u=ne[u];
    }
}
int main()
{
    memset(h,-1,sizeof h);
    scanf("%d",&n);
    for (int i=1;i<=n;++i)
    {scanf("%d",&to[i]);add(to[i],i);}
    for (int i=1;i<=n;++i)
    {
        if (!vis[i])
        {
            dfs(i);
            b=to[p];
            solve(p,0);
            int temp=f[p];
            b=0;
            solve(p,0);
            ans+=min(temp,g[p]);
        }
    }
    printf("%d\n",n-ans);
}