• 一道不知道哪里来的数论题


    Pro:
    在一张(n*m)的台球桌上上打台球,初始位置(x,y),初速度(v_x,v_y),打到边缘会无能量损失的反弹。
    问最少反弹多少次可以进洞。
    Sol:
    考虑把台球桌无限延伸。
    然后发现可以用剩余系来搞。
    大概就是

    [egin{align*} x+v_x*t&=0 (mod n) \ y+v_y*t&=0 (mod m) end{align*} ]

    即:

    [egin{align*} v_x*t&=k_1*n-x \ v_y*t&=k_2*m-y end{align*} ]

    联立消一下(t)

    [egin{align*} (n*v_y)*k_1-(m*v_x)*k_2=x*v_y-y*v_x end{align*} ]

    然后需要最小化(k_1+k_2-2)
    (exgcd)冲一下就行了

    #include<bits/stdc++.h>
    #define N 110000
    #define db double
    #define ll long long
    #define ldb long double
    #define ull unsigned long long
    using namespace std;
    const ll h=3,ki=149,mo=998244353;
    const double pi=acos(-1),eps=1e-7,inf=1e18+7;
    ll mod(ll x){return (x%mo+mo)%mo;}
    ll inc(ll x,ll k){x+=k;return x<mo?x:x-mo;}
    ll dec(ll x,ll k){x-=k;return x>=0?x:x+mo;}
    ll ksm(ll x,ll k)
    {
    	ll ans=1;
    	while(k){if(k&1)ans=1ll*ans*x%mo;k>>=1;x=1ll*x*x%mo;}
    	return mod(ans);
    }
    ll inv(ll x){return ksm(x,mo-2);}
    ll read()
    {
    	char ch=0;ll x=0,flag=1;
    	while(!isdigit(ch)){ch=getchar();if(ch=='-')flag=-1;}
    	while(isdigit(ch)){x=(x<<3)+(x<<1)+ch-'0',ch=getchar();}
    	return x*flag;
    }
    void write(ll x)
    {
    	if(!x)return (void)putchar(48);
    	if(x<0)putchar(45),x=-x;
    	ll len=0,p[20];
    	while(x)p[++len]=x%10,x/=10;
    	for(ll i=len;i>=1;i--)putchar(p[i]+48);
    }
    void writeln(ll x){write(x),putchar('
    ');}
    ll X,Y;
    ll exgcd(ll a,ll b)
    {
    	if(!b){X=1;Y=0;return a;}
    	ll d=exgcd(b,a%b);
    	ll t=X;X=Y;Y=t-(a/b)*Y;
    	return d;
    }
    void work()
    {
    	ll n=read(),m=read(),x=read(),y=read(),vx=read(),vy=read(); 
    	if(!vx||!vy){printf("NO
    ");return;}
    	ll a=n*vy,b=m*vx,c=x*vy-y*vx;
    	ll d=exgcd(a,b),p1=b/d,p2=a/d,k=c/d;
    	
    	if(c%d)printf("NO
    ");
    	else
    	{
    		printf("YES
    ");
    		X*=k;Y*=k;X=(X%p1+p1)%p1;if(!X)X+=p1;Y=(c-a*X)/b;
    		while(Y>0)X+=p1,Y-=p2;
    		writeln(X-Y-2); 
    	}
    }
    int main()
    {
    	ll t=read();
    	for(ll i=1;i<=t;i++)work();
    	return 0;
    } 
    
  • 相关阅读:
    bzoj2064分裂(dp)
    关于逆元的学习笔记(尚未完成)
    线性基学习笔记
    bzoj2460元素(线性基,贪心)
    CF280C Game on tree(期望dp)
    XJOI contest800
    关闭极域电子教室学生端的各种方法
    大坑!有网,电脑qq登不上去!!
    4.25 ZJOI2017 10AM
    bzoj1621
  • 原文地址:https://www.cnblogs.com/Creed-qwq/p/14818899.html
Copyright © 2020-2023  润新知