欧拉函数O(sqrt(n))

利用的是公式:φ(x)=x(1-1/p1)(1-1/p2)(1-1/p3)(1-1/p4)…..(1-1/pn),其中p1, p2……pn为x的所有质因数,x是不为0的整数。φ(1)=1(唯一和1互质的数(小于等于1)就是1本身)。 (注意:每种质因数只一个。比如12=2*2*3那么φ(12)=12*(1-1/2)*(1-1/3)=4)。

program pro;
var
        oula:array[0..10000]of longint;
        n,i,j:longint;

procedure makeoula(x:longint);
var
        ii,jj:longint;
begin
        for ii:=1 to x do oula[ii]:=ii;

        for ii:=2 to x do
        begin
                if oula[ii]=ii then
                begin
                        jj:=ii;
                        while jj<=x do
                        begin
                                oula[jj]:=oula[jj]div ii*(ii-1);
                                inc(jj,ii);
                        end;
                end;
        end;
        for ii:=1 to x do writeln(ii,' ',oula[ii]);
end;



begin
assign(input,'oula.in'); reset(input);
assign(output,'oula.out'); rewrite(output);
        readln(n);
        makeoula(n);
close(input);
close(output);
end.



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

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值