Cod sursa(job #118015)

Utilizator bogdan88Bogdan Popescu bogdan88 Data 22 decembrie 2007 22:29:36
Problema Sum Scor 45
Compilator fpc Status done
Runda Arhiva de probleme Marime 2.06 kb
var fi,fo:text;
    n,ct,m,l,phi:int64;
    prime:array[1..50000]of longint;
    puteri:array[1..100000]of longint;
    prime2:array[1..100000]of byte;
    primee:array[1..90000]of longint;
    ap:array[1..100000]of int64;
procedure gen(ba:longint);
var i,j:int64;
begin
  i:=3;
  primee[1]:=2;
  l:=1;
  prime2[1]:=1;
  while i<=trunc(sqrt(ba)) do
    begin
      if prime2[i]=0 then
        begin
          j:=i*i;
          inc(l);
          primee[l]:=i;
          while j<=ba do
             begin
                prime2[j]:=1; inc(j,i);
             end;
        end;
      inc(i,2);
    end;
end;
procedure desc(n:longint);
var i:longint;
    r,nenner,zahler:int64;
begin
  {ct:=0;} nenner:=1; zahler:=1;
  for i:=1 to l do
    begin
      if n mod primee[i]=0 then
        begin
          nenner:=nenner*primee[i];
          zahler:=(primee[i]-1)*zahler;
          {inc(ct);
          prime[ct]:=primee[i];}
          {r:=1;
          while n mod primee[i]=0 do}
            {begin
              n:=n div primee[i];
              r:=r*primee[i];
            end;
          puteri[primee[i]]:=r;}
        end;
        if i>=n then break;
    {  if n=1 then exit;
      if prime2[n]=0 then begin inc(ct); prime[ct]:=n; puteri[n]:=n; exit; end;}
    end;
  {if n>1 then
    begin
      {inc(ct);
      prime[ct]:=n;
      puteri[n]:=n;}
      nenner:=nenner*n;
      zahler:=zahler*(n-1); end;}
    phi:=n*zahler div nenner;
end;
procedure solv(n:longint);
var rez:int64;
    i:longint;
begin
  desc(n);
  {phi:=1;
  for i:=1 to ct do
     phi:=phi*(prime[i]-1)*puteri[prime[i]] div prime[i];}
  rez:=n*phi shl 1;
  writeln(fo,rez);
  ap[n]:=rez;
end;
var i:longint;
begin
  assign(fi,'sum.in'); reset(fi);
  assign(fo,'sum.out'); rewrite(fo);
  read(fi,m);
  gen(100000);
  for i:=1 to m do
    begin
      read(fi,n);
      if ap[n]<>0 then writeln(fo,ap[n])
        else
          if (n and 1=1)and(prime2[n]=0)then writeln(fo,(n-1)*n shl 1)
            else solv(n);
    end;
  close(fi);
  close(fo);
end.