結果

問題 No.62 リベリオン(Extra)
ユーザー fura
提出日時 2020-06-14 15:15:18
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 116 ms / 5,000 ms
コード長 2,022 bytes
コンパイル時間 3,757 ms
コンパイル使用メモリ 198,712 KB
最終ジャッジ日時 2025-01-11 04:15:37
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 3
権限があれば一括ダウンロードができます

ソースコード

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;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0