PKU1061 青蛙的约会

其实就是求解一个线性同余方程

用扩展gcd做:

  1. function gcd(a,b:int64;var x,y:int64):int64;
  2.   var
  3.     t:int64;
  4.   begin
  5.     if b=0 then begin
  6.       gcd:=a;
  7.       x:=1;
  8.       y:=0;
  9.     end else begin
  10.       gcd:=gcd(b,a mod b,x,y);
  11.       t:=x;
  12.       x:=y;
  13.       y:=t-(a div b)*y;
  14.     end;
  15.   end;

然后求最小正整数解时直接mod一下就可以了~

  1.     ans:=gcd(n,l,xx,yy);
  2.     if x mod ans<>0 then begin
  3.       writeln('Impossible');
  4.       exit;
  5.     end;
  6.     xx:=xx*x div ans;
  7.     if xx>0 then
  8.       xx:=xx mod (l div ans);
  9.     if xx<0 then begin
  10.       xx:=xx mod (l div ans);
  11.       if xx<>0 then xx:=xx+l div ans;
  12.     end;
  13.     writeln(xx);
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值