线性筛

为了写更优美的筛法,也为纪念chp和xjh。

var
 prime,notprime,phi:array [1..100000] of longint;
 n,i,j,cnt:longint;
begin
 readln(n);
 for i:=2 to n do
  begin
   if notprime[i]=0 then
    begin
     cnt:=cnt+1;
     prime[cnt]:=i;
     phi[i]:=i-1;
    end;
     for j:=1 to cnt do
      begin
       if prime[j]*i>n then break;
       notprime[i*prime[j]]:=1;
       if i mod prime[j]=0 then
        begin
         phi[i*prime[j]]:=phi[i]*prime[j];
         break;
        end;
       phi[i*prime[j]]:=phi[i]*(prime[j]-1);
      end;
  end;
 for i:=1 to cnt do
  writeln(prime[i]);
 for i:=1 to n do
  writeln(phi[i]);
end.

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值