http://www.lydsy.com/JudgeOnline/problem.php?id=4424

图是二分图的条件:没有奇环

所以,如果图不存在奇环,删除任意一条边都可以

如果存在奇环,

对于树边来说:

那么可能可以删除的边一定在所有奇环的交集内

而且这条边不能在偶环内

因为如果一条边既是奇环上的一条边,又是偶环上的一条边

删除这条边后,这个奇环和偶环会合并成一个新的奇环

所以最终的答案= 奇环的交集-偶环的并集

对于非树边来说:

如果只有一个奇环,那么可以删除构成环的这条非树边

树边和非树边的判断:并查集

奇环交与偶环并:差分,统计树上后缀和

#include<cmath>
#include<cstdio>
#include<iostream>
#include<algorithm> using namespace std; #define N 1000001 int e[N][],ty[N]; int tot=;
int front[N],nxt[N<<],to[N<<],id[N<<]; int F[N]; int maxd;
int fa[N][]; int dep[N]; int odd_cir,sum[N]; int ans_tot,ans[N]; void read(int &x)
{
x=; char c=getchar();
while(!isdigit(c)) c=getchar();
while(isdigit(c)) { x=x*+c-''; c=getchar(); }
} int find(int i) { return F[i]==i ? i : F[i]=find(F[i]); } void add(int u,int v,int num)
{
to[++tot]=v; nxt[tot]=front[u]; front[u]=tot; id[tot]=num;
to[++tot]=u; nxt[tot]=front[v]; front[v]=tot; id[tot]=num;
} void dfs(int u)
{
for(int i=;i<=maxd;++i) fa[u][i]=fa[fa[u][i-]][i-];
int t;
for(int i=front[u];i;i=nxt[i])
{
t=to[i];
if(t==fa[u][]) continue;
dep[t]=dep[u]+;
fa[t][]=u;
dfs(t);
}
} int cal(int x)
{
for(int i=front[x];i;i=nxt[i])
if(to[i]!=fa[x][])
{
cal(to[i]);
if(sum[to[i]]==odd_cir) ans[++ans_tot]=id[i];
sum[x]+=sum[to[i]];
}
} int get_lca(int u,int v)
{
if(dep[u]<dep[v]) swap(u,v);
int cnt=dep[u]-dep[v];
for(int i=maxd;i>=;--i)
if(cnt&(<<i)) u=fa[u][i];
if(u==v) return u;
for(int i=maxd;i>=;--i)
if(fa[u][i]!=fa[v][i]) u=fa[u][i],v=fa[v][i];
return fa[u][];
} int main()
{
int n,m;
int u,v,fu,fv;
read(n); read(m);
for(int i=;i<=n;++i) F[i]=i;
for(int i=;i<=m;++i)
{
read(u); read(v);
fu=find(u); fv=find(v);
if(fu!=fv)
{
ty[i]=;
F[fu]=fv;
add(u,v,i);
}
e[i][]=u; e[i][]=v;
}
maxd=log(n)/log();
for(int i=;i<=n;++i)
if(find(i)==i) dep[i]=,dfs(i);
int self=;
int lca;
for(int i=;i<=m;++i)
{
u=e[i][]; v=e[i][];
if(u==v && self) { printf(""); return ; }
if(u==v) { self=i; continue; }
if(!ty[i])
{
lca=get_lca(u,v);
if((dep[u]-dep[lca]+dep[v]-dep[lca])&)
{
sum[v]--;
sum[u]--;
sum[lca]+=;
}
else
{
ty[i]=;
odd_cir++;
sum[v]++;
sum[u]++;
sum[lca]-=;
}
}
}
if(self)
{
if(!odd_cir) printf("1\n%d",self);
else printf("");
return ;
}
if(!odd_cir)
{
printf("%d\n",m);
for(int i=;i<m;++i) printf("%d ",i);
if(m) printf("%d",m);
return ;
}
for(int i=;i<=n;++i)
if(find(i)==i) cal(i);
if(odd_cir==)
for(int i=;i<=m;++i)
if(ty[i]==)
{
ans[++ans_tot]=i;
break;
}
sort(ans+,ans+ans_tot+);
printf("%d\n",ans_tot);
for(int i=;i<ans_tot;++i) printf("%d ",ans[i]);
if(ans_tot) printf("%d",ans[ans_tot]);
}

4424: Cf19E Fairy

Time Limit: 10 Sec  Memory Limit: 256 MB
Submit: 740  Solved: 187
[Submit][Status][Discuss]

Description

给定 n 个点,m 条边的无向图,可以从图中删除一条边,问删除哪些边可以使图变成
一个二分图。

Input

第 1 行包含两个整数 n,m。分别表示点数和边数。
第 2 到 m+1 行每行两个数 x,y 表示有一条(x,y)的边。

Output

输出第一行一个整数,表示能删除的边的个数。
接下来一行按照从小到大的顺序输出边的序号。

Sample Input

4 4
1 2
1 3
2 4
3 4

Sample Output

4
1 2 3 4

HINT

