【扩展欧几里德】Vijos P1009 清帝之惑之康熙

44 篇文章 0 订阅
28 篇文章 0 订阅

题目链接:

  https://vijos.org/p/1009

题目大意

  两个人,一个在坐标x,每天走m,一个在坐标y,每天走n,坐标长L,问几天后碰面。

题目思路:

  【扩展欧几里德】

  根据同余方程的ax+by=c,可以得出a=n-m,b=l,c=x-y 之后扩展欧几里得求线性同余方程。

  题目数据较大需要使用long long 


//
//by coolxxx
//
#include<iostream>
#include<algorithm>
#include<string>
#include<iomanip>
#include<memory.h>
#include<time.h>
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<stdbool.h>
#include<math.h>
#define min(a,b) ((a)<(b)?(a):(b))
#define max(a,b) ((a)>(b)?(a):(b))
#define abs(a) ((a)>0?(a):(-(a)))
#define lowbit(a) (a&(-a))
#define sqr(a) (a)*(a)
#define swap(a,b) (a)^=(b),(b)^=(a),(a)^=(b)
#define eps 1e-8
#define MAX 0x7f7f7f7f
#define INF 20000
#define PI 3.1415926535897
#define N 100004
using namespace std;
int n,m,lll,ans,cas;
long long x,y,a,b,c,d;
long long exgcd(long long a,long long b)
{
	if(!b)
	{
		x=1;y=0;
		return a;
	}
	d=exgcd(b,a%b);
	cas=x;
	x=y;
	y=cas-(a/b)*y;
	return d;
}
int main()
{
	#ifndef ONLINE_JUDGE
//	freopen("1.txt","r",stdin);
//	freopen("2.txt","w",stdout);
	#endif
	int i,j,k;
	while(~scanf("%d",&n) && n)
	{
		x=n;
		scanf("%lld%lld%lld%lld",&y,&m,&n,&b);//(n-m)X mod b = x-y; aX mod b = c
		a=((m-n)%b+b)%b;
		c=((y-x)%b+b)%b;
		d=exgcd(a,b);
		if(c%d!=0)
		{
			puts("Impossible");
			continue;
		}
		a/=d;b/=d;c/=d;
		exgcd(a,b);
		x=(x%b+b)%b;
		x=(x*c)%b;
		printf("%lld\n",x);
	}
	return 0;
}


/*
//
	
//
*/


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值