题目描述
对N!进行质因子分解。
输入输出格式
输入格式:
输入数据仅有一行包含一个正整数N,N<=10000。
输出格式:
输出数据包含若干行,每行两个正整数p,a,中间用一个空格隔开。表示N!包含a个质因子p,要求按p的值从小到大输出。
输入输出样例
输入样例#1:
10
输出样例#1:
2 8
3 4
5 2
7 1
说明
10!=3628800=(2^8)(3^4)(5^2)*7
分析:不用高精度,先扫描一遍,打标记,一层一层除,最后把不是0的输出
var
a:array[0..100000] of longint;
i,j,x,max,n:longint;
function prime(k:longint):boolean;
var i:longint;
begin
for i:=2 to trunc(sqrt(k)) do if k mod i=0 then exit(false);
exit(true);
end;
procedure work(x,y:longint);
begin
while(x mod y=0) do begin
x:=x div y;
inc(a[y]);
end;
end;
begin
readln(n);
for i:=1 to n do begin
x:=i;
if prime(x) then inc(a[x]);
max:=x;
for j:=2 to i div 2 do begin
if (x mod j=0)and(prime(j)) then begin
work(x,j);
if j>max then max:=j;
end;
end;
end;
for i:=2 to max do if a[i]<>0 then writeln(i,' ',a[i]);
end.