【图论】Codeforces 711D Directed Roads

题目链接：

　　http://codeforces.com/problemset/problem/711/D

题目大意：

　　给一张N个点N条有向边的图，边可以逆向。问任意逆向若干条边使得这张图无环的方案数(mod 1e9+7)。

题目思路：

　　【图论】

　　因为是N条边所以不会有复杂的环，最多只会有若干组一个环加一条链。

　　推算得到，一个大小为k的环对答案的贡献是*(2^k-2)，而长度为k的链对答案的贡献是2^k（链不包括环上的）

　　用dfs找出每一组环的大小和链的长度，计算答案即可。

//
//by coolxxx
//#include<bits/stdc++.h>
#include<iostream>
#include<algorithm>
#include<string>
#include<iomanip>
#include<map>
#include<stack>
#include<queue>
#include<set>
#include<bitset>
#include<memory.h>
#include<time.h>
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
//#include<stdbool.h>
#include<math.h>
#define min(a,b) ((a)<(b)?(a):(b))
#define max(a,b) ((a)>(b)?(a):(b))
#define abs(a) ((a)>0?(a):(-(a)))
#define lowbit(a) (a&(-a))
#define sqr(a) ((a)*(a))
#define swap(a,b) ((a)^=(b),(b)^=(a),(a)^=(b))
#define mem(a,b) memset(a,b,sizeof(a))
#define eps (1e-8)
#define J 10
#define mod 1000000007
#define MAX 0x7f7f7f7f
#define PI 3.14159265358979323
#define N 200014
using namespace std;
typedef long long LL;
int cas,cass;
int n,m,lll,ans;
LL aans;
LL e[N];
int to[N];
int t[N];
bool mark[N];
LL mi(int x,int y)
{
	LL sum=1;
	while(y)
	{
		if(y&1)sum=(sum*x)%mod;
		x=(x*x)%mod;y>>=1;
	}
	return sum;
}
void dfs(int u)
{
	while(!mark[u])
	{
		mark[u]=1;t[u]=++cas;
		u=to[u];
	}
	if(t[u]<=cass)
		aans=(aans*e[cas-cass])%mod;
	else
	{
		aans=(aans*((e[cas-t[u]+1]-2+mod)%mod))%mod;
		aans=(aans*e[t[u]-1-cass])%mod;
	}
	return;
}
void init()
{
	int i;
	e[0]=1;
	for(i=1;i<N;i++)e[i]=(e[i-1]*2)%mod;
}
int main()
{
	#ifndef ONLINE_JUDGE
//	freopen("1.txt","r",stdin);
//	freopen("2.txt","w",stdout);
	#endif
	int i,j,k;
	init();
//	for(scanf("%d",&cass);cass;cass--)
//	for(scanf("%d",&cas),cass=1;cass<=cas;cass++)
//	while(~scanf("%s",s+1))
	while(~scanf("%d",&n))
	{
		mem(mark,0);
		for(i=1;i<=n;i++)scanf("%d",&to[i]);
		aans=1;cas=0;
		for(i=1;i<=n;i++)
		{
			if(mark[i])continue;
			cass=cas;
			dfs(i);
		}
		printf("%I64d\n",aans);
	}
	return 0;
}
/*
//

//
*/