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);
發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章