I. Rise of Shadows

本文介绍了一种算法,通过计算时针和分针的速度差,将时针问题转化为追及相遇问题,利用同余性质简化求解,以求得在H小时M分钟内,时针与分针角度小于等于给定角度α的整数次。关键步骤包括找到最大公约数并应用到时间计算中。
摘要由CSDN通过智能技术生成

LINK
题意:给定H和M,每天有H小时M分钟,求时针和分针角度<= α \alpha α的整分针数。
思路:
α = 2 π A H M \alpha=\frac{2\pi A}{HM} α=HM2πA,可以发现每分钟时针的速率是 v 时 = 2 π H M v_{时}=\frac{2\pi }{HM} v=HM2π,分针是 v 分 = 2 π M v_{分}=\frac{2\pi }{ M} v=M2π,他俩的差值是 2 π ( H − 1 ) H M \frac{2\pi(H-1)}{HM} HM2π(H1)我们会发现两个指针一块动的时候不太好计算,我们尝试固定时针,那么现在就变成了追及相遇问题,也就是
t ∗ ( H − 1 ) % ( H M ) < = ∣ A ∣ − A < = t ∗ ( H − 1 ) % ( H M ) < = A    t ∈ [ 0 , H ∗ M − 1 ] 两 边 同 时 除 以 g − A g < = t ∗ H − 1 g < = A g   % ( H M g ) 由 于 同 余 的 性 质 , a ∗ c ≡ b ∗ c ( % m ) 等 价 于 a ≡ b ( % ( m g c d ( c , m ) ) 我 们 先 算 正 数 的 , 算 出 正 数 以 后 负 数 是 对 称 的 t < = A g ( % ( H ∗ M / g ) )    g = gcd ⁡ ( H M , H − 1 )   t ∈ [ 0 , H ∗ M / g − 1 ] 那 么 这 个 数 量 就 是 2 ∗ A g + 1 , 这 里 的 1 就 是 0 点 的 贡 献 t*(H-1) \% (HM)<=|A|\\ -A<=t*(H-1)\%(HM)<=A ~~t\in[0,H*M-1]\\ 两边同时除以g\\ -\frac{A}{g}<=t*\frac{H-1}{g}<=\frac{A}{g}~\%(\frac{HM}{g})\\ 由于同余的性质,a*c\equiv b*c(\% m) \\等价于a\equiv b(\%(\frac{m}{gcd(c,m)})\\ 我们先算正数的,算出正数以后负数是对称的\\ t<=\frac{A}{g}(\%(H*M/g))~~g=\gcd(HM,H-1) ~t\in [0,H*M/g-1]\\ 那么这个数量就是2*\frac{A}{g}+1,这里的1就是0点的贡献 t(H1)%(HM)<=AA<=t(H1)%(HM)<=A  t[0,HM1]ggA<=tgH1<=gA %(gHM)acbc(%m)ab(%(gcd(c,m)m)t<=gA(%(HM/g))  g=gcd(HM,H1) t[0,HM/g1]2gA+1,10
AC代码:

//#pragma GCC target("avx")
//#pragma GCC optimize(2)
//#pragma GCC optimize(3)
//#pragma GCC optimize("Ofast")
// created by myq
#include<iostream>
#include<cstdlib>
#include<string>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<climits>
#include<cmath>
#include<cctype>
#include<stack>
#include<queue>
#include<list>
#include<vector>
#include<set>
#include<map>
#include<sstream>
#include<unordered_map>
#include<unordered_set>
using namespace std;
typedef long long ll;
#define x first
#define y second
typedef pair<int,int> pii;
const int N = 400010;
const int mod=998244353;
inline int read()
{
	int res=0;
	int f=1;
	char c=getchar();
	while(c>'9' ||c<'0')
	{
		if(c=='-')	f=-1;
		c=getchar();
	}
	while(c>='0'&&c<='9')
	{
		res=(res<<3)+(res<<1)+c-'0';
	}
	return res;
 }
#define int long long
signed main()
{
	int h,m;
	cin>>h>>m;
	int a;
		cin>>a;
	int g=__gcd(h-1,h*m);
	cout<<min(h*m,g*(1+2*(a/g)));
	return 0;

}
/**
* In every life we have some trouble
* When you worry you make it double
* Don't worry,be happy.
**/



  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值