文章标题

题目描述

对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.