15 lines
248 B
Text
15 lines
248 B
Text
do(lim)={
|
|
my(prim,total,P);
|
|
lim\=1;
|
|
for(m=2,sqrtint(lim\2),
|
|
forstep(n=1+m%2,min(sqrtint(lim-m^2),m-1),2,
|
|
P=2*m*(m+n);
|
|
if(gcd(m,n)==1 && P<=lim,
|
|
prim++;
|
|
total+=lim\P
|
|
)
|
|
)
|
|
);
|
|
[prim,total]
|
|
};
|
|
do(100)
|