100%的数据,n,m<=1000000

bzoj千题计划229:bzoj4424: Cf19E Fairy的更多相关文章

  1. bzoj千题计划300:bzoj4823: [Cqoi2017]老C的方块

    http://www.lydsy.com/JudgeOnline/problem.php?id=4823 讨厌的形状就是四联通图 且左右各连一个方块 那么破坏所有满足条件的四联通就好了 按上图方式染色 ...

  2. bzoj千题计划196:bzoj4826: [Hnoi2017]影魔

    http://www.lydsy.com/JudgeOnline/problem.php?id=4826 吐槽一下bzoj这道题的排版是真丑... 我还是粘洛谷的题面吧... 提供p1的攻击力:i,j ...

  3. bzoj千题计划280:bzoj4592: [Shoi2015]脑洞治疗仪

    http://www.lydsy.com/JudgeOnline/problem.php?id=4592 注意操作1 先挖再补,就是补的范围可以包含挖的范围 SHOI2015 的题 略水啊(逃) #i ...

  4. bzoj千题计划177:bzoj1858: [Scoi2010]序列操作

    http://www.lydsy.com/JudgeOnline/problem.php?id=1858 2018 自己写的第1题,一遍过 ^_^ 元旦快乐 #include<cstdio> ...

  5. bzoj千题计划317:bzoj4650: [Noi2016]优秀的拆分(后缀数组+差分)

    https://www.lydsy.com/JudgeOnline/problem.php?id=4650 如果能够预处理出 suf[i] 以i结尾的形式为AA的子串个数 pre[i] 以i开头的形式 ...

  6. bzoj千题计划304:bzoj3676: [Apio2014]回文串(回文自动机)

    https://www.lydsy.com/JudgeOnline/problem.php?id=3676 回文自动机模板题 4年前的APIO如今竟沦为模板,,,╮(╯▽╰)╭,唉 #include& ...

  7. bzoj千题计划292:bzoj2244: [SDOI2011]拦截导弹

    http://www.lydsy.com/JudgeOnline/problem.php?id=2244 每枚导弹成功拦截的概率 = 包含它的最长上升子序列个数/最长上升子序列总个数 pre_len ...

  8. bzoj千题计划278:bzoj4590: [Shoi2015]自动刷题机

    http://www.lydsy.com/JudgeOnline/problem.php?id=4590 二分 这么道水题 没long long WA了两发,没判-1WA了一发,二分写错WA了一发 最 ...

  9. bzoj千题计划250:bzoj3670: [Noi2014]动物园

    http://www.lydsy.com/JudgeOnline/problem.php?id=3670 法一:KMP+st表 抽离nxt数组,构成一棵树 若nxt[i]=j,则i作为j的子节点 那么 ...

随机推荐

  1. docker之搭建私有仓库

    一.私有仓库 1.防止网络原因:下载慢,访问不到的情况,需要在内网搭建一个私有仓库. 二.仓库镜像下载 [root@node03 ~]# docker pull registry 三.创建私有仓库容器 ...

  2. 我的小游戏上线海外AppStore完整流程心得

    1,购买一台Mac或者用VMWare 安装Mac OS流程,笔者使用VMWare. 先安装Mac OS 10.13,教程,成功后不要着急安装vmtools, 首先更新系统至最新版,因为真机测试往往需要 ...

  3. PAT甲题题解-1030. Travel Plan (30)-最短路+输出路径

    模板题最短路+输出路径如果最短路不唯一,输出cost最小的 #include <iostream> #include <cstdio> #include <algorit ...

  4. 云平台项目--学习经验--BootstrapValidate表单验证插件

    使用前提,需要加载jquery和bootstrap库.并且引入bootstrapValidator.js和bootstrapValidator.css文件然后建立一个form表单,添加表单控件,表单控 ...

  5. Docker(四)-Dcoker镜像

    Docker 运行容器前需要本地存在对应的镜像,如果镜像不存在本地, Docker 会从镜像仓库下载(默认是Docker Hub公共注册服务器中的仓库). Docker Hub:https://hub ...

  6. [windows]转帖 windows 版本的含义

    1.N是None的意思,由于欧洲反垄断法不让系统捆绑浏览器ie和播放器.LTSB全称为Long Term Support Branch,意思为长期支持分支,是给Windows 10企业用户的一个更新选 ...

  7. angular浏览器滚动条滚动到指定element 触发事件

    angular.module('app').directive('ScrollTrigger', () => { return { restrict: "A", link:f ...

  8. python 按照固定长度分割字符串

    >>> import re >>> string = '123456789abcdefg' >>> re.findall(r'.{3}', str ...

  9. java通过dom读写xml文件

    java通过dom读写xml文件 要读的xml文件 <?xml version="1.0" encoding="GB2312"?><学生花名册 ...

  10. eclipse对maven项目进行打war包

    1. 利用pom.xml进行打包 右键项目名称或者pom.xml文件,选择debugAS或者runAS,选择maven install执行成功之后日志会打印输出位置 2. 使用命令打包 右键项目名称, ...