【BZOJ1064 】[NOI2008] 假面舞会(图上DFS)

点此看题面

大致题意:有k种面具(k是一个未知数且k≥3,每种面具可能有多个),已知戴第i种面具的人能看到第i+1种面具上的编号,特殊的,戴第k种面具的人能看到第1种面具上的编号,现在用x和y来表示戴着第x号的面具的人能看到第y号面具的编号,给你m组x和y(信息可能并不完整),请你求出至多至少有多少个面具。

这道题可以近似地看作一个有向图,但是有向图在这道题目中是极难操作的,因此我们可以用一个简(xuan)单(xue)的小技巧:

add(x,y,1),add(y,x,-1);//将从x到y的有向边分成从x到y的权值为1的边和从y到x的权值为-1的边,虽说依然是有向图,但操作起来与无向图差不多,可以从两个方向走

有了这个铺垫,后面的过程就会省力许多。

其实,我们可以对这道题中的图进行一个分类讨论

第一种情况是图中只存在环

如果是个有向图,找环是个很麻烦的过程,但由于我们之前已经把这张图改成了无向图,找环就非常方便啦!

在找环的过程中,我们可以轻松计算出环的长度(用当前值减去上一次访问该节点时的值,然后取绝对值即可,具体实现见代码)。

此时,我们可以得出一个十分显然的结论:面具的种类数是所有环长的gcd的一个因数(证明?我也不知道,感性理解一下即可)。

因此,最终面具的种数的最大值应为所有环长的gcd,最小值应为该gcd大于等于3的最小因数

第二种情况则是图中不存在环,而是由若干棵树(这里我们把链也当作一棵树)组成

这个时候的答案就更好推了,最大值就是最大的树的深度,最小值就是3,当然要注意判断最大值是否小于3,小于3要输出-1。

那不就好了吗?直接上代码!

等等。。。如果原图中既有环又有树(链)呢?那该怎么办?

答案是没关系!当作只有环来做就可以了,因为树在这种情况中可以忽略不计!

代码见下:

#include
#define N 100000
#define M 1000000
using namespace std;
int n,m,ans=0,ee=0,Min,Max,sum,lnk[N+5],vis[N+5],s[N+5];
struct edge
{
    int to,nxt,val;
}e[2*M+5];
inline char tc() 
{
    static char ff[100000],*A=ff,*B=ff;
    return A==B&&(B=(A=ff)+fread(ff,1,100000,stdin),A==B)?EOF:*A++;
} 
inline void read(int &x)
{
    x=0;int f=1;char ch;
    while(!isdigit(ch=tc())) if(ch=='-') f=-1;
    while(x=(x<<3)+(x<<1)+ch-'0',isdigit(ch=tc()));
    x*=f;
}
inline void write(int x)
{
    if(x<0) putchar('-'),x=-x;
    if(x>9) write(x/10);
    putchar(x%10+'0'); 
} 
inline void add(int x,int y,int z) 
{
    e[++ee].to=y,e[ee].nxt=lnk[x],e[ee].val=z,lnk[x]=ee;
}
inline int gcd(int x,int y) 
{
    return y?gcd(y,x%y):x;
}
inline void dfs(int x)//用dfs对原图进行遍历,记录是环长gcd和最大树的深度 
{
    vis[x]=1;//标记已访问 
    for(register int i=lnk[x];i;i=e[i].nxt) 
    {
    	if(vis[e[i].to]) ans=gcd(s[x]-s[e[i].to]+e[i].val,ans);//说明有环,并更新环长gcd 
		else s[e[i].to]=s[x]+e[i].val,Min=min(Min,s[e[i].to]),Max=max(Max,s[e[i].to]),dfs(e[i].to);//更新最大树的深度,并继续往下dfs 
	}
}
int main()
{
	register int i;int x,y; 
    for(read(n),read(m),i=1;i<=m;++i) 
        read(x),read(y),add(x,y,1),add(y,x,-1);
    for(i=1;i<=n;++i) 
        if(!vis[i]) Min=Max=0,dfs(i),sum+=Max-Min+1;//若当前节点没有访问过,则对其进行dfs 
    if(ans<0) ans=-ans;//ans在一波操作后可能会小于0,若其小于0则要将其改为正数 
    if(ans)//ans不为0说明有环 
    {
        if(ans<3) return puts("-1 -1"),0;//ans小于3,说明无解,直接输出-1并退出程序 
        write(ans),putchar(' ');
        for(i=3;i<=ans;i++) 
            if(!(ans%i)) return write(i),0;//寻找环长gcd大于3的最小因数 
    }
    if(sum<3) puts("-1 -1");//判断最大深度是否小于3 
    else write(sum),putchar(' '),putchar('3');//输出答案 
    return 0;
}

 

你可能感兴趣的:(BZOJ,DFS)