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+xn≡i+xs⋅d≡i(modd),i+ym≡i+yt⋅d≡i(modd)。

所以我们将所有男生 0∼n−10\sim n-10∼n−1 和女生 0∼m−10\sim m-10∼m−1 按编号取模 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,y,y−1y-1y−1 女生时间一定最大。

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

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

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

需要注意的是,初始就快乐的女生真正获得快乐的时间是 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;
}

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

  1. Codeforces 516D Drazil and Morning Exercise (栈、二分)

    题目链接 https://codeforces.com/contest/516/problem/D 题解 我还是数据结构水平太低了啊--连一个点子树内距离不超过\(l\)的点数都不会求 首先有一个熟知 ...

  2. codeforces 264 B. Good Sequences(dp+数学的一点思想)

    题目链接:http://codeforces.com/problemset/problem/264/B 题意:给出一个严格递增的一串数字,求最长的相邻两个数的gcd不为1的序列长度 其实这题可以考虑一 ...

  3. CodeForces - 1646E Power Board (思维,数学)

    题目链接:点击这里 题目大意: 给定一个 n×mn\times mn×m 的矩阵,其中 a[i][j]=ija[i][j]=i^ja[i][j]=ij ,求这个 n×mn\times mn×m 的矩阵 ...

  4. 枚举 ---- B. RPG Protagonist[Educational Codeforces Round 94 (Rated for Div. 2)]数学枚举

    B. RPG Protagonist 题目大意:就是你有两个人,每个人都有一个最大的体力值p和f,这两个人要去搬运剑和盾牌,剑的数量是cnts,盾的数量是cntw,每个剑的重量是是s,每个盾的重量是w ...

  5. Codeforces Educational 38 C. Constructing Tests ( 数学公式推导+暴力)

    C. Constructing Tests time limit per test 1 second memory limit per test 256 megabytes input standar ...

  6. codeforces 516c// Drazil and Park// Codeforces Round #292(Div. 1)

    题意:一个圆环上有树,猴子上下其中一棵树,再沿着换跑,再上下另一棵树.给出一个区间,问最大的运动距离是. 给出区间大小dst,和数高数组arr. 设区间[x,y],a[x]=2*arr[x]+dst[ ...

  7. CodeForces - 1330E Drazil Likes Heap(贪心)

    题目链接:点击查看 题目大意:首先定义 n 阶完全大顶堆的概念,就是满足节点个数恰好为个,序号从 1 开始排列,对于每个节点 i 都满足 接下来给出一段伪代码,教你如何从大顶堆中删去节点 i 最后给你 ...

  8. 【CodeForces - 1152C 】Neko does Maths(数学数论,lcm,gcd性质)

    题干: 给出a,b<=1e9,你要找到最小的k使得lcm(a+k,b+k)尽可能小,如果有多个k给出同样的最小公倍数,输出最小的一个k. 解题报告: 因为题目中k太多了,先化简一下公式,假设a& ...

  9. 【CodeForces - 1062C】Banh-mi (贪心,数学,找规律,快速幂)

    题干: JATC loves Banh-mi (a Vietnamese food). His affection for Banh-mi is so much that he always has ...

最新文章

  1. 多线程还是多进程的区别
  2. 粉丝该从哪里来?搞定这8种吸粉渠道
  3. 四、Delphi 2009 泛型容器单元(Generics.Collections)[3]: TStackT
  4. SDUT oj 选拔赛1 迷之好奇
  5. XDOC云服务API(二)
  6. 两个字符串的最长公共子序列长度_算法学习笔记(58): 最长公共子序列
  7. C# winform 弹出输入框
  8. MarkDown说明
  9. linux基础命令入门到精通
  10. 微信小程序从云开发到上线
  11. 飞鸽传书2011下载(飞鸽传书)
  12. 桌面版docker安装搭建
  13. URL中文参数错误解决方法
  14. 2021阿里淘系工程师推荐书单
  15. 辽宁移动_E900V21E_S905L3_线刷固件包
  16. word如何添加下滑线
  17. win7 IIS Web.config节点锁定问题
  18. Python Flask框架-开发简单博客-项目布局、应用设置
  19. js获取当月第一天年月日和最后一天年月日
  20. Android 控制台读取sqlite数据库以及adb获取root权限问题

热门文章

  1. 监督学习:KNN(K-近邻)算法实现手写数字识别的三种方法
  2. 【汇总推荐】深度学习、自然语言处理干货笔记汇总
  3. 看国外女神级程序员,直播写代码一年的感悟
  4. linux实时进程优先级rt,Linux实时性- PREEMPT_RT实时抢占实现
  5. android sdk中添加自定义api,android SDK中添加自定义api【转】
  6. 用java写注册表单_利用HTML表单标签编写一个注册页面
  7. 短文件名漏洞如何修复_短文件名漏洞修复
  8. 求职学习笔记|并发编程知识点常见问题总结(一)
  9. linux 7.0查看防火墙状态,centos7查看防火墙状态
  10. android photopicker怎么修改状态栏,有没有办法阻止UIImagePickerController更改状态栏样式?...