3520. 【NOIP2013模拟11.7B组】原根(math)

Description

Input

有且只有一个正整数m。

Output

以递增序依次输出模m的所有原根，每行输出一个原根。
如果不存在模m的原根，输出-1。

Sample Input

7

Sample Output

3
5

Data Constraint

50%的数据，m≤ 200。
100%的数据，m ≤ 10000。

思路

先算出m以内与m互质（gcd(m,i)=1）的个数x，即为m的欧拉函数，再在这些互质的数中找出满足a^1到a^x-1每个数mod m都不为1，且a^x mod m为1，那么为m的原根。

function gcd(x,y:longint):longint;
var
  z:longint;
begin
  z:=x mod y;
  while z<>0 do
    begin
      x:=y;
      y:=z;
      z:=x mod y;
    end;
  gcd:=y;
end;
var
  l,i,j,x,m:longint;
  c,f:boolean;
  a:array[1..20000] of longint;
begin
  assign(input,'math.in');reset(input);
  assign(output,'math.out');rewrite(output);
  readln(m);
  l:=0;
  if m=1 then begin writeln(1);exit;end;
  f:=false;
  for i:=1 to m do
    if gcd(m,i)=1 then begin inc(l);a[l]:=i;end;
  for i:=1 to l do
    begin
      x:=1;
      c:=true;
      for j:=1 to l-1 do
        begin
          x:=(x*a[i]) mod m;
          if x=1 then begin c:=false;break;end;
        end;
      if (x*a[i] mod m=1)and(c) then begin writeln(a[i]);f:=true;end;
    end;
  if not f then writeln(-1);
  close(input);close(output);
end.