sumaD=0; liczbaD = 0;
dzielniki( n ) = { local( j, d=2, g=n );
sumaD = 1; liczbaD = 1;
while( d < g,
j = n\d;
if( d * j == n,
sumaD += d + j;
liczbaD += 2 );
d++ ; g = j
);
if( g*j == n, \\ n == x^2
liczbaD--; sumaD -= j )
};