CodeForces 516E Drazil and His Happy Friends(数学+最短路)

problem

洛谷链接

solution

gcd⁡(n,m)=d\gcd(n,m)=dgcd(n,m)=d,则快乐只会在同余 ddd 的关系中传递。

i+xn≡i+xs⋅d≡i(modd),i+ym≡i+yt⋅d≡i(modd)i+xn\equiv i+xs·d\equiv i\pmod d,i+ym\equiv i+yt·d\equiv i\pmod di+xni+xsdi(modd),i+ymi+ytdi(modd)

所以我们将所有男生 0∼n−10\sim n-10n1 和女生 0∼m−10\sim m-10m1 按编号取模 ddd 的结果分组 [0,d)[0,d)[0,d)

因为每一类的快乐只能在组内传递,所以无解的情况就是某个组内没有一个人初始是快乐的。

答案显然是取所有人获得快乐的最早时间的最大值。

我们分组计算答案,考虑任意一组内的答案如何计算?

假设这里面有个男生 iii 一开始就是快乐的,那么他会在第 iii 天将快乐传递给 imodmi\mod mimodm 女生。

再过 nnn 天,这个男生又把自己的快乐传递给了 (i+n)modm(i+n)\mod m(i+n)modm 女生。

那我们能否看成是 imodmi\mod mimodm 女生经过 nnn 天将快乐传递给 (i+n)modm(i+n)\mod m(i+n)modm 女生的呢?——答案是肯定的。

于是我们就又有了将同组的男女生分开计算的想法。

不妨以女孩子为例。(女孩纸就是香香的),将初始快乐的男生 iii 信息附属在女生 imodmi\mod mimodm 身上,即把她当成初始就快乐的人。

考虑如何组内每个女孩纸最早获得快乐的时间?

我们将这个快乐传递当成边,时间就是边权。你会发现每条边的边权都是 nnn

跑个最短路就行了?——太遗憾了,点可能会很多。

我们关注到初始快乐的女生总数不会超过 1e51e51e5,很容易想到能否通过这些关键女生计算出每个女生获得快乐的时间?

在这里还获得了一个粗略判无解的条件,那就是分的组数 d>2e5d>2e5d>2e5,每组只放个男生或女生都不能让每组里有初始即快乐的人。

将本组所有女生编号处理出来,并将边连出来,肯定能刚好连出一个环来。

让所有女生都多加一个属性——环的重编号。

则任意两个女生快乐传递时间就是环上编号差 ×n\times n×n(注意环头和环尾的特殊计算方式)。

现在将组内所有女生按环的重编号递增排列,则相邻两个女生间快乐传递时间均为 nnn

继续考虑重新排列后,连续两个关键女生中夹着的一些普通女生。

disi:idis_i:idisi:i 最早获得快乐的时间,且对于所有关键女生 kkk 初始化 disk=kdis_k=kdisk=k

我们没有必要去算这中间每个女生的时间,因为最后答案只取最大值。

所有这段女生中时间最大值肯定是编号最大的,即重编号的连续两个关键女生 x,yx,yx,yy−1y-1y1 女生时间一定最大。

那我们只用算这一段的这一个特殊的普通女生时间即可。

然后将每一段的答案取较大值。

线性扫一遍关键女生即可求解。

需要注意的是,初始就快乐的女生真正获得快乐的时间是 000 而不是 disidis_idisi。(所以不能选这些女生的时间比最大值)

还有我们把某些男生快乐的信息放在了某些女生身上,让她们成为关键女生,但这是虚假的关键女生,所以这些女生的 disidis_idisi 又是需要计算的。

男生同理,交换一下 n,mn,mn,m 等信息即可。

最后在每组的女生/男生最大值中再取最大值。

code

