ACT初代奥特曼_ac自动机上dp

大家好，又见面了，我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全系列IDE使用 1年只要46元售后保障童叟无欺

上帝手中有 N 种世界元素，每种元素可以限制另外 1 种元素，把第 i 种世界元素能够限制的那种世界元素记为 A[i]。

现在，上帝要把它们中的一部分投放到一个新的空间中去建造世界。

为了世界的和平与安宁，上帝希望所有被投放的世界元素都有至少一个没有被投放的世界元素限制它。

上帝希望知道，在此前提下，他最多可以投放多少种世界元素？

输入格式
第一行是一个整数 N，表示世界元素的数目。

第二行有 N 个整数 A[1],A[2],…,A[N]。A[i] 表示第 i 个世界元素能够限制的世界元素的编号。

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

数据范围
N≤106,1≤A[i]≤N
输入样例：
6
2 3 1 3 6 5
输出样例：
3

#include<bits/stdc++.h>
using namespace std;
const int N = 2e6 + 10;
const int M = 2 * N;
const int INF = 0x0f3f3f3f;
struct Edge{ 

int v,next;
}edge[M];
int head[N],cnt;
void add(int u,int v){ 

edge[cnt].v = v;
edge[cnt].next = head[u];
head[u] = cnt ++;
}
int ins[N],st[N],cir[N],pox,pre[N],pree[N];
int f[N][2];
bool rm[M],mark[N];
int ans;
void dfs_c(int u){ 

st[u] = ins[u] = true;
for(int i = head[u];~i;i = edge[i].next){ 

if((i & 1) == 1)continue;
int v = edge[i].v;
pre[v] = u;
pree[v] = i;
if(!st[v])dfs_c(v);
else if(ins[v]){ 

for(int k = u;k != v;k = pre[k]){ 

cir[pox ++] = k;
}
cir[pox ++] = v;
}
}
ins[u] = false;
}
void dfs(int u){ 

f[u][0] = 0,f[u][1] = -INF;
int m = -INF;
for(int i = head[u];~i;i = edge[i].next){ 

int v = edge[i].v;
if((i & 1) == 0 || rm[i])continue;
dfs(v);
f[u][0] += max(f[v][1],f[v][0]);
m = max(m,-max(f[v][1],f[v][0]) + f[v][0]);
}
if(mark[u])f[u][1] = max(f[u][0] + m,max(f[u][0] + 1,1));
else f[u][1] = f[u][0] + m + 1;
}
int main(){ 

int n;
memset(head,-1,sizeof head);
scanf("%d",&n);
int x;
for(int i = 1;i <= n;i ++){ 

scanf("%d",&x);
add(i,x),add(x,i);
}
int ans = 0;
for(int i = 1;i <= n;i ++){ 

int res = 0;
if(!st[i]){ 

pox = 0;
dfs_c(i);
if(pox >= 1){ 

int u = cir[0];
int p = pre[u];
rm[pree[u] ^ 1] = true;
mark[u] = true;
memset(f,0,sizeof f);
dfs(p);
res = max(res,f[p][0]);
mark[u] = false;
memset(f,0,sizeof f);
dfs(p);
res = max(res,f[p][1]);
rm[pree[u] ^ 1] = false;
}
}
ans += res;
}
cout<<ans<<endl;
return 0;
}