标签：le 洛谷 int 荒岛 color 野人 P2421 xp 山洞
前言
\(\color{red}{\texttt{WYCdalao}}\)在\(1.542\)年之前切的题现在才写掉\(qwq\)，\(\color{grey}{\texttt{stoorz}}\)怕不是要被\(\color{red}{\texttt{WYCdalao}}\)碾爆\(qwqwq\)。
\[\color{white}{\texttt{上图不可以点QWQ}}\]
题目
克里特岛以野人群居而著称。岛上有排列成环行的M个山洞。这些山洞顺时针编号为1,2,…,M。岛上住着N个野人，一开始依次住在山洞C1,C2,…,CN中，以后每年，第i个野人会沿顺时针向前走Pi个洞住下来。
每个野人i有一个寿命值Li，即生存的年数。
下面四幅图描述了一个有6个山洞，住有三个野人的岛上前四年的情况。三个野人初始的洞穴编号依次为1，2，3；每年要走过的洞穴数依次为3，7，2；寿命值依次为4，3，1。
奇怪的是，虽然野人有很多，但没有任何两个野人在有生之年处在同一个山洞中，使得小岛一直保持和平与宁静，这让科学家们很是惊奇。他们想知道，至少有多少个山洞，才能维持岛上的和平呢？
思路
这道题的实际数据范围是\(1\le n\le 15\)，\(1\le C_i,P_i\le 100\)，\(0\le L_i\le 10^6\)。
其实题目说了那么多，题意其实就是给出\(c_i,p_i,l_i\)，求一个最小的\(m\)使得所有方程组
\[\left\{\begin{matrix}c_i+xp_i\equiv c_j+xp_j\ (\rm{mod}\ m)
\\ x\leq l_i
\\ x\leq l_j
\end{matrix}\right.\]
无解。
由于\(m\)较小而且题目保证一定存在一个\(m\leq 10^6\)的解。所以考虑枚举\(m\)，然后枚举每一个野人。
\[c_i+xp_i\equiv c_j+xp_j\ (mod\ m) \Leftrightarrow (p_i-p_j)x+m·y=x_j-c_i\]
我们发现可以用\(Exgcd\)来求出满足要求的\(x\)。
那么要使得上述方程组无解，必须满足一下条件之一：
- \(c_j-c_i\)不是\(gcd(p_i-p_j,m)\)的倍数
- 求出来的\(x>min(l_i,l_j)\)
时间复杂度\(O(mn^2\log p_i)\)。
代码
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int N=20;
int n,m,c[N],p[N],l[N];
bool flag;
int exgcd(int a,int b,int &x,int &y)
{
if (!b)
{
x=1; y=0;
return a;
}
int d=exgcd(b,a%b,x,y);
int temp=x; x=y; y=temp-a/b*y;
return d;
}
void check(int i,int j)
{
if (p[i]<p[j]) swap(i,j);
int x,y,d=exgcd(p[i]-p[j],m,x,y);
if ((c[j]-c[i])%d) return;
x=(x*(c[j]-c[i])/d%(m/d)+(m/d))%(m/d);
if (x<=min(l[i],l[j])) flag=0;
}
int main()
{
scanf("%d",&n);
for (int i=1;i<=n;i++)
{
scanf("%d%d%d",&c[i],&p[i],&l[i]);
m=max(m,c[i]);
}
for (;;m++)
{
flag=1;
for (int i=1;i<=n && flag;i++)
for (int j=i+1;j<=n && flag;j++)
check(i,j);
if (flag) return !printf("%d",m);
}
printf("WYC AK IOI OrzOrz");
return 0;
}标签：le,洛谷,int,荒岛,color,野人,P2421,xp,山洞 来源： https://www.cnblogs.com/stoorz/p/12230804.html
本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享； 2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关； 3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关； 4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除； 5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。