仪仗队(points)

仪仗队(points)

【题意】

求存在多少个i,j互质的二元组(i,j)(0<=i,j<=n)

【输入】

一个数字表示n(n<=40000)

【输出】

一个数字,表示有多少个满足条件的二元组


枚举2到4000,对于每一个求其欧拉函数,求和乘二+3


program points;
var
  n,i,j,k,tot,ans:longint;
  prime:array [0..10001] of longint;
  yes:array [0..40001] of boolean;

begin
  assign(input,'points.in');
  reset(input);
  assign(output,'points.out');
  rewrite(output);
  read(n);
  dec(n);
  tot:=0;
  fillchar(yes,sizeof(yes),false);
  for i:=2 to n do
    if not yes[i] then
      begin
        inc(tot);
        prime[tot]:=i;
        yes[i]:=true;
        for j:=i to n div i do
          yes[i*j]:=true;
      end;
  ans:=0;
  for i:=2 to n do
    begin
      k:=i;
      j:=1;
      while (j<=tot)and(prime[j]<=i) do
        begin
          if i mod prime[j] = 0 then k:=k*(prime[j]-1) div (prime[j]);
          inc(j);
        end;
      ans:=ans+k;
    end;
  ans:=ans*2+3;
  writeln(ans);
  close(input);
  close(output);
end.





评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值