HDU3844Mining Your Own Business

目测某年HNOI，（其实这个题是2011年的WF，hdu上找到的，HNOI2012那个中文题在bzoj和loj上都有，叫矿场搭建，题意几乎一样，数据比较弱，交这份代码也能A）。

先讲题解，然后说一些有（e）趣（xin）的事情。

首先肯定是个点双，先求一下点双。

然后我们分析一下。

当这个点双里没割点时（一个大号的孤立点双），那么我们要建两个特殊节点，因为一个塌了可以去另一个，对第1问贡献为2，第2问贡献贡献为C(n，2)=（n-1）*n/2。（n是点双大小哦）。

当这个点双有一个割点（一个叶子），那么我们在除割点以外的地方建特殊节点，因为割点塌了去自家的，自家的塌了走割点去隔壁的。对第1问贡献为1，第2问贡献为n-1。（n是点双大小哦）。

当这个点双有两个及以上割点时，别建，天塌下来有两个及以上隔壁顶着。无贡献。

多测，注意清空。

下面说一些题外话。

这个题码完以后去bzoj交那个中文题，WA了，交loj，70，一看，错的是那三个小点，大点都过了……后来发现T清零了……case一个1，若干0，赶紧改，一交A了，回bzoj，A了。

然后交自家oj，RE，一看，数组可能要2倍，开完一交，又RE，一直交，一直RE，我去，这数据也TQL，一想，可能挂在外网上，min和max都改手打，还是RE，到处窜座去问，没人能给出解答（因为我是倒着刷的），然后去找啾啾，他上hdu给我找了一手，找到了，交，TLE，3000，一个测试点？WF?这么那啥。静态调错，init里边数没清零……明白了，那两个A的是因为数据太小，我的数组开的太大，不清空边数都没炸……那WF的题5e4的大范围，不出5个就卡死了，赶紧改，交，WA了。

静态调错，ans2*=写的赋值，点双大小写的n（知道我上文为什么要强调了吗？我推的时候拿n推的，打的时候忘了），改，交，A了。

现在又明白一个问题，bz和loj上的那道题，只有一个点双就是一个大的点双，没有孤立的，所以上面的错误就是正确的。

所以，做完这个题就不要去做那个题了，没意思了。（试了试，ans2*=C（n,2）也能A，有点迷，难不成WF也默认了？）

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<vector>
#include<queue>
#include<stack>
#include<set>
#include<map>
using namespace std;
int read(){
    int sum=,f=;char x=getchar();
    while(x<''||x>''){
        if(x=='-') f=-;
        x=getchar();
    }while(x>=''&&x<=''){
        sum=sum*+x-'';
        x=getchar();
    }return sum*f;
}
struct EDGE{
    int ed,nex;
}edge[];
int first[],num;
int n,m,root;
int dfn[],low[],sta[];
int ord,top,vccnum,T;
vector<int>vcc[];
long long ans1=,ans2=;
bool cut[];
inline int max(int a,int b){
    return a>b?a:b;
}
inline int min(int a,int b){
    return a<b?a:b;
}
void init(){
    memset(edge,,sizeof(edge));
    memset(first,,sizeof(first));
    memset(dfn,,sizeof(dfn));
    memset(low,,sizeof(low));
    memset(cut,,sizeof(cut));
    for(int i=;i<=vccnum;i++) vcc[i].clear();
    ord=top=num=vccnum=ans1=root=n=;
    ans2=;
}

void tarjan(int x){
    dfn[x]=low[x]=++ord;
    sta[++top]=x;
    if(x==root&&first[x]==){
        vcc[++vccnum].push_back(x);
        return ;
    }
    int child=;
    for(int i=first[x];i;i=edge[i].nex){
        int y=edge[i].ed;
        if(!dfn[y]){
            tarjan(y);
            low[x]=min(low[x],low[y]);
            if(low[y]>=dfn[x]){
                child++;
                if((x==root&&child>=)||(x!=root&&child>)) cut[x]=;
                vcc[++vccnum].push_back(x);
                int p;
                do{
                    p=sta[top--];
                    vcc[vccnum].push_back(p);
                }while(p!=y);
            }
        }else low[x]=min(low[x],dfn[y]);
    }
}
void add(int st,int ed){
    edge[++num].ed=ed;
    edge[num].nex=first[st];
    first[st]=num;
}
int main(){
//    freopen("c.in","r",stdin);
    m=read();
    while(m!=){
        T++;
        for(int i=,x,y;i<=m;i++){
            x=read();y=read();
            add(x,y);add(y,x);
            n=max(n,x);n=max(n,y);
        }
        for(int i=;i<=n;i++)
            if(!dfn[i]) root=i,tarjan(i);
        for(int i=;i<=vccnum;i++){
            int sum=;
            for(int j=;j<vcc[i].size();j++)
                if(cut[vcc[i][j]])
                    sum++;
            if(!sum){
                ans1+=;
                ans2*=vcc[i].size()*(vcc[i].size()-)/;
            }
            else if(sum==){
                ans1++;
                ans2*=1ll*(vcc[i].size()-);
            }
        }
        printf("Case %d: %lld %lld\n",T,ans1,ans2);
        m=read();init();
    }
}

带着浓重的戾气。