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[..] 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 ];
     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-;
       end;
     until i>j;
     if l<j then sort(l,j);
     if i<r then sort(i,r);
   end;

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