一道非常神奇的图论题解法无比新奇清新

我们首先把图分成三种情况:

  1. 有环的,此时答案一定是环长的因数(否则不能满足题意)
  2. 存在入度大于1的DAG图的
  3. 一棵树/一条链

很容易发现,最后一种情况想怎么取就怎么取,那么我们只需要考虑1,2的影响即可

但是如果我们傻乎乎地直接跑好像也是可以的那就太烦了。

我们考虑这样的一个图:

我们把原图中的边分成两类:红色(顺时针)和绿色(逆时针)。然后我们发现这种图对应上面讲的情况2。然后手推颜色发现要两种。

而这个2是怎么来的,很简单,\(color=| total_{red}-total_{green} |\)。为什么,我们仔细观察一下对于每一个被一条顺时针边和一条逆时针边连接的点,与它相连的点颜色一定相同

然后我们就可以把问题转化成无向边并找环了。但是图中的边并没有顺时针/逆时针,于是我们分成正向边/反向边考虑即可。

可以设为相反的边权(如\(1 \&\& -1\)等)。然后取绝对值即可。

然后对于有环的情况(无论时1还是2),我们都可以得出\(ans=gcd(|len|)\),其中\(len\)表示环长,那么最小值就是\(min(ans|u)\)

值得注意的还有森林的情况,稍加推到此时的\(ans=\sum maxdis-mindis+1 \ (for\ each\ Unicom\ blocks\ of \ the\ gragh)\)

最后当得到的\(ans<3\)时无解。找环当然是BFS/DFS了(注意不要把SCC搞混了)

BFS版CODE

#include<cstdio>
#include<cctype>
#include<cstring>
using namespace std;
const int N=100005,M=1000005;
struct edge
{
int to,next,v;
}e[M<<1];
int head[N],dis[N],q[N],n,m,x,y,ans,lans,sum,cnt,MIN,MAX;
bool vis[N];
inline char tc(void)
{
static char fl[100000],*A=fl,*B=fl;
return A==B&&(B=(A=fl)+fread(fl,1,100000,stdin),A==B)?EOF:*A++;
}
inline void read(int &x)
{
x=0; char ch; while (!isdigit(ch=tc()));
while (x=(x<<3)+(x<<1)+ch-'0',isdigit(ch=tc()));
}
inline void double_add(int x,int y)
{
e[++cnt].to=y; e[cnt].v=1; e[cnt].next=head[x]; head[x]=cnt;
e[++cnt].to=x; e[cnt].v=-1; e[cnt].next=head[y]; head[y]=cnt;
}
inline int gcd(int n,int m)
{
return m?gcd(m,n%m):n;
}
inline int min(int a,int b)
{
return a<b?a:b;
}
inline int max(int a,int b)
{
return a>b?a:b;
}
inline void BFS(int s)
{
register int i,H=0,T=1; vis[s]=1; q[1]=s;
while (H<T)
{
int now=q[++H];
for (i=head[now];~i;i=e[i].next)
if (vis[e[i].to]) ans=gcd(ans,dis[now]-dis[e[i].to]+e[i].v); else
{
dis[e[i].to]=dis[now]+e[i].v; vis[e[i].to]=1; q[++T]=e[i].to;
MIN=min(MIN,dis[e[i].to]); MAX=max(MAX,dis[e[i].to]);
}
}
}
int main()
{
//freopen("CODE.in","r",stdin); freopen("CODE.out","w",stdout);
register int i; read(n); read(m);
memset(head,-1,sizeof(head)); memset(e,-1,sizeof(e));
for (i=1;i<=m;++i)
read(x),read(y),double_add(x,y);
for (i=1;i<=n;++i)
if (!vis[i])
{
MAX=MIN=0; BFS(i);
sum+=MAX-MIN+1;
}
if (ans<0) ans=-ans;
if (ans)
{
if (ans<3) return puts("-1 -1"),0;
for (lans=3;lans<ans&&ans%lans;++lans);
return printf("%d %d",ans,lans),0;
}
if (sum<3) return puts("-1 -1"),0;
return printf("%d 3",sum),0;
}

