bzoj2503&poj3387[NEERC2006]IdealFrame
其实只是把别人的题解强行扩写了
写这篇题解之前我不会的预备知识:
欧拉通路:从图中一个点出发不重复地遍历所有边的路径(可以停在另一个点)
欧拉回路:从图中一个点出发不重复地遍历所有边的回路(必须回到出发点)
欧拉图:存在欧拉回路的图。判断无向图为欧拉图的充要条件是所有点的度数均为偶数。
半欧拉图:存在欧拉通路的图。判断无向图为欧拉图的充要条件是所有点的度数均为偶数或只有两个点的度数为奇数。
一个图中如果存在度数为奇数的点,那么这样的点一定有偶数个。
从一些简单的例子考虑到复杂的例子。
- 给定的图是一个环,答案为0
- 给定的图是一个欧拉图,那么我们只需要把其中的边按照欧拉回路展开即可,答案为度数大于2的点的个数。因为最后的环中所有点的度数均为2,所以度数大于2的点必须被熔断,而且存在只熔断度数大于2的点的方案使得欧拉回路展开成一个环。
例如左边”8”字形的图:从中间度数为4的点中间上下分开即可,得到右图,紫色节点为被熔断的节点形成的新节点。

3.给定的图是一个连通图但不是欧拉图。那么这个图中一定包含偶数个奇度点。因此我们可以将奇度点两两“配对”,每一对点之间连一条“边”,这样就把原图变成了一个所有点度数均为偶的欧拉图,可以将其展开成一个环。这个地方比较难以理解,因为原题没有加边操作,那么怎么连“边”呢?关键在于:虽然我们认为先“加边”再“展开”,但真正操作时,需要先进行熔断操作,再进行连接操作,因为一个度数大于1的点无法进行连接操作。因此连“边”对应到题中的操作,是将两个奇度点进行熔断后再将它们各自产生的一个度数为1的节点连接起来。计算答案的时候,首先统计出原图中奇度点数目和偶度点数目,需要花费奇度点数/2的花费将图转化成欧拉图。再统计出“连边”后图中有多少点的度数大于2。只需用总的点数减去最后度数为2的点数。注意原先度数为1 的点“连边”后会变成度数为2的点。
接下来考虑多个连通块。
4.给定的图不是一个连通图,但每个连通块都是一条链。首尾相接即可,答案=链数
5.给定的图不是一个连通图,但每个连通块都是一个环。每个环熔断一次拆成链再首尾相接,答案=链数×2
6. 给定的图不是一个连通图,每个连通块是一个欧拉图但不是环。我们是否需要先通过几次熔断把欧拉回路展开成环,然后再把环用一次熔断拆成链?不需要。事实上我们在把欧拉回路展开成环的过程中可以“顺手”把环拆成链。还是上面那个8字形的图,我们可以通过一步熔断完成“拆欧拉回路为环”和“拆环为链”。因此,如果一个欧拉回路不是环,那么在把它展开成环的过程中我们就可以完成“拆环为链”,不需要额外熔断一次。对于更加复杂的图,也是这样的。只有本来就是一个环的图才需要单独一次熔断“拆环为链”。

