[bzoj1124][树形DP]枪战Maf

Description

有n个人，每个人手里有一把手枪。一开始所有人都选定一个人瞄准（有可能瞄准自己）。然后他们按某个顺序开枪，且任意时刻只有一个人开枪。因此，对于不同的开枪顺序，最后死的人也不同。

Input

输入n人数<1000000 每个人的aim

Output

你要求最后死亡数目的最小和最大可能

Sample Input

8

2 3 2 2 6 7 8 5

Sample Output

3 5

题解

显然是一个基环森林
最多能存活的人每次就是一个没有上司的舞会
在基环树上dp一下就可以了
最少能存活的
如果某个树是一个简单环，显然至少都有一个人存活
否则就只有没有入度的人能活
还卡常…
有点恶心…

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
#include<queue>
#include<vector>
#include<ctime>
#include<map>
#include<bitset>
#define LL long long
#define mp(x,y) make_pair(x,y)
using namespace std;
inline int read()
{
	int f=1,x=0;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;
}
int stack[20];
inline void write(int x)
{
    if(!x){putchar('0');return;}
    int top=0;
    while(x)stack[++top]=x%10,x/=10;
    while(top)putchar(stack[top--]+'0');
}
inline void pr1(int x){write(x);printf(" ");}
inline void pr2(int x){write(x);puts("");}
inline int _max(int x,int y){return x>y?x:y;}
int fa[1100000],rx[1100000],ry[1100000];
int findfa(int x){return fa[x]==x?fa[x]:fa[x]=findfa(fa[x]);}
struct node{int x,y,next;}a[2100000];int len,last[1100000];
void ins(int x,int y){len++;a[len].x=x;a[len].y=y;a[len].next=last[x];last[x]=len;}
int n,in[1100000],ot[1100000];
int f[1100000][2][2],ed,opt;
void dp(int x,int fa)
{
	f[x][1][opt]=1;
	if(x==ed&&opt)f[x][1][opt]=-999999999;
	for(int k=last[x];k;k=a[k].next)
	{
		int y=a[k].y;
		if(y!=fa)
		{
			dp(y,x);
			if(y!=ed)
			{
				f[x][1][opt]+=f[y][0][opt];
				f[x][0][opt]+=_max(f[y][1][opt],f[y][0][opt]);
			}
			else
			{
				if(opt)f[x][1][opt]+=f[y][0][opt];
				else f[x][1][opt]=-999999999;
				if(opt)f[x][0][opt]+=f[y][0][opt];
				else f[x][0][opt]+=_max(f[y][0][opt],f[y][1][opt]);
			}
		}
	}
}
int solve(int x)
{
	ed=ry[x];opt=1;dp(rx[x],0);
	ed=ry[x];opt=0;dp(rx[x],0);
	return _max(_max(f[rx[x]][0][1],f[rx[x]][1][1]),f[rx[x]][0][0]);
}
int ok2[1100000],ok1[1100000],tot[1100000];
int main()
{
	n=read();
	for(int i=1;i<=n;i++)fa[i]=i;
	for(int i=1;i<=n;i++)
	{
		int x=read();ot[i]++;in[x]++;
		int p=findfa(i),q=findfa(x);
		if(p!=q)
		{
			fa[p]=q,ins(i,x),ins(x,i);
			rx[q]=_max(rx[q],rx[p]);ry[q]=_max(ry[q],ry[p]);
		}
		else rx[p]=i,ry[p]=x;
	}
	for(int i=1;i<=n;i++)
	{
		int p=findfa(i);tot[p]++;
		ok1[p]+=(in[i]==0);
		if(in[i]!=1||ot[i]!=1)ok2[p]=1;
	}
	int a1=0,a2=0;
	for(int i=1;i<=n;i++)if(fa[i]==i)
	{
		int ln=tot[i];
		if(!ok2[i]&&ln!=1)a2+=1;
		else if(ln!=1)a2+=ok1[i];
		if(ln!=1)a1+=solve(i);
	}
	pr1(n-a1);pr2(n-a2);
	return 0;
}