关闭

[基环] codeforces 711D. Directed Roads

标签: codeforces
207人阅读 评论(0) 收藏 举报
分类:

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);
}
0
0

查看评论
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场
    个人资料
    • 访问:31805次
    • 积分:1056
    • 等级:
    • 排名:千里之外
    • 原创:77篇
    • 转载:0篇
    • 译文:0篇
    • 评论:9条
    文章分类
    最新评论