数论之题,多半乃只可意会不可言传者,然吾必试之,盖练吾文——喵喵喵
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=∑i=1n(k−⌊ki⌋∗i)
=nk−∑i=1n(⌊ki⌋∗i)
差不多画到这里就ok了,然后我们还可以知道求和上方的n是要和k取min的(因为超过k的除出来肯定是0没有意义)
但是O(n)的效率我们接受不了,可以发现 ⌊ki⌋ 在一段连续的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
欲求
∑i=1ngcd(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);
}