数论小练

数论之题，多半乃只可意会不可言传者，然吾必试之，盖练吾文——喵喵喵

bzoj1968

n/k也是求1~n中含k为约数 的个数的意思（因为每k个数字就会有一个数含k为约数）
比如10/3=3（3,6,9）都是含3为约数
这道题让求前缀和，就从1-n枚举，求n/i的累加和

代码：

#include <cstdio>
using namespace std;
int main()
{
    int ans=0,n,i;
    scanf("%d",&n);
    for (i=1;i<=n;i++) 
      ans+=n/i;
    printf("%d",ans);   
}

bzoj1257

先画一波柿子

$$ans=\sum_{i=1}^{n}(k-\lfloor \frac ki \rfloor*i)$$
$$=nk-\sum_{i=1}^{n}(\lfloor \frac ki \rfloor*i)$$
差不多画到这里就ok了，然后我们还可以知道求和上方的n是要和k取min的（因为超过k的除出来肯定是0没有意义）
但是O(n)的效率我们接受不了，可以发现 $\lfloor \frac ki \rfloor$ 在一段连续的i中是一样的，根据这个，我们可以采用等差数列求和的方法实现阶段求和。k/(k/i)是一个阶段的终止点

代码：

#include <cstdio>
#include <iostream>
#define LL long long
using namespace std;
int n,k;
int main()
{
    scanf("%d%d",&n,&k);
    LL ans=(LL)n*k;
    n=min(n,k);
    int i,j=1;
    for (i=1;i<=n;i=j+1)
    {
        j=min(n,k/(k/i));
        ans-=(LL)(k/i)*(j-i+1)*(i+j)/2;
    }
    printf("%lld",ans);
}

bzoj2705

欲求

$$\sum_{i=1}^{n}gcd(i,n)$$
我们知道gcd可能出现的数字是n的因数，那我们枚举n的因数k，只要再乘上[以k为最大公约数的数字个数]s(k)，累加一下就是答案
$gcd(i,n)=k$ 变个形就是 $gcd(i/k,n/k)=1$ 那么i/k的个数就是s(k)，这个gcd的柿子看起来眼熟？这不就是个欧拉函数phi嘛，求一下就ok啦

代码：

#include <cstdio>
#include <cmath>
#include <cstring>
#include <iostream>
#include <algorithm>
#define LL long long
using namespace std;
int i;LL n,ans;
LL phi(int n)
{
    LL ans=1;
    for (int i=2;i*i<=n;i++)
      if (n%i==0)
    {
        ans*=i-1; n/=i;
        while (n%i==0) n/=i,ans*=i;
    }
    if (n>1) ans*=n-1;
    return ans;
}
int main()
{
    scanf("%d",&n);
    int m=sqrt(n);
    for (i=1;i<=m;i++)
      if (n%i==0)
      {
        ans+=(LL)phi(n/i)*i;
        ans+=(LL)phi(n/(n/i))*(n/i);
      }

    printf("%lld",ans);
}