PKU1061 青蛙的约会

来源:互联网 发布:淘宝好评返现被投诉 编辑:程序博客网 时间:2024/04/29 15:05

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

用扩展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);
原创粉丝点击