当前位置: 首页 > news >正文

CF512E. Cycling City

题目传送门

十分有趣的题。

思路

三条路径,本质上其实就是 \(x,y\) 同时属于两个有交集(至少交一条边)的简单环,这个肯定没问题。

套路的跑一遍 dfs,然后就有了返祖边树边和横叉边,然后朴素的分讨然后用个数据结构之类的维护一下就可以 \(n\log\) 解决了。

判有没有解倒是可以树上差分做,但这道题要求方案数,我们用一个暴力的思想,考虑每次对于一条非树边,我们暴力去覆盖树边,如果树边已经被覆盖,那么一定有解,取出这两条非树边然后去求答案即可。

具体求答案我简单叙述一下,具体还是自行画图看,首先把交集的点拿出来,然后按深度排序从小到大输出就完成了一部分。

然后对于 \(x,y\) 这条边,如果是横叉边,一定有一个点是结束点,然后输出就是形如 \(S,x,E\),其中 \(S,E\) 分别是开始点和结束点。否则是返祖边,那么假设 \(dep_x < dep_y\),那么就是 \(S,fa_S...,x,y,fa_y,...,E\),构造即可。

code

#include<bits/stdc++.h>
using namespace std;
#define int long long
namespace IO
{template<typename T>void read(T &_x){_x=0;int _f=1;char ch=getchar();while(!isdigit(ch)) _f=(ch=='-'?-1:_f),ch=getchar();while(isdigit(ch)) _x=_x*10+(ch^48),ch=getchar();_x*=_f;}template<typename T,typename... Args>void read(T &_x,Args&...others){Read(_x);Read(others...);}const int BUF=20000000;char buf[BUF],top,stk[32];int plen;#define pc(x) buf[plen++]=x#define flush(); fwrite(buf,1,plen,stdout),plen=0;template<typename T>inline void print(T x){if(!x){pc(48);return;}if(x<0) x=-x,pc('-');for(;x;x/=10) stk[++top]=48+x%10;while(top) pc(stk[top--]);}
}
using namespace IO;
const int N = 2e5+10;
int n,m,x,y,head[N],cnt,v[N],dep[N],son[N],op,fa[N],ly,ny,v1[N],v2[N],v3[N],len,len1,T[N],T1[N];
int X,Y,X1,Y1; 
struct w
{int to,nxt;
}b[N<<1];
inline void add(int x,int y)
{b[++cnt].nxt = head[x];b[cnt].to = y;head[x] = cnt;
}
void dfs(int x,int y)
{if(op == 1) return; v[x] = 1,dep[x] = dep[y]+1; fa[x] = y;for(int i = head[x];i;i = b[i].nxt){if(op == 1) return; if(!v[b[i].to]) dfs(b[i].to,x);else if(b[i].to != y)//不是树边{if(dep[b[i].to] == dep[x])//横叉边{ly = x; ny = b[i].to;while(ly != ny){//边转点 if(v1[ly]){X = x,Y = b[i].to;X1 = v1[ly],Y1 = v2[ly]; op = 1;break;}if(v1[ny]){X = x,Y = b[i].to;X1 = v1[ny],Y1 = v2[ny]; op = 1;break;}v1[ly] = x,v2[ly] = b[i].to;ly = fa[ly];v1[ny] = x,v2[ny] = b[i].to;ny = fa[ny];}}else if(dep[b[i].to] < dep[x])//返祖边{ly = x;while(ly != b[i].to){//边转点 if(v1[ly]){X = x,Y = b[i].to;X1 = v1[ly],Y1 = v2[ly]; op = 1;break;}v1[ly] = x,v2[ly] = b[i].to;ly = fa[ly];}}} }
}
inline void C(int x,int y)//染色,注意这里是染点 
{if(dep[x] == dep[y]){ly = x,ny = y;while(ly != ny) v1[ny] = x,v2[ny] = y,v3[ly]++,v3[ny]++,ly++,ny++; v3[ly]++;}else{ly = x;while(ly != y) v3[ly]++,ly = fa[ly];v3[ly]++;}
} 
inline bool cmp(int x,int y){ return dep[x] < dep[y]; }
inline void solve(int X,int Y)
{len1 = 0; ly = T[1]; T1[++len1] = ly;		if(dep[X] == dep[Y] && dep[Y] == ly) swap(X,Y);while(dep[Y] < dep[ly]) ly = fa[ly],T1[++len1] = ly;//S一直到Y去 if(dep[X] == dep[Y]) T1[++len1] = X,T1[++len1] = Y;//画图可以看出,这种一定是S->X->E这么连,dep相等一定有一个是Eelse {T1[++len1] = X;while(X != T[len]) X = fa[X],T1[++len1] = X;//一直走到E为止 }print(len1),pc(' '); for(int i = 1;i <= len1;i++) print(T1[i]),pc(' '); pc('\n');
}
signed main()
{//freopen(".in","r",stdin);//freopen(".out","w",stdout);read(n),read(m);for(int i = 1;i <= m;i++) read(x),read(y),add(x,y),add(y,x);//图不一定联通for(int i = 1;i <= n;i++)if(!v[i]){dfs(i,0);if(op){pc('Y'),pc('E'),pc('S'),pc('\n');C(X,Y),C(X1,Y1);// cout<<X<<" * "<<Y<<" "<<X1<<" "<<Y1<<" "<<dep[X]<<" "<<dep[Y]<<" "<<dep[X1]<<" "<<dep[Y1]<<'\n';for(int i = 1;i <= n;i++)if(v3[i] == 2) T[++len] = i;sort(T+1,T+1+len,cmp);print(len),pc(' '); for(int i = 1;i <= len;i++) print(T[i]),pc(' '); pc('\n');solve(X,Y); solve(X1,Y1);flush();return 0;}}pc('N'),pc('O'),pc('\n');flush(); return 0;
}
/*
10 20
5 1
10 5
2 10
6 4
10 6
9 6
1 7
3 10
3 2
6 2
5 4
7 10
3 9
9 1
5 3
7 9
8 4
4 7
6 3
10 8
注意到有三条路径
本质上就是有两个简单环交于同一条线段(交多部分也行)
考虑先跑dfs,这样只有返祖边,横叉边和树边
然后大致有四种情况
嗯,分别去看一下就好了,具体见代码
实现,emm目前我只会单log的,先写了来吧
我去我知道了,我们只需要看合不合法,然后构造任意一组方案
我们直接对于每条非树边暴力覆盖树边
然后有一条覆盖次数大于一那么一定有解了,输出
具体实现,请见代码 
*/
http://www.hskmm.com/?act=detail&tid=38421

相关文章:

  • 好想成为人类啊——2025 . 10 . 24
  • 10 24(+第14场补题)
  • 详细介绍:C++ 位运算 高频面试考点 力扣 268. 丢失的数字 题解 每日一题
  • 详细介绍:第十六届蓝桥杯软件赛C组省赛C++题解(京津冀)
  • OOP实验二
  • ABP - 缓存(Caching)[IDistributedCache、ICacheManager、ICacheKeyNormalizer、[Cache]、[CacheInvalidate]]
  • 《打造自己的 DeepSeek》第 1 期:为什么要打造自己的 DeepSeek?
  • ret2text
  • ABP - 异常处理(Exception Handling)[AbpExceptionFilter、UserFriendlyException、IExceptionSubscriber]
  • 2025年沸腾干燥机厂家权威推荐榜单:专业直销与高效节能技术深度解析,提供优质沸腾干燥设备及定制方案
  • CF Round 1046(#2135) 总结
  • 重组蛋白表达的几种类型介绍
  • ABP - 接口授权 [Authorize、AllowAnonymous、IPermissionChecker]
  • 日总结 17
  • Luogu P5479 [BJOI2015] 隐身术 题解 [ 紫 ] [ 多维 DP ] [ 交换维度 ] [ 后缀数组 ] [ 哈希 ]
  • 2025年10月23日
  • 杂题选做-3
  • 10.24每日总结
  • 利用Eval Villain挖掘CSPT漏洞的完整指南
  • Button按钮插入图片后仍有白色边框的解决办法
  • Hugo主题的修改和配置
  • 多元生成函数+多项式方程组——[AGC058D] Yet Another ABC String
  • ABP - JWT 鉴权(JWT Authentication)[AbpJwtBearerModule、JwtBearerOptions]
  • 最小生成树 kruskal算法
  • 【Java】Synchronized-你知道Java是如何上锁的吗?
  • Java中的字符串及相关类的介绍
  • ABP - 工作单元(Unit of Work)[UnitOfWorkAttribute、IUnitOfWorkManager、UnitOfWorkOptions]
  • LeetCode刷题笔记
  • [NOIP2023] 双序列拓展 题解
  • 洛谷 P9530 Fish 2