BZOJ 4173: 数学

4173: 数学

Time Limit: 10 Sec  Memory Limit: 256 MB
Submit: 462  Solved: 227
[Submit][Status][Discuss]

Description

Input

输入文件的第一行输入两个正整数 。

Output

如题

Sample Input

5 6

Sample Output

240

HINT

N,M<=10^15

Source

 

[Submit][Status][Discuss]

据说单个欧拉函数$phi(N)$可以在$O(\sqrt{N})$的时间内求出来，所以只关心$\sum_{k\in S(n,m)}{\phi{k}}$是什么。看下PoPoQQQ。

n%k+m%k≥k等价于⌊n+mk⌋−⌊nk⌋−⌊mk⌋=1 
无视掉前面的φ(n)∗φ(m)的话答案就是 
∑n%k+m%k≥kφ(k) 
=∑n+mk=1φ(k)∗⌊n+mk⌋−∑nk=1φ(k)∗⌊nk⌋−∑mk=1φ(k)∗⌊mk⌋

那么∑nk=1φ(k)∗⌊nk⌋又是什么呢？ 
∑ni=1i=∑ni=1∑k|iφ(k)=∑nk=1φ(k)∗⌊nk⌋

因此答案为φ(n)∗φ(m)∗(∑n+mi=1i−∑ni=1i−∑mi=1i)=φ(

 #include <cstdio>

 typedef long long lnt;

 __inline lnt phi(lnt n)
 {
     lnt r = n;

     for (lnt i = ; i*i <= n; ++i)
         if (n % i == )
         {
             r = (r / i) * (i - );
             while (n % i == )n /= i;
         }

     if (n != )r = (r / n) * (n - );

     return r;
 }

 const lnt mod = 998244353LL;

 lnt n, m;

 signed main(void)
 {
     scanf("%lld%lld", &n, &m);
     printf("%lld\n",
         ((phi(n) % mod) * (phi(m) % mod) % mod)
     *    ((n % mod) * (m % mod) % mod) % mod);
 }

@Author: YouSiki