DFS版CODE

#include<cstdio>
#include<cctype>
#include<cstring>
using namespace std;
const int N=100005,M=1000005;
struct edge
{
int to,next,v;
}e[M<<1];
int head[N],dis[N],n,m,x,y,ans,lans,sum,cnt,MIN,MAX;
bool vis[N];
inline char tc(void)
{
static char fl[100000],*A=fl,*B=fl;
return A==B&&(B=(A=fl)+fread(fl,1,100000,stdin),A==B)?EOF:*A++;
}
inline void read(int &x)
{
x=0; char ch; while (!isdigit(ch=tc()));
while (x=(x<<3)+(x<<1)+ch-'0',isdigit(ch=tc()));
}
inline void double_add(int x,int y)
{
e[++cnt].to=y; e[cnt].v=1; e[cnt].next=head[x]; head[x]=cnt;
e[++cnt].to=x; e[cnt].v=-1; e[cnt].next=head[y]; head[y]=cnt;
}
inline int gcd(int n,int m)
{
return m?gcd(m,n%m):n;
}
inline int min(int a,int b)
{
return a<b?a:b;
}
inline int max(int a,int b)
{
return a>b?a:b;
}
inline void DFS(int now)
{
register int i; vis[now]=1;
for (i=head[now];~i;i=e[i].next)
if (vis[e[i].to]) ans=gcd(ans,dis[now]-dis[e[i].to]+e[i].v);
else MIN=min(MIN,dis[e[i].to]=(dis[now]+e[i].v)),MAX=max(MAX,dis[e[i].to]),DFS(e[i].to);
}
int main()
{
//freopen("CODE.in","r",stdin); freopen("CODE.out","w",stdout);
register int i; read(n); read(m);
memset(head,-1,sizeof(head)); memset(e,-1,sizeof(e));
for (i=1;i<=m;++i)
read(x),read(y),double_add(x,y);
for (i=1;i<=n;++i)
if (!vis[i])
{
MAX=MIN=0; DFS(i);
sum+=MAX-MIN+1;
}
if (ans<0) ans=-ans;
if (ans)
{
if (ans<3) return puts("-1 -1"),0;
for (lans=3;lans<ans&&ans%lans;++lans);
return printf("%d %d",ans,lans),0;
}
if (sum<3) return puts("-1 -1"),0;
return printf("%d 3",sum),0;
}

