【图论】Codeforces 711D Directed Roads

题目链接：

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


题目大意：

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

题目思路：

　　【图论】

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

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

　　用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
102f2
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
;
int to
;
int t
;
bool mark
;
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;
}
/*
//

//
*/