7.多个连通块,每个连通块没有特殊性质(即原题)。
我们希望将每个连通块都变成链,然后把每个链首尾相接即可。那么首先统计连接所有链的花费,把答案加上连通块数目。
然后我们算出每个连通块拆成链的最小费用。实质上我们需要把每个连通块变成半欧拉图(有一条欧拉通路的图)再展开。将半欧拉图展开成一条链的花费和将欧拉图展开成一个环的花费可以相似地计算,等于度数大于2的点数。接下来按照奇度点的个数讨论。
如果一个连通块的奇度点个数恰好为2,我们就可以直接计算展开成链的费用了。
如果奇度点个数大于2,我们需要用之前把“一些相关联的熔断和连接操作”等效为“加边”的思路,将奇度点的个数减少到2,并统计“加边”所需的连接次数(熔断次数不在这里统计)。加完“边”,就可以统计展开成链的费用了。
如果奇度点个数为0(是一个欧拉图),我们可以按照6计算这个连通块的费用。
然而考试写代码的时候并没有完全想清楚上面这些性质…..所以我写成了一大坨的分类讨论和手动枚举情况,还记录每个连通块中度数为1,2,3,4的点的个数,最后没有搞出“拆环为链”有时可以不用额外花费的性质,再加上一堆莫名其妙的细节错误,只有70...如果搞清楚上面这些性质,应该代码不会这样恶心。这个代码是考试代码基础上改的,不建议看,只是让读题解的你感受一下不想透性质就写代码会导致什么后果->_->
#include<cstdio>
const int maxn=,maxm=;
int u[maxm],v[maxm];
int ufs[maxn+maxm*];
int tot=;
int deg[maxn+maxm*];
int find(int x){
return x==ufs[x]?x:ufs[x]=find(ufs[x]);
}
void link(int a,int b){
ufs[find(a)]=find(b);
}
int odd[maxn+maxm*],even[maxn+maxm*],one[maxn+maxm*],two[maxn+maxm*],three[maxn+maxm*],four[maxn+maxm*];
int scc[maxn+maxm*];int cnt;
int maxdeg[maxn+maxm*];
int main(){
// freopen("frame.in","r",stdin);
int n,m;scanf("%d%d",&n,&m);
tot=n;
for(int i=;i<=m;++i){
scanf("%d%d",u+i,v+i);
if(u[i]==)u[i]=++tot;
if(v[i]==)v[i]=++tot;
deg[u[i]]++;deg[v[i]]++;
}
for(int i=;i<=tot;++i)ufs[i]=i;
for(int i=;i<=m;++i){
if(find(u[i])!=find(v[i]))link(u[i],v[i]);
}
for(int i=;i<=tot;++i){
if(deg[i]!=){
if(deg[i]&)odd[find(i)]++;
else even[find(i)]++;
if(deg[i]==)one[find(i)]++;
else if(deg[i]==)two[find(i)]++;
else if(deg[i]==)three[find(i)]++;
else if(deg[i]==)four[find(i)]++;
}
}
for(int i=;i<=tot;++i){
if(ufs[i]==i&°[i]!=){
scc[++cnt]=i;
}
}/*
for(int i=1;i<=tot;++i){
if(find(i)==scc[4])printf("%d\n",i);
}*/
// printf("%d\n",cnt);
if(cnt==){
int x=scc[];
two[x]+=one[x];
even[x]+=odd[x];
int ans=odd[x]/+(even[x]-two[x]);
printf("%d\n",ans);
}else{
int toteven=,tottwo=;
int ans=,x;
for(int i=;i<=cnt;++i){
x=scc[i];//printf("%4d %4d %4d %4d %4d %4d ",odd[x],even[x],one[x],two[x],three[x],four[x]);
if(odd[x]==){
if(four[x]!=){
two[x]++;odd[x]=one[x]=;
four[x]--;
}
else if(two[x]!=){
two[x]--;even[x]--;
one[x]=odd[x]=;
}else{
odd[x]=one[x]=;
}
ans++;
tottwo+=two[x];
toteven+=even[x];
}else{//竟然还有可能拆偶点...
if(one[x]>=){//不需要再拆1了
odd[x]-=;
one[x]-=;
two[x]+=one[x];
ans+=odd[x]/;
}else if(one[x]==){//再拆一个1
if(three[x]!=){
three[x]--;ans++;
two[x]++;
odd[x]-=;
even[x]++;
ans+=odd[x]/;
}else{//这里是拆奇点,能否拆偶点?
odd[x]-=;
even[x]++;
ans++;
ans+=odd[x]/;
}
}else if(one[x]==){//拆出两个1
int g1=0x7f7f7f7f,g2=0x7f7f7f7f,g3=0x7f7f7f7f,g4=0x7f7f7f7f,g5=0x7f7f7f7f;
if(three[x]>=){//拆两个3
g1=+(odd[x]-)/;
}
if(two[x]!=){
g2=+odd[x]/;
}
if(four[x]!=){
g3=odd[x]/-;
}
if(three[x]!=odd[x]&&three[x]!=){
g4=+(odd[x]-)/;
}
if(three[x]==){
g5=+(odd[x]-)/;
}
if(g1<=g2&&g1<=g3&&g1<=g4&&g1<=g5){
ans+=;
odd[x]-=;even[x]+=;two[x]+=;
ans+=(odd[x]/);
}else if(g2<=g1&&g2<=g3&&g2<=g4&&g2<=g5){
ans+=;
two[x]--;even[x]--;
ans+=odd[x]/;
}else if(g3<=g1&&g3<=g2&&g3<=g4&&g3<=g5){
two[x]++;//ans+=1;
ans+=odd[x]/;
}else if(g4<=g1&&g4<=g2&&g4<=g3&&g4<=g5){
two[x]++;even[x]+=;//ans+=2;
odd[x]-=;
ans+=odd[x]/;
}else{//printf("!");
even[x]+=;//ans+=2;
odd[x]-=;
ans+=odd[x]/;
}
}
tottwo+=two[x];
toteven+=even[x]+odd[x];
}//printf("ansnow %d\n",ans);
}
ans+=(toteven-tottwo);ans+=cnt;
// printf("%d %d\n",toteven,tottwo);
printf("%d\n",ans);
}//while(1);
return ;
}
bzoj2503&poj3387[NEERC2006]IdealFrame的更多相关文章
- 【BZOJ-2503】相框 并查集 + 分类讨论
2503: 相框 Time Limit: 3 Sec Memory Limit: 128 MBSubmit: 71 Solved: 31[Submit][Status][Discuss] Desc ...
- [BZOJ2503][HAOI2006]均分数据
BZOJ Luogu sol 如果已经确定了一个序列,现要求把这个序列分成m个连续段作为答案,那么就可以用一个显而易见的DP DP显然可以得到当前序列下的最优解. 所以模拟退火瞎JB改一改序列每次DP ...
- BZOJ.1312.[Neerc2006]Hard Life(分数规划 最大权闭合子图)
BZOJ 最大密度子图. 二分答案\(x\),转为求是否存在方案满足:\(边数-x*点数\geq 0\). 选一条边就必须选两个点,所以可以转成最大权闭合子图.边有\(1\)的正权,点有\(x\)的负 ...
- Bzoj1312 / POJ3155 Neerc2006 Hard Life
Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 459 Solved: 114 Description 在一家公司中,人事部经理与业务部经理不和.一次 ...
- bzoj2503 相框——思路
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=2503 思路题: 首先,这种问题应该注意到答案只跟度数有关,跟其他什么连接方法之类的完全无关: ...
- bzoj AC倒序
Search GO 说明:输入题号直接进入相应题目,如需搜索含数字的题目,请在关键词前加单引号 Problem ID Title Source AC Submit Y 1000 A+B Problem ...
随机推荐
- The operation couldn’t be completed. (LaunchServicesError error 0.)
问题描述: 当运行Xcode时,编译代码成功,但是登陆模拟器失败,显示错误:The Operation couldn't be completed.(LaunchServicesError error ...
- YourSQLDba将数据库置于紧急模式的原因浅析
从SQL SERVER 2000 上迁移了一个数据库到SQL SERVER 2008 R2上,暂且用DataBaseName代替迁移的真实的数据库名(后面的资料也会将数据库真实的名字用DataBase ...
- Ant:build.xml 结构
Ant build.xml 结构 project target task data property datatype v\:* {behavior:url(#default#VML);} o\:* ...
- jQuery 中的 39 个技巧
1.当document文档就绪时执行JavaScript代码. 我们为什么使用jQuery库呢?原因之一就在于我们可以使jQuery代码在各种不同的浏览器和存在bug的浏览器上完美运行. <sc ...
- W3School-CSS 尺寸 (Dimension) 实例
CSS 尺寸 (Dimension) 实例 CSS 实例 CSS 背景实例 CSS 文本实例 CSS 字体(font)实例 CSS 边框(border)实例 CSS 外边距 (margin) 实例 C ...
- 【转】关于Java的Daemon线程的理解
原文地址:http://www.cnblogs.com/ChrisWang/archive/2009/11/28/1612815.html 关于Java的Daemon线程的理解 网上对Java的Dae ...
- IP分片重组的分析和常见碎片攻击 v0.2
IP分片重组的分析和常见碎片攻击 v0.2http://www.nsfocus.net/index.php?act=magazine&do=view&mid=584 作者:yawl ( ...
- strstr 函数的实现
strstr函数:返回主串中子字符串的位置后的所有字符. #include <stdio.h> const char *my_strstr(const char *str, const c ...
- Centos下yum安装PHP
添加 yum 源 Centos 6.x 的源 # rpm -Uvh http://download.fedoraproject.org/pub/epel/6/x86_64/epel-release-6 ...
- [原]经典bootstrap模态框使用文章
1,Bootstrap 模态对话框和简单使用 <div id="myModal" class="modal hide fade"> <div ...