关闭

[Codeforces Round #369 (Div. 2)D. Directed Roads]Tarjan强连通分量+组合计数

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

[Codeforces Round #369 (Div. 2)D. Directed Roads]Tarjan强连通分量+组合计数

题目链接[Codeforces Round #369 (Div. 2)D. Directed Roads]
题意描述:给定N个点,N条边的有向图。 (2N2105)。可以选择某些顶点,构成一个顶点集合,然后将连接顶点集合中的所有有向边进行翻转操作。在满足翻转之后的有向图不包含环的情况下,求翻转之后的图有多少种。答案 对109+7取模。
解题思路:首先对原有向图求强连通分量,求出强连通分量的个数scc,以及每个强连通分量大小(即包含顶点的个数),令numi表示第i个强连通分量的大小(1iscc)。并求出连接各个强连通分量之间的有向边的条数t
那么答案

ans=2ti=1scc(2numi2)

为什么呢?
首先我们对第i个强连通分量进行分析。对于一个大小为numi的强连通分量。改变一些边之后让他不再构成强连通分量的方案数是2numi2,减2是因为所有的总方案数 减去 构成强连通图的方案数。
将所有的强连通图变成不强连通之后,然后我们考虑原图中连接各个强连通分量的边,边数是t,显然,这些边的方向可以任意变化。那么方案数就是2t
所以,答案就是将上面的式子进行累乘。
Tarjan算法直接套kuangbin的模板了~~。

#include <bits/stdc++.h>
using namespace std;

#define FIN             freopen("input.txt","r",stdin)
#define FOUT            freopen("output.txt","w",stdout)
#define fst             first
#define snd             second

typedef __int64 LL;
//typedef long long LL;
typedef pair<int, int> PII;

const int MAXN = 2e5 + 5;
const LL MOD = 1e9 + 7;

int N, M, K;
struct Edge {
    int v, next;
    Edge() {}
    Edge(int v, int next) : v(v), next(next) {}
} edges[MAXN];
int head[MAXN], ESZ;
int Low[MAXN], DFN[MAXN], Stack[MAXN], Belong[MAXN];
int Index, top;
int scc;
bool Instack[MAXN];
int num[MAXN];
int par[MAXN], blocks;

void init() {
    ESZ = 0;
    blocks = N;
    memset(head, -1, sizeof(head));
    memset(par, -1, sizeof(par));
}
void add_edge(int u, int v)  {
    edges[ESZ] = Edge(v, head[u]);
    head[u] = ESZ ++;
}
int Find(int x) {
    return -1 == par[x] ? x : (par[x] = Find(par[x]));
}
void Union(int x, int y) {
    x = Find(x),  y = Find(y);
    if(x != y) par[y] = x, blocks --;
}
LL quick_pow(LL a, LL b, LL mod) {
    LL ret = 1;
    while(b) {
        if(b & 1) ret = ret * a % mod;
        a  = a * a % mod;
        b >>= 1;
    }
    return ret;
}
void Tarjan(int u) {
    int v;
    Low[u] = DFN[u] = ++Index;
    Stack[top++] = u;
    Instack[u] = true;
    for(int i = head[u]; ~i; i = edges[i].next) {
        v = edges[i].v;
        if( !DFN[v] ) {
            Tarjan(v);
            if( Low[u] > Low[v] )Low[u] = Low[v];
        } else if(Instack[v] && Low[u] > DFN[v])
            Low[u] = DFN[v];
    }
    if(Low[u] == DFN[u]) {
        scc++;
        do {
            v = Stack[--top];
            Instack[v] = false;
            Belong[v] = scc;
            num[scc]++;
        } while(v != u);
    }
}
void solve() {
    memset(DFN, 0, sizeof(DFN));
    memset(Instack, false, sizeof(Instack));
    memset(num, 0, sizeof(num));
    Index = scc = top = 0;
    for(int i = 1; i <= N; i++)
        if(!DFN[i]) Tarjan(i);
    LL ans = 1;
    for(int i = 1; i <= scc; i++) {
        LL x = (quick_pow(2LL, (LL)num[i],  MOD) -  2 +  MOD) % MOD;
        x = max(x, 1LL);
        ans = x * ans % MOD;
    }
    ans = ans * quick_pow(2LL, (LL)scc - blocks, MOD)  % MOD;
    printf("%I64d\n", ans);
}
int main() {
#ifndef ONLINE_JUDGE
    FIN;
#endif // ONLINE_JUDGE
    int u, v;
    while(~scanf("%d", &N)) {
        init();
        for (u = 1; u <= N; u++) {
            scanf("%d", &v);
            add_edge(u, v);
            Union(u, v);
        }
        solve();
    }
    return 0;
}
1
0
查看评论

Codeforces 842D Vitya and Strange Lesson【逆向思维+字典树查询亦或最小值】

D. Vitya and Strange Lesson time limit per test 2 seconds memory limit per test 256 megabytes input standard input output stan...
  • mengxiang000000
  • mengxiang000000
  • 2017-08-30 16:45
  • 807

Codeforces 842 D Vitya and Strange Lesson

题目地址 题意:告诉你一个序列,然后给你m个操作,每个操作就是把这个序列异或上x,然后再得出最小的没有在这个序列中出现的数。 思路:因为异或的定义就是转成2进制的形式,每一位相同的为0,不同的为1。这样的话我们就可以把数转化为二进制存入字典树(二叉树,向左为0向右为1)中,异或的话就是把异或的数...
  • xjh_shin
  • xjh_shin
  • 2017-09-01 17:47
  • 321

Codeforces 842 D. Vitya and Strange Lesson (trie)

Description Today at the lesson Vitya learned a very interesting function — mex. Mex of a sequence of numbers is the minimum non-negative number tha...
  • qq_28954601
  • qq_28954601
  • 2017-08-30 16:17
  • 542

codeforces 842D. Vitya and Strange Lesson

http://codeforces.com/contest/842/problem/D 思路: 我把数组中的数扔进字典树里,之后我记录下每个状态出现的次数。随后我们不需要去每次都把数组异或,我们把所有查询的数异或,#include <bits/stdc++.h> #defin...
  • My_stage
  • My_stage
  • 2017-08-30 20:43
  • 256

Codeforces 842D Vitya and Strange Lesson

Codeforces 842D Vitya and Strange Lesson01字典树,数据结构题意给出长度为 n 的非负整数序列,求该序列异或 x 以后的 mex 值。思路求mex可以用trie来搞。把整个序列的trie建出来后,假如要全局异或一个x,若x的从高到低第i位为1,则trie的第i...
  • xzxxzx401
  • xzxxzx401
  • 2017-08-30 23:56
  • 165

codeforces 842D Trie 树异或不存在的最小值

D. Vitya and Strange Lesson time limit per test2 seconds memory limit per test256 megabytes inputstandard input outputstandard output Today at th...
  • Littlewhite520
  • Littlewhite520
  • 2017-09-01 15:50
  • 155

Codeforces 842C - Ilya And The Tree - 树形DP或DFS

链接:  http://codeforces.com/contest/842/problem/C题目:Ilya is very fond of graphs, especially trees. During his last trip to the forest Ilya found a very...
  • xs18952904
  • xs18952904
  • 2017-08-31 00:55
  • 278

Codeforces 842 D Vitya and Strange Lesson 线段树

传送门:Codeforces 842D 题意:给出一个长度为n的序列和m次操作,每次操作将序列中所有数异或上给定的数x,然后求序列的mex。 mex:给定序列中未出现的最小的非负数。 思路:首先不要被题目中说的每次序列都会变化迷惑,因为异或满足结合律,因此我们只要把所有的x异或起来,就相当于对...
  • lxy767087094
  • lxy767087094
  • 2017-08-30 18:03
  • 124

Codeforces 842C Ilya And The Tree(树上DP+因子个数估计)

C. Ilya And The Tree time limit per test2 seconds memory limit per test256 megabytes inputstandard input outputstandard output Ilya is very fond ...
  • Yasola
  • Yasola
  • 2017-08-30 11:05
  • 711

codeforces D. Maximum Value

D. Maximum Value time limit per test 1 second memory limit per test 256 megabytes input standard input output standard output ...
  • liyunlong41
  • liyunlong41
  • 2016-05-18 23:30
  • 237
    个人资料
    • 访问:307827次
    • 积分:5872
    • 等级:
    • 排名:第5194名
    • 原创:273篇
    • 转载:8篇
    • 译文:0篇
    • 评论:45条
    友情链接
    最新评论