题目描述:

Mouse Hunt

time limit per test

2 seconds

memory limit per test

256 megabytes

input

standard input

output

standard output

Medicine faculty of Berland State University has just finished their admission campaign. As usual, about 80%of applicants are girls and majority of them are going to live in the university dormitory for the next

4(hopefully) years.

The dormitory consists of n rooms and a single mouse! Girls decided to set mouse traps in some rooms to get rid of the horrible monster. Setting a trap in room number costs \(c_i\) burles. Rooms are numbered from 1 to n

Mouse doesn't sit in place all the time, it constantly runs. If it is in room i in second then it will run to room \(a_i\) without visiting any other rooms inbetween (i=\(a_i\)means that mouse won't leave room i) It 's second 0 in the start.If the mouse is in some room with a mouse trap in it, then the mouse get caught into this trap.

That would have been so easy if the girls actually knew where the mouse at. Unfortunately, that's not the case, mouse can be in any room from 1 to at second 0

What it the minimal total amount of burles girls can spend to set the traps in order to guarantee that the mouse will eventually be caught no matter the room it started from?

Input

The first line contains as single integers n(\(1\leq n\leq 2*10^5\)) — the number of rooms in the dormitory.

The second line contains nintegers \(c_1,c_2,...,c_n\) (1≤\(c_i\)≤\(10^4\))is the cost of setting the trap in room number i

The third line contains nintegers \(a_1,a_2,...,a_n\)(1≤\(a_i\)n) is the room the mouse will run to the next second after being in room i

Output

Print a single integer — the minimal total amount of burles girls can spend to set the traps in order to guarantee that the mouse will eventually be caught no matter the room it started from.

Examples

Input

Copy

51 2 3 2 101 3 4 3 3

Output

Copy

3

Input

Copy

41 10 2 102 4 2 2

Output

Copy

10

Input

Copy

71 1 1 1 1 1 12 2 2 3 6 7 6

Output

Copy

2

Note

In the first example it is enough to set mouse trap in rooms 1 and 4. If mouse starts in room 1.then it gets caught immideately. If mouse starts in any other room then it eventually comes to room 4.

In the second example it is enough to set mouse trap in room 2. If mouse starts in room then it gets caught immideately. If mouse starts in any other room then it runs to room 2 in second 1.

Here are the paths of the mouse from different starts from the third example:

1→2→2→…;

2→2→…;

3→2→2→…;

4→3→2→2→…;

5→6→7→6→…;

6→7→6→…;

7→6→7→…;

So it's enough to set traps in rooms 2 and 6.

思路:

题目是说各一个数组a,其中i指向ai,也就是给了一个图,这个图是老鼠行程,根据这个图我们可以看到老鼠最终会在某个点停下,要么会一直的兜圈子。要求不管老鼠一开始在哪个房间,以某种代价最小的方式放置捕鼠夹,使得最后定能捕住老鼠。求这个最小代价。

我们发现只要在老鼠停下来的地方或者环上代价最小的一个地方放置捕鼠夹就能捕住老鼠。其实这道题就可以看成是图的强连通分解,通过tarjan算法求这个图有多少个强连通块,最后求出的强连通块再进行缩点,把一个强连通块看成一个点,并存储这个强连通块上的最小代价。最后要求的点是出度为0的点,为什么?因为出度为零的点就是最后老鼠停下来的点和环(已经缩点)。

