[基环] codeforces 711D. Directed Roads

Directed Roads

题意：

给n个点，和P[i]，表示从第i个点到P[i]有一条有向边，每条边可以至多翻转方向一次，问有多少种方法翻转后没有环。

思路：

n个点n条边，无自环，每点出度为1，这个模型应该也出过很多次了。
这个图要么就是一整个环，要么就是一个环加一棵树，要么就是前面的复制几次。
考虑一个连通块，将它分为环上的边x和环外的边y，可以看出环外的边随便翻转都不形成环，所以任意翻转，答案是$2^{|y|}$，环上的任意翻转会有顺时针和逆时针两种环，所以答案是$2^{|x|}-2$，一个连通块的答案就是$2^{|y|} * (2^{|x|}-2)$，多个连通块相互独立，直接乘起来就行。

#include<bits/stdc++.h>
using namespace std;
const int N = 2e5+5;
typedef long long ll;
const ll mod = 1e9+7;
vector<int>G[N];
int in[N];
bool hs[N];
ll qpow(ll a, ll k){
    ll res = 1;
    while(k){
        if(k&1) res = (res*a)%mod;
        a = (a*a)%mod;
        k >>= 1;
    }
    return res;
}
ll solve(int rt){
    int all = 0;
    queue<int>q, qq;
    q.push(rt); hs[rt] = 1;
    while(!q.empty()){
        int u = q.front(); q.pop();
        all++;
        if(in[u] == 1) qq.push(u);
        for(int v : G[u]){
            if(!hs[v]){
                hs[v] = 1;
                q.push(v);
            }
        }
    }
    int ciy = all;
    while(!qq.empty()){
        int u = qq.front(); qq.pop();
        --ciy;
        for(int v : G[u]){
            if((--in[v]) == 1) qq.push(v);
        }
    }
    return qpow(2, all-ciy)*((qpow(2, ciy)-2+mod)%mod)%mod;
}
int main(){
    int n;
    scanf("%d", &n);
    ll ans = 1;
    for(int x, i = 1; i <= n; ++i){
        scanf("%d", &x);
        G[x].push_back(i);
        G[i].push_back(x);
        in[x]++, in[i]++;
    }
    for(int i = 1; i <= n; ++i){
        if(hs[i]) continue;
        ans = ans*solve(i)%mod;
    }
    printf("%lld\n", ans);
}