Cod sursa(job #743154)

Utilizator vasile_pojogaPojoga Vasile vasile_pojoga Data 3 mai 2012 15:54:25
Problema Suma si numarul divizorilor Scor 10
Compilator fpc Status done
Runda Arhiva educationala Marime 2.21 kb
program suma_dv_nr;
var t:integer;
    n,i,s,k:longint;
    p:array[1..1000005] of boolean;
    x:array[1..1000005] of longint;
procedure ciur;
const nmax=1000005;
var i,j:longint;
begin
x[1]:=1; k:=1;
for i:=2 to nmax do
         if not p[i] then
                     begin
                     j:=i; k:=k+1; x[k]:=i;
                     while j<nmax do begin
                                     p[j]:=true;
                                     j:=j+i;
                                     end;
                     p[i]:=false;
                     end;
p[2]:=true;
end;

procedure solve(n:longint);
const md=9973;
var i,s,nr,d,s1:longint;

function putere(a,b:longint):longint;
var i:longint;
begin
putere:=a;
for i:=2 to b do putere:=putere*a;
end;

begin
nr:=1;s:=1;s1:=1;i:=2;d:=0;
while n>0 do begin
             if n mod x[i]=0 then begin
                                  d:=d+1;
                                  n:=n div x[i];
                                  end
                             else begin
                                  if d>0 then begin
                                              nr:=nr*(d+1);
                                              s:=s*(putere(x[i],d+1)-1);
                                              s1:=s1*(x[i]-1);
                                              d:=0;
                                              end;
                                  if p[n]=false then begin
                                                     if n=1 then break;
                                                     nr:=nr*2;
                                                     s:=s*(putere(n,2)-1);
                                                     s1:=s1*(n-1);
                                                     n:=0;
                                                     end;
                                  i:=i+1;
                                  end;
             end;
s:=(s div s1) mod md;
writeln(nr,' ',s);
end;

begin
ciur;
assign(input,'ssnd.in'); reset(input);
assign(output,'ssnd.out'); rewrite(output);
readln(t);
for k:=1 to t do begin
                 readln(n);
                 solve(n);
                 end;
close(output);
end.