結果

問題 No.62 リベリオン(Extra)
ユーザー furafura
提出日時 2020-06-14 15:15:18
言語 C++17
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 86 ms / 5,000 ms
コード長 2,022 bytes
コンパイル時間 2,045 ms
コンパイル使用メモリ 204,824 KB
実行使用メモリ 4,380 KB
最終ジャッジ日時 2023-09-09 18:00:45
合計ジャッジ時間 2,850 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,380 KB
testcase_01 AC 1 ms
4,380 KB
testcase_02 AC 43 ms
4,376 KB
testcase_03 AC 43 ms
4,380 KB
testcase_04 AC 86 ms
4,380 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

#define rep(i,n) for(int i=0;i<(n);i++)

using namespace std;
using lint=long long;

template<class T>
T xgcd(T a,T b,T& x,T& y){
	if(b==0){ x=1; y=0; return a; }
	T g=xgcd(b,a%b,y,x);
	y-=a/b*x;
	return g;
}

pair<lint,lint> Chinese_remainder(vector<lint> a,vector<lint> b,vector<lint> m){
	int n=a.size();
	assert(n==b.size() && n==m.size());

	rep(i,n){
		a[i]%=m[i]; if(a[i]<0) a[i]+=m[i];
		b[i]%=m[i]; if(b[i]<0) b[i]+=m[i];

		lint x,y,g;
		g=xgcd(a[i],m[i],x,y);
		if(b[i]%g!=0) return {-1,-1};
		m[i]/=g;
		b[i]=(b[i]/g)*x%m[i];
	}

	lint x=0,M=1;
	rep(i,n){
		lint p,q,g;
		g=xgcd(M,m[i],p,q);
		if((b[i]-x)%g!=0) return {-1,-1};
		x+=M*(((b[i]-x)/g*p)%(m[i]/g));
		M*=m[i]/g;
	}
	x%=M; if(x<0) x+=M;
	return {x,M};
}

/*
	w * h の二次元トーラス上で, (sx, sy) からスタートして
	(dx, dy) 先にジャンプすることを d 回繰り返す間に (gx, gy) にたどり着けるか?
*/
bool solve(lint w,lint h,lint d,lint sx,lint sy,lint gx,lint gy,lint dx,lint dy){
	if(sx!=0 || sy!=0){
		return solve(w,h,d,0,0,(gx-sx+w)%w,(gy-sy+h)%h,dx,dy);
	}

	if(dx==0 && dy==0) return gx==0 && gy==0;
	if(dx==0){
		if(gx!=0) return false;
		auto [r,M]=Chinese_remainder({dy},{gy},{h});
		if(r==-1) return false;
		return r<=d;
	}
	if(dy==0){
		if(gy!=0) return false;
		auto [r,M]=Chinese_remainder({dx},{gx},{w});
		if(r==-1) return false;
		return r<=d;
	}

	vector<lint> a={dx,dy},b={gx,gy},m={w,h};
	auto [r,M]=Chinese_remainder(a,b,m);
	if(r==-1) return false;
	return r<=d;
}

int main(){
	int q; cin>>q;
	rep(_,q){
		lint w,h,d,sx,sy,gx,gy,vx,vy;
		cin>>w>>h>>d>>gx>>gy>>sx>>sy>>vx>>vy;

		lint g=gcd(abs(vx),abs(vy));
		d*=g;
		lint dx=(vx/g)%(2*w); if(dx<0) dx+=2*w;
		lint dy=(vy/g)%(2*h); if(dy<0) dy+=2*h;

		if(solve(2*w,2*h,d,sx,sy,gx,gy,dx,dy)
		|| solve(2*w,2*h,d,sx,sy,(2*w-gx)%(2*w),gy,dx,dy)
		|| solve(2*w,2*h,d,sx,sy,gx,(2*h-gy)%(2*h),dx,dy)
		|| solve(2*w,2*h,d,sx,sy,(2*w-gx)%(2*w),(2*h-gy)%(2*h),dx,dy)) puts("Hit");
		else puts("Miss");
	}

	return 0;
}
0