bzoj1406

这道题很有意思

我们解过线性同余方程，也解过同余方程

这个则是求x^2≡1 (mod p)

可以将问题转化为(x-1)(x+1)≡0 (mod p)

然后我们穷举一下p的约数i，

看i|x-1,p/i|x+1 或者i|x+1, p/i|x-1是否可行解

然后排序去重即可

var ans:array[0..1000010] of longint;
    a,n,i,b,t:longint;
    j:int64;

procedure sort(l,r:longint);
  var i,j,x,y:longint;
  begin
    i:=l;
    j:=r;
    x:=ans[(l+r) shr 1];
    repeat
      while ans[i]<x do inc(i);
      while x<ans[j] do dec(j);
      if not(i>j) then
      begin
        y:=ans[i];
        ans[i]:=ans[j];
        ans[j]:=y;
        inc(i);
        j:=j-1;
      end;
    until i>j;
    if l<j then sort(l,j);
    if i<r then sort(i,r);
  end;

begin
  readln(n);
  for i:=1 to trunc(sqrt(n)) do
    if n mod i=0 then
    begin
      b:=n div i;
      j:=1;
      while j<=n do
      begin
        if (j+1) mod i=0 then
        begin
          inc(t);
          ans[t]:=j;
        end;
        j:=j+b;  //保证x-1是b的倍数
      end;
      j:=b-1;
      while j<=n do
      begin
        if (j-1) mod i=0 then
        begin
          inc(t);
          ans[t]:=j;
        end;
        j:=j+b;
      end;
    end;
  sort(1,t);
  for i:=1 to t do
    if ans[i]<>ans[i-1] then writeln(ans[i]);
end.