需要注意的是缩点的一些细节,记录好强连通块的数目,每个强连通块建立对应的最小代价数组,遍历一遍点求出每个强连通块对应的最小代价,和每个强连通块的出度(只要有边相连的两个点不在一个强连通块上出度就减小一。最后遍历每个强连通块,求出出度为零的强连通块的最小代价和。

代码:

#include <iostream>
#include <cstdio>
#include <stack>
#define max_n 200005
#define INF 0x3f3f3f3f
using namespace std;
//题目数据
int a[max_n];
int c[max_n];
int n;
//链式前向星
int cnt = 0;
int head[max_n];
struct
{
int v;
int next;
}e[max_n<<1];
void add(int u,int v)
{
++cnt;
e[cnt].v = v;
e[cnt].next = head[u];
head[u] = cnt;
}
//读入优化
template<typename T>
inline void read(T& x)
{
x = 0;int f = 0;char ch = getchar();
while(ch<'0'||ch>'9') {f |= (ch=='-');ch = getchar();}
while(ch>='0'&&ch<='9') {x = x*10+ch-'0';ch = getchar();}
x = f?-x:x;
}
//tarjan算法
int Bcnt = 0;
int idx = 0;
int dfn[max_n];
int low[max_n];
int instack[max_n];
int Belong[max_n];
stack<int> s;
int Deg[max_n];
int cost[max_n];
void tarjan(int u)
{
dfn[u] = low[u] = ++idx;
s.push(u);
instack[u] = 1;
int v;
for(int i = head[u];i;i=e[i].next)
{
v = e[i].v;
if(!dfn[v])
{
tarjan(v);
low[u] = min(low[u],low[v]);
}
else if(instack[v])
{
low[u] = min(low[u],dfn[v]);
}
}
if(low[u]==dfn[u])
{
Bcnt++;
do
{
v = s.top();
s.pop();
instack[v] = 0;
Belong[v] = Bcnt;
}while(u!=v);
}
}
void solve()
{
for(int i = 1;i<=n;i++)
{
cout << "d " << dfn[i] << " l " << low[i] << endl;
}
cout << "共有 " << Bcnt << " 个连通分量" << endl;
for(int i = 1;i<=Bcnt;i++)
{
cout << "第 " << i << "个 " << endl;
for(int j = 1;j<=n;j++)
{
if(Belong[j]==i)
{
cout << j << " ";
}
}
cout << endl;
}
}
int main()
{
read(n);
for(int i = 1;i<=n;i++)
{
read(c[i]);
}
for(int i = 1;i<=n;i++)
{
read(a[i]);
add(i,a[i]);//建边
}
for(int i = 1;i<=n;i++)
{
if(!dfn[i])
{
tarjan(i);
}
}
//solve();
//缩点
for(int i = 1;i<=Bcnt;i++)
{
cost[i] = INF;
}
for(int i = 1;i<=n;i++)
{
cost[Belong[i]] = min(cost[Belong[i]],c[i]);//求出每个强连通分量的最小cost
}
for(int i = 1;i<=n;i++)
{
if(Belong[i]!=Belong[a[i]])
{
Deg[Belong[i]]++;//如果连个相连的点不在一个强连通分量里面,那么他所属的强连通分量的出度加一
}
}
/*for(int i =1;i<=n;i++)
{
cout << "Deg " << i << " : " << Deg[i] << endl;
}*/
int ans = 0;
for(int i = 1;i<=Bcnt;i++) //遍历所有的强连通分量
{
//cout << "deg " << i << ":" << Deg[i] << endl;
if(Deg[i]==0)
{
ans += cost[i];//如果该连通分量的出度为一,则比在这个连通分量上放一个捕鼠夹
}
}
cout << ans << endl;
return 0;
}

参考文章:

悠悠呦~,Codeforces 1027D Mouse Hunt (强连通缩点 || DFS+并查集),https://www.cnblogs.com/00isok/p/10394721.html

最新文章

  1. 关于Ruby常用语法案例累积
  2. JSON格式互转集合
  3. Startup key combinations for Intel-based Macs
  4. cocos2d-x之多个移动的小球
  5. System.OutOfMemoryException: 内存不足。(转)
  6. IOS-简单计时器的使用
  7. java13 InputStream,Reader
  8. ESMOD北京高级时装艺术学校_百度百科
  9. switch_case,&amp;&amp;,||,条件操作符和逗号操作符,循环语句
  10. ASP.NET文件上传和下载
  11. TCP/IP传输层,你懂多少?
  12. bash报错./mq.sh: line 15: warning: here-document at line 10 delimited by end-of-file (wanted `eof&#39;)
  13. STM32F10XX存储器细节
  14. 异常驱动的开发(Exception-Driven Development)
  15. CentOS/Linux开放某些端口
  16. 基于C/S 结构的IM即时通讯软件--上篇
  17. MySQL优化--INSERT ON DUPLICATE UPDATE死锁
  18. 【CH6802】车的放置
  19. k-近邻算法-优化约会网站的配对效果
  20. python---手动实现两个有序列表的合并

热门文章

  1. ubuntu 16 typora 安装 ,14系统的不管用。。
  2. webapi+swagger ui 文档描述
  3. 百度AI文本审核API使用说明
  4. [Windows] - 在 Windows Server 2019 找不到无线网卡 之解决
  5. 手撕面试官系列(三):微服务架构Dubbo+Spring Boot+Spring Cloud
  6. BJFU-216-基于链式存储结构的图书信息表的修改
  7. 使用码云或GitHub搭建简单的个人网站
  8. nmon2influxdb+grafana:服务监控可视化部署
  9. C# 简单的定时器使用
  10. React 语法