luogu P3225 [HNOI2012]矿场搭建

传送门

刷这个题之前就觉得自己点双不太会QAQ

最后果然是点双写跪了……

题意与连通性有关+无向图 考虑tarjan->割点

自动想到分情况

对于每一个点双进行考虑

只有一个点肯定要放 方案数1

同理 如果一个割点都没有 需要放两个出口防止其中一个塌了 方案数C(sze,2)

如果只有有一个割点 那么一定要在里面放一个而且除了割点随便放（防止割点塌）方案数sze-1

如果有两个以上割点 那么一个都不用放 因为一个割点塌另一个也可以用

注意：

1.tot初值设为1……

2.我万年不用vector这次都用了(存点双所有点真恶心…)

3.tarjan点双弹栈谈到儿子（Debug的时候样例直接出一个sze=50000…）

Code:

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
#include<vector>
#define ms(a,b) memset(a,b,sizeof a)
#define rep(i,a,n) for(int i = a;i <= n;i++)
#define per(i,n,a) for(int i = n;i >= a;i--)
#define inf 2147483647
using namespace std;
typedef long long ll;
ll read() {
    ll as = 0,fu = 1;
    char c = getchar();
    while(c < '0' || c > '9') {
        if(c == '-') fu = -1;
        c = getchar();
    }
    while(c >= '0' && c <= '9') {
        as = as * 10 + c - '0';
        c = getchar();
    }
    return as * fu;
}
const int N = 505;
//head
int n,m,T;
int head[N],nxt[N<<1],mo[N<<1],cnt;
void _add(int x,int y) {
    mo[++cnt] = y;
    nxt[cnt] = head[x];
    head[x] = cnt;
}
void add(int x,int y){if(x^y)_add(x,y),_add(y,x);}

int dfn[N],low[N],stk[N],top,idx,scc;
bool cut[N];

vector<int> mem[N];

void tarjan(int x,bool rt) {
    dfn[x] = low[x] = ++idx;
    stk[++top] = x;
    int cld = 0;
    for(int i = head[x];i;i = nxt[i]) {
    int sn = mo[i];
    if(!dfn[sn]) {
        tarjan(sn,0);
        low[x] = min(low[x],low[sn]);
        if(low[sn] == dfn[x]) {
        mem[++scc].push_back(x);
        int t = -1;
        while(t != sn) {
            t = stk[top--];
            mem[scc].push_back(t);
        }
        cld++;
        cut[x] = 1;
        }
    } else low[x] = min(low[x],dfn[sn]);
    }
    if(rt && cld < 2) cut[rt] = 0;
}


void clr() {
    rep(i,1,scc) mem[i].clear();
    ms(low,0),ms(dfn,0),ms(head,0),ms(cut,0);
    cnt = top = scc = idx = n = m = 0;
}

ll C2(ll n) {return n*(n-1)/2;}

void solve() {
    clr();
    m = read();
    if(!m) exit(0);
    rep(i,1,m) {
    int x = read();
    int y = read();
    n = max(n,max(x,y));
    add(x,y);
    }
    int ans = 0;
    ll tot = 1;
    rep(i,1,n) if(!dfn[i]) top = 0,tarjan(i,1);
    rep(i,1,scc) {
    int num = 0;
    ll sze = (ll)mem[i].size();
    rep(j,0,sze-1) if(cut[mem[i][j]]) num++;
    if(sze == 1ll) ans++;
    else if(!num) ans += 2,tot *= C2(sze);
    else if(num == 1) ans++,tot *= sze - 1;
    }
    printf("Case %d: %d %lld
",++T,ans,tot);
}
    

int main() {
    while(1) solve();
}