洛谷P3601签到题（欧拉函数）

题目背景

这是一道签到题！

建议做题之前仔细阅读数据范围！

 

题目描述

我们定义一个函数：qiandao(x)为小于等于x的数中与x不互质的数的个数。

这题作为签到题，给出l和r，要求求。

 

输入输出格式

输入格式：

一行两个整数，l、r。

 

输出格式：

一行一个整数表示答案。

 

输入输出样例

 

输入样例#1：

233 2333

输出样例#1：

1056499

输入样例#2：

2333333333 2333666666

输出样例#2：

153096296

说明

对于30%的数据，。

对于60%的数据，。

对于100%的数据，，。

 

 

qiandao(x)=x-phi(x)，转化为求欧拉函数，而x只可能有一个大于n^0.5的素因子，所以只要求到n^0.5的素数就行了。还是批量处理欧拉函数的值，否则TLE啊。

program rrr(input,output);
const
  cs=666623333;
var
  a:array[0..1000010]of boolean;
  f,c:array[0..1000010]of int64;
  i,j,n:longint;
  l,r,k,ans:int64;
begin
   assign(input,'r.in');assign(output,'r.out');reset(input);rewrite(output);
   readln(l,r);
   n:=trunc(sqrt(r));
   fillchar(a,sizeof(a),true);
   k:=l;while k<=r do begin f[k-l]:=k;c[k-l]:=k;inc(k); end;
   for i:=2 to n do
      if a[i] then
         begin
            k:=l div i*i;if k<l then k:=k+i;
            while k<=r do
               begin
                  f[k-l]:=f[k-l] div i*(i-1);
                  while c[k-l] mod i=0 do c[k-l]:=c[k-l] div i;
                  k:=k+i;
               end;
            j:=i<<1;while j<=n do begin a[j]:=false;j:=j+i; end;
         end;
   k:=l;while k<=r do begin if c[k-l]>1 then f[k-l]:=f[k-l] div c[k-l]*(c[k-l]-1);inc(k); end;
   ans:=0;
   k:=l;while k<=r do begin ans:=(ans+k-f[k-l]) mod cs;inc(k); end;
   write(ans);
   close(input);close(output);
end.