Luogu P1477 [NOI2008]假面舞会的更多相关文章

  1. 洛谷 P1477 [NOI2008]假面舞会

    题目链接 题目描述 一年一度的假面舞会又开始了,栋栋也兴致勃勃的参加了今年的舞会. 今年的面具都是主办方特别定制的.每个参加舞会的人都可以在入场时选择一 个自己喜欢的面具.每个面具都有一个编号,主办方 ...

  2. 【洛谷】1477:[NOI2008]假面舞会【图论】

    P1477 [NOI2008]假面舞会 题目描述 一年一度的假面舞会又开始了,栋栋也兴致勃勃的参加了今年的舞会. 今年的面具都是主办方特别定制的.每个参加舞会的人都可以在入场时选择一 个自己喜欢的面具 ...

  3. 【做题记录】[NOI2008] 假面舞会—有向图上的环与最长链

    luogu 1477 [NOI2008] 假面舞会 容易发现: 如果图中没有环,那么面具种数一定是所有联通块内最长链之和,最少为 \(3\) . 如果有环,则面具种数一定是所有环的大小的最大公约数. ...

  4. 图论 公约数 找环和链 BZOJ [NOI2008 假面舞会]

    BZOJ 1064: [Noi2008]假面舞会 Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 1655  Solved: 798[Submit][S ...

  5. [BZOJ1064][Noi2008]假面舞会

    [BZOJ1064][Noi2008]假面舞会 试题描述 一年一度的假面舞会又开始了,栋栋也兴致勃勃的参加了今年的舞会.今年的面具都是主办方特别定制的.每个参加舞会的人都可以在入场时选择一 个自己喜欢 ...

  6. NOI2008假面舞会

    1064: [Noi2008]假面舞会 Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 883  Solved: 462[Submit][Status] ...

  7. 【BZOJ1064】[Noi2008]假面舞会 DFS树

    [BZOJ1064][Noi2008]假面舞会 Description 一年一度的假面舞会又开始了,栋栋也兴致勃勃的参加了今年的舞会.今年的面具都是主办方特别定制的.每个参加舞会的人都可以在入场时选择 ...

  8. 1064: [Noi2008]假面舞会 - BZOJ

    Description 一年一度的假面舞会又开始了,栋栋也兴致勃勃的参加了今年的舞会.今年的面具都是主办方特别定制的.每个参加舞会的人都可以在入场时选择一 个自己喜欢的面具.每个面具都有一个编号,主办 ...

  9. 【BZOJ】1064: [Noi2008]假面舞会(判环+gcd+特殊的技巧)

    http://www.lydsy.com/JudgeOnline/problem.php?id=1064 表示想到某一种情况就不敢写下去了.... 就是找环的gcd...好可怕.. 于是膜拜了题解.. ...

随机推荐

  1. HBuilder离线打包启用Chrome Inspect调试

    解决方法: 修改这个文件 assets/data/dcloud_control.xml <msc version="1.9.9.39354" debug="true ...

  2. Loading加载页面

    一般页面有四种情况 加载中 :就是滚动页面,后台获取加载的数据,每个页面的数据不同所以就让子类来实现,直接抽象abstract了. 加载失败 :一般都需要点击后重新加载 空页面 :也需要点击后重新加载 ...

  3. maven(二):创建一个可用的maven项目,完整过程

    环境:eclipse4.5 (内置maven插件) 创建maven项目 文件菜单--新建--其他-- maven project 下一步 选择web 结构 group id:  指项目在maven本地 ...

  4. 洗礼灵魂,修炼python(30)--装饰器(2)—>装饰器总结+进阶使用

    在上一篇博文的经典案例中,我想你应该对装饰器有很好的了解了,不过光有那些还不够真的,还需要总结和进阶一下,所以本篇博文解析装饰器进阶. 装饰器 1.什么是装饰器? 个人理解:装饰器又叫语法糖,指的是对 ...

  5. python第三十六天-----类中的特殊成员方法

    __doc__ 查看尖的描述信息 __module__表示当前操作的对象所在的模块 __class__表示当前操作的对象所属的类 __init__构造方法 通过类创建对象自动执行 __del__析构方 ...

  6. css权重 vs 浏览器渲染 -- css之弊病

    昨日,突现一个bug,令人十分恼火. 基本场景 自己实现一多选日历,可多选多天(相连或不相连均可)."贵司"的需求真心有些小复杂了,"市面"上没有这样的相似的东 ...

  7. (四) 天猫精灵接入Home Assistant-ESP-WIFI模块通过mqtt协议接入HASS

    总过程 1 ESP8266上电后,初始化 连接MQTT服务器 发布自身配置信息----hass自动发现该设备 订阅hass的命令话题---接收命令 发布hass的状态话题---返回自身状态 2 ESP ...

  8. linux shell脚本调用java main方法 代码

    #!/bin/sh # #该脚本为Linux下启动java程序的通用脚本.即可以作为开机自启动service脚本被调用, #也可以作为启动java程序的独立脚本来使用. # #Author: tuda ...

  9. Redis漏洞,远程攻击

    文章转自http://blog.csdn.net/whs_321/article/details/51734602 http://blog.knownsec.com/2015/11/analysis- ...

  10. mybatis多数据源切换

    文章转自https://yq.aliyun.com/articles/188540?t=t1 https://www.liaoxuefeng.com/article/00151054582348974 ...