标签：frac gcd int 题解 bmod P7835 times Wdoi 1005

题目的意思在题目中有化简版。

首先看到数据，这么大直接放弃吧发现可以从 \(k\) 入手，可以做到 \(O(k^2)\)。

我们就可以想到枚举每两个 \(x_i,y_i\) 和 \(x_j,y_j\)，可能在同一个时刻第一次必然是 \(L=\text{lcm}(t_i , t_j)\)。这个时候，两个球的编号是 \(x_i + ( \frac{L}{t_i} - 1 ) \times y_i \bmod n\) 和 \(x_j + ( \frac{L}{t_j} - 1 ) \times y_j \bmod n\)。如果他们是不一样的，我们把所有这种可能的 \(L\) 取一个最小值。

如果他们一样，可能正好是巧合，下次就不一定了。所以我们考虑 $ x_i + ( 2 \times \frac{L}{t_i} - 1 ) \times y_i \bmod n$ 和 $ x_j + ( 2 \times \frac{L}{t_j} - 1 ) \times y_j \bmod n$。如果他们不相等，\(2 \times L\) 也要取最小值。

如果还不相等，就是永远不会相等了。赛场上我以防万一写了个3倍，不过没啥用（

最后，如果没有找到符合的情况，输出 Mystia will cook forever...。

有一些微不足道的优化，好像没啥用。

代码：

#include<iostream>
#define int unsigned long long
using namespace std;
int t[1005],x[1005],y[1005],at;//at是保存最小公倍数的
int gcd(int a,int b){
	return b==0?a:gcd(b,a%b);
}
int lcm(int a,int b){//基本数论常识
	return a/gcd(a,b)*b;
}
signed main(){
	int n,k,ans=-1;//ans是最小值，不可能会有-1的时间
	cin>>n>>k;
	if(k==1){//可以不用判断
		cout<<"Mystia will cook forever..."<<endl;
		return 0;
	}
	for(int i=1;i<=k;i++){
		cin>>t[i]>>x[i]>>y[i];
		x[i]%=n,y[i]%=n;//小优化
		for(int j=1;j<i;j++){//寻找每一对i，j
			at=lcm(t[i],t[j]);//求出lcm
			if((x[i]+(at/t[i]-1)%n*y[i])%n!=(x[j]+(at/t[j]-1)%n*y[j])%n) ans=min(ans,at-1);//注意，问的是存活时间，at是死亡时间
			else if((x[i]+(at/t[i]*2-1)%n*y[i])%n!=(x[j]+(at/t[j]*2-1)%n*y[j])%n) ans=min(ans,2*at-1);//同理
			else if((x[i]+(at/t[i]*3-1)%n*y[i])%n!=(x[j]+(at/t[j]*3-1)%n*y[j])%n) ans=min(ans,3*at-1);//好像可以删掉
		}
	}
	if(ans==-1) cout<<"Mystia will cook forever..."<<endl;//没有找到
	else cout<<ans<<endl;
	return 0;
}

标签：frac,gcd,int,题解,bmod,P7835,times,Wdoi,1005
来源： https://www.cnblogs.com/tigerchen/p/15868872.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。