#include <bits/stdc++.h>
using namespace std;
#define inf 1e18
#define int long long
#define maxn 200005
vector < int > b[maxn], g[maxn];
int n, m, B, G;
int d[maxn], id[maxn], s[maxn];int exgcd( int a, int b, int &x, int &y ) {if( ! b ) { x = 1, y = 0; return a; }else { int d = exgcd( b, a % b, y, x  ); y -= a / b * x; return d; }
}int solve( int gcd, int n, int m, vector < int > g, vector < int > b ) {//每组的男女生个数就是 n/gcd,m/gcdif( m == g.size() ) return -1; //要是return 0外层还要+i有可能还是被误判成最大值int cnt = 0, ans = 0, dis = inf;/*对于每一组,我们怎么将男生 $i\mod n$ 的快乐移加到女生 $i\mod m$ 的快乐上呢?我们按 $i\mod d$ 分组的时候,扔的关键男女生编号就直接是 $\lfloor\frac id\rfloor$。那么将所有组内编号再都 $\times d\pmod m$ 取出来。相当于是把模 $d$ 相同的数移到 $d$ 的整数倍,$r+id\rightarrow id$。然后再一一对应到模 $m$ 的位置上。离散化的赶脚 因为编号太大不可能直接以编号为下标做*/for( int i : g ) s[++ cnt] = i * gcd % m, d[cnt] = i, id[cnt] = cnt;//真的关键女生点for( int i : b ) s[++ cnt] = i * gcd % m, d[cnt] = i, id[cnt] = cnt;//假的关键女生点sort( id + 1, id + cnt + 1, []( int x, int y ) { return s[x] < s[y]; } );s[id[0] = 0] = s[id[cnt]] - m, s[id[cnt + 1] = cnt + 1] = s[id[1]] + m;//把环拆成链for( int i = 1;i <= cnt;i ++ )dis = min( d[id[i]], dis + ( s[id[i]] - s[id[i - 1]] ) * n );for( int i = 1;i <= cnt;i ++ ) {dis = min( d[id[i]], dis + ( s[id[i]] - s[id[i - 1]] ) * n );if( s[id[i]] == s[id[i + 1]] ) continue;if( s[id[i]] + 1 == s[id[i + 1]] and id[i] <= g.size() ) continue;//不是虚假特殊女生才能跳过ans = max( ans, dis + ( s[id[i + 1]] - s[id[i]] - 1 ) * n );}return ans;
}signed main() {scanf( "%lld %lld", &n, &m );int x, y, w, gcd = exgcd( n, m, x, y );if( gcd > 2e5 ) return ! puts("-1");scanf( "%lld", &B );for( int i = 1;i <= B;i ++ ) scanf( "%lld", &w ), b[w % gcd].push_back( w / gcd );scanf( "%lld", &G );for( int i = 1;i <= G;i ++ ) scanf( "%lld", &w ), g[w % gcd].push_back( w / gcd );int ans = 0;( x += m ) %= m, ( y += n ) %= n;for( int i = 0;i < gcd;i ++ ) {//第i天才会开始第i组信息的传递if( g[i].empty() and b[i].empty() ) { puts("-1"); return 0; }ans = max( ans, solve( x, n / gcd, m / gcd, g[i], b[i] ) * gcd + i );ans = max( ans, solve( y, m / gcd, n / gcd, b[i], g[i] ) * gcd + i );}printf( "%lld\n", ans );return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/316496.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

.NetCore使用skywalking实现实时性能监控

一、简介很久之前写了一篇 《.Net Core 2.0 InfluxDBGrafanaApp Metrics 实现跨平台的实时性能监控》关于NetCore性能监控的文章&#xff0c;使用InfluxdbAppMetrics进行项目性能监控&#xff0c;由于技术有限&#xff0c;在正式环境使用一段时间后&#xff0c;莫名的AppMetric…

Hamburger Steak

Hamburger Steak 题意&#xff1a; 有m个锅&#xff0c;n块肉&#xff0c;每个肉需要煎ai分钟才能熟&#xff0c;每个锅同时最多煎一块肉&#xff0c;一块肉最多可以被两个锅煎(但不能同时)&#xff0c;问最少多长时间能全部煎熟&#xff0c;并输出方案 题解&#xff1a; 我…

「ROI 2017 Day 2」反物质(单调队列优化dp)

problem LOJ2772 solution 这道题求的是“保证最多”&#xff0c;这个“保证”真的屑啊&#xff01; 因为我们无法确定实验生成的克数&#xff0c;所以我们应当计算的是最坏情况。 又因为我们可以选择做哪些实验&#xff0c;所以要的是所有实验组合最坏情况下的最大价值。…

netcore开发windows普通服务(非Web)并一键发布到服务器

netcore下开发windows服务如果是web项目的话&#xff0c;由于aspnetcore本身是支持的&#xff0c;把默认的host.Run改为host.RunAsService就可以了。但是普通的netcore的控制台项目我终于找到了如下方式来实现&#xff1a;Microsoft.Extensions.HostingSystem.ServiceProcess.S…

Hopping Rabbit

Hopping Rabbit 题意&#xff1a; 给你n个矩阵&#xff0c;每个矩阵(给出左上标和右下标)&#xff0c;现在让你给出一个点的位置&#xff0c;这个点每次只能上下左右四个方向移动&#xff0c;且移动距离为d&#xff0c;是否存在一个这样的点&#xff0c;其所有落点都不在矩阵…

删好串(区间dp)

problem 定义长度为 nnn 的“好”的串 aaa 满足&#xff1a; ∣ai−ai−1∣1,i∈[2,n]|a_i-a_{i-1}|1,i\in[2,n]∣ai​−ai−1​∣1,i∈[2,n]。ai≥ai−1ai12,i∈[2,n−1]a_i\ge \frac{a_{i-1}a_{i1}}{2},i\in[2,n-1]ai​≥2ai−1​ai1​​,i∈[2,n−1]。 给定长度为 nnn 的序…

直播预告 - 微软MVP为你揭秘Visual Studio 2019新特性

作为"宇宙第一IDE“的Visual Studio集成开发环境&#xff0c;已经经历了超过十几年的迭代成为一款功能丰富且高效的开发工具&#xff0c;微软自己给Visual Studio 的定位是 “更快、更可靠&#xff0c;对个人和团队更具生产力&#xff0c;更易于使用&#xff0c;并且更容易…

Delete Edges

Delete Edges 题意&#xff1a; 给出一个n个点的完全图&#xff0c;删一些三元环使得边数<n 输出所删的三元环 3<n<2000 题解&#xff1a; 肯定是结论题&#xff0c;但是我不会。。 结论&#xff1a;xyz0(mod n)1<x<y<z<n所有解即可 证明过程可以看看…

.NET 机器学习生态调查

机器学习是一种允许计算机使用现有数据预测未来行为、结果和趋势的数据科学方法。 使用机器学习&#xff0c;计算机可以在未显式编程的情况下进行学习。机器学习的预测可以使得应用和设备更智能。 在线购物时&#xff0c;机器学习基于历史购买推荐你可能喜欢的其他产品。 刷信用…

图中异色点对最短距离(最小生成树+线段树)

problem 给定一个 nnn 个点&#xff0c;mmm 条边的无向连通图&#xff0c;图有边权&#xff0c;每个点有一个颜色。 有 qqq 次操作&#xff0c;每次操作可更改某一个点颜色。 求每次操作后图中不同颜色点之间的最短距离。 若图中点颜色全相同&#xff0c;输出 000。 一行给…

Defend Your Country

Defend Your Country 题意&#xff1a; n个点&#xff0c;m条边的简单无向连通图&#xff0c;每个点一个权值ai,一个连通块的贡献&#xff1a;(−1)块内点数∗∑ai[点i在该连通块内](-1)^{块内点数}*\sum a_{i}[点i在该连通块内](−1)块内点数∗∑ai​[点i在该连通块内] 可以…

SkyWalking Liunx 环境搭建NetCore接入

背景前两天看见有小哥介绍windows下安装skywalking的介绍地址。正好最近也在搭建linux环境的SkyWalking&#xff0c;顺便把linux环境搭建的经验分享下&#xff0c;帮助下使用linux部署DotNetCore项目的同学。介绍SkyWalking是开源的apm工具&#xff0c;服务器端使用java编写&am…

AtCoder 2305 [AGC010D] Decrementing(博弈)

problem Luogu链接 solution case1&#xff1a;aaa 中如果有 111 的存在&#xff0c;那么相当于第二步失效&#xff0c;最后的结果仅由 ∑ai−1\sum a_i-1∑ai​−1 的奇偶性决定&#xff08;奇数先手赢&#xff0c;偶数后手赢&#xff09; 奇数先手赢&#xff0c;偶数后手…

P4568 [JLOI2011]飞行路线

P4568 [JLOI2011]飞行路线 题意&#xff1a; n个城市&#xff0c;m个航班&#xff0c;你有k次免费坐航班的机会&#xff0c;问从s到t最少花费是多少&#xff1f; 题解&#xff1a; 分层图问题 简单说说什么是分层图&#xff1a; 其实就是将一个平面的图重新建图&#xff0c…

.NET 基金会完成第一次全面改选

.NET基金会是一个独立的组织&#xff0c;支持.NET社区和开源&#xff0c;旨在拓宽和加强.NET生态系统和社区。这可以通过多种方式完成&#xff0c;包括项目指导&#xff0c;指导&#xff0c;法律和营销帮助&#xff0c;技术和财务支持设置等&#xff0c;2014年微软组织成立.NET…

一般图带权多重匹配(欧拉图+最小费用流)

problem 给定 nnn 个数 {ai}\{a_i\}{ai​}&#xff0c;其中 kkk 个 aia_iai​ 是奇数&#xff0c;再给一个 nnn\times nnn 的矩阵 {ci,j}\{c_{i,j}\}{ci,j​}&#xff0c;无论是 aaa 还是 ccc&#xff0c;都保证是非负整数。 现在可以做一类操作&#xff1a;将 ai−1,aj−1a_…

xay loves trees

xay loves trees 题意&#xff1a; 有两棵树&#xff0c;现在让你找到一个最大的点集合S&#xff0c;要求S中的点在第一棵树中任意两点存在祖先儿子关系且所有点是连接的&#xff0c;在第二棵树中任意两点都不存在祖先儿子关系&#xff0c;问S集合的最大是多少&#xff1f; …

特来电混沌工程实践

一、导语随着大型分布式系统架构的演进和广泛应用&#xff0c;软件工程的最佳实践也随之改变。我们通过分布式、服务化、DevOps、敏捷开发&#xff0c;快速响应业务的需求变化&#xff0c;支持大规模分布式应用。但这些做法带来效益的同时&#xff0c;也带来了另一个紧迫问题&a…

[ZJOI2007] 棋盘制作(单调栈 / DP悬线法)

problem 洛谷链接 solution1-单调栈 很容易想到&#xff0c;预处理出每个点向上最大能延伸的长度&#xff0c;然后对每个点求一个矩阵面积。 然后思考优化&#xff0c;不难想到每次对一行进行求解。 每一行的所有列一起构成了一个直方图。 直方图直接经典笛卡尔树。笛卡尔…

P2567 [SCOI2010]幸运数字

P2567 [SCOI2010]幸运数字 题意&#xff1a; 我们规定只含6或8的数字为幸运号码&#xff0c;而幸运号码的倍数我们也认为是幸运号码&#xff0c;问[l,r]中有多少个幸运号码&#xff1f; 题解&#xff1a; 第一反应是数位dp&#xff0c;但其实不是&#xff0c;我们可以打表观…