Codeforces 492E Vanya and Field(拓展欧几里得)

题目链接:Codeforces 492E Vanya and Field

通过拓展欧几里得算法求出每个位置在移动过程中,在x为0时,y的位置。统计相应y坐标最多的即为答案。

#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;
const int maxn = 1e6 + 5;
const int maxm = 1e5 + 5;
typedef long long ll;

int N, M, DX, DY, C[maxn], X[maxm], Y[maxm], V[maxm];

inline void exgcd (int a, int b, int& d, int& x, int& y) {
	if (!b) {
		d = a; x = 1; y = 0;
	} else {
		exgcd(b, a%b, d, y, x);
		y -= x * (a / b);
	}
}

int main () {
	scanf("%d%d%d%d", &N, &M, &DX, &DY);

	int x, y, t;
	exgcd(N, -DX, t, x, y);
	ll T = (y * t % N + N) % N;

	for (int i = 0; i < M; i++) {
		scanf("%d%d", &X[i], &Y[i]);
		V[i] = (Y[i] + T * X[i] % N * DY % N + N) % N;
		C[V[i]]++;
	}

	int ans = 0;
	for (int i = 0; i < N; i++)
		if (C[i] > ans) {
			ans = C[i];
			y = i;
		}

	for (int i = 0; i < M; i++)
		if (V[i] == y) {
			printf("%d %d\n", X[i], Y[i]);
			break;
		}

	return 0;
}
时间: 2024-10-12 14:34:20

Codeforces 492E Vanya and Field(拓展欧几里得)的相关文章

codeforces 492E. Vanya and Field(exgcd求逆元)

题目链接:codeforces 492e vanya and field 留个扩展gcd求逆元的板子. 设i,j为每颗苹果树的位置,因为gcd(n,dx) = 1,gcd(n,dy) = 1,所以当走了n步后,x从0~n-1,y从0~n-1都访问过,但x,y不相同. 所以,x肯定要经过0点,所以我只需要求y点就可以了. i,j为每颗苹果树的位置,设在经过了a步后,i到达了0,j到达了M. 则有 1----------------------(i + b * dx) % n = 0 2------

CodeForces 492E Vanya and Field

题意: n*n的矩形内有m(10^5)个点  从任一点开始以速度向量(dx,dy)移动  直到走到曾经走过的位置  问  从哪里开始  可以经过最多的点数  dx和dy均与n互素 思路: 很容易想到以dx,dy移动的话  走的一定是一个圈  不会是出现"蝌蚪形"然后循环  注意题意最后一句  而且每个x坐标一定只经过一次 那么对于m个点  我们可以求出它是(0,y)这个点可以到达的  求解方法就是用拓展欧几里德解下面方程 (x+dx*k)%mod=xi , (y+dy*k)%mod=y

Codeforces 492E Vanya and Field 规律题

题目链接:点击打开链接 给定n*n的矩阵(0,0)->(n-1, n-1) m个苹果(下面m行给出苹果坐标)(dx, dy) 向量. 任选一个起点,用这个向量在矩阵里跑,问最多能采摘多少个苹果(坐标是%n, 即超过矩阵时 (x%n, y%n)) 输出起点. 思路: 把向量所在的点集写出来会发现一个起点一定经过了n个点,即至多只有n种起点 所以把点分成n个组即可. #pragma comment(linker, "/STACK:1024000000,1024000000") #in

bzoj4517: [Sdoi2016]排列计数--数学+拓展欧几里得

这道题是数学题,由题目可知,m个稳定数的取法是Cnm 然后剩下n-m本书,由于编号为i的书不能放在i位置,因此其方法数应由错排公式决定,即D(n-m) 错排公式:D[i]=(i-1)*(D[i-1]+D[i-2]); 所以根据乘法原理,答案就是Cnm * D(n-m) 接下来就是怎么求组合数的问题了 由于n≤1000000,因此只能用O(n)的算法求组合,这里用乘法逆元(inv[])来辅助求组合数 即 Cnm = n! / ((n-m)! * m!) = fac[n]*inv[n-m]*inv[

数论之拓展欧几里得求解不定方程和同余方程组(一)

今天接到scy的压缩包,开始做数论专题.那今天就总结一下拓展欧几里得求解不定方程和同余方程组. 首先我们复习一下欧几里得算法: 1 int gcd(int a,int b){ 2 if(b==0) return a; 3 return gcd(b,a%b);4 } 拓展欧几里得算法: 推导过程: 给出A和B,求它们的最大公约数,并且求出x和y,满足Ax+By=gcd(A,B). 当A=0时,x=0,y=1; 当A>0时, 因为exgcd(A,B,x,y)表示Ax+By=gcd(A,B) 而且ex

uva 10413 - Crazy Savages(拓展欧几里得)

题目链接:uva 10413 - Crazy Savages 题目大意:一座山有m个山洞,形成一个圈,现在有n个部落的人,每个部落一开始住在ci山洞,第2天会向后面移动pi个位置,一共会在这座山住li天.现在如果两个部落在同一个山洞相遇,则会发生战争,问说m最小时多少的时候,保证不会发生争斗. 解题思路:因为每个部落都有自己的存在时间,所以枚举m,然后枚举两个部落,判断他们有没有可能相遇. 假设在第k天: ci+k?pi=a?m???(1) cj+k?pj=b?m???(2) 由(2)-(1)得

uva 10548 - Find the Right Changes(拓展欧几里得)

题目链接:uva 10548 - Find the Right Changes 题目大意:给定A,B,C,求x,y,使得xA+yB=C,求有多少种解. 解题思路:拓展欧几里得,保证x,y均大于等于0,确定通解中t的取值. #include <cstdio> #include <cstring> #include <cmath> #include <algorithm> using namespace std; typedef long long ll; co

POJ 1061 青蛙的约会(拓展欧几里得)

青蛙的约会 Time Limit: 1000MS   Memory Limit: 10000KB   64bit IO Format: %I64d & %I64u Submit Status Description 两只青蛙在网上相识了,它们聊得很开心,于是觉得很有必要见一面.它们很高兴地发现它们住在同一条纬度线上,于是它们约定各自朝西跳,直到碰面为止.可是它们出发之前忘记了一件很重要的事情,既没有问清楚对方的特征,也没有约定见面的具体位置.不过青蛙们都是很乐观的,它们觉得只要一直朝着某个方向跳

poj2891 拓展欧几里得

1 //Accepted 164 KB 16 ms 2 //拓展欧几里得 3 //m=a1*x+b1 --(1) 4 //m=a2*(-y)+b2 --(2) 5 //->a1*x+a2*y=b2-b1 6 //由欧几里得算法可得上式的解 7 //由a*x+b*y=gcd(a,b) 8 //可得a(x+b)+b(y-a)=gcd(a,b) 9 //所以最小正整数解x=(x%b+b)%b; 10 //现考虑由(1)(2)两式得到的解m 11 //有x=m mod (a1*a2/gcd(a1,a2)