欧拉函数

欧拉函数
n<=2311

const
 maxn=1006350;
var
 check:array[0..maxn]of boolean;
 prime:array[0..maxn]of longint;
 i,j:longint;
 n,s,ans,len,t:longint;
begin
 readln(n);
 s:=n; ans:=n; n:=trunc(sqrt(n)); len:=0;
 for i:=2 to n do
  begin
   if check[i]=false
   then begin inc(len); prime[len]:=i; end;
   for j:=1 to len do
    begin
     if prime[j]*i>n
     then break;
     check[prime[j]*i]:=true;
     if i mod prime[j]=0
     then break;
    end;
  end;
 for i:=1 to len do
  begin
   if s mod prime[i]=0
   then
    begin
     ans:=(ans div prime[i])*(prime[i]-1);
     while s mod prime[i]=0 do
      s:=s div prime[i];
    end;
  end;
 if s<>1
 then ans:=(ans div s)*(s-1);
 writeln(ans);
end.
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值