【bzoj2154】Crash的数字表格

  • 本文为博主原创,未经许可不得转载

题目大意

一句话题意:给定n,m,求$\sum_{i=1}^{n}\sum_{j=1}^{m}lcm(i,j)(n,m<=10^7)$

-
-


题解

$ans=\sum_{i=1}^{n}\sum_{j=1}^{m}lcm(i,j)=\sum_{i=1}^{n}\sum_{j=1}^{m}\frac{i\times j}{gcd(i,j)}$

令$d=gcd(i,j)$,$F(x,y)=\sum_{1<=i<=x,1<=j<=y,gcd(i,j)=1}^{ }i \times j$

然后我们枚举$d$,则有$ans=\sum_{d=1}^{min(n,m)}\frac{d^2 \times F(\left [ \frac{n}{d} \right ],\left [ \frac{m}{d} \right ])}{d}=\sum_{d=1}^{min(n,m)}d\times F(\left [ \frac{n}{d} \right ],\left [ \frac{m}{d} \right ])$

令$Sum(x,y)=\sum_{i=1}^{x}\sum_{j=1}^{y}i \times j=\frac{x(x-1)}{2}\times\frac{y(y-1)}{2}$

显然$Sum(x,y)=\sum_{i=1}^{min(x,y)}i^2 \times F(\left [ \frac{x}{i} \right ],\left [ \frac{y}{i} \right ])$

通过莫比乌斯反演可得$F(x,y)=\sum_{i=1}^{min(x,y)}i^2\times \mu(i)\times Sum(\left [ \frac{x}{i} \right ],\left [ \frac{y}{i} \right ])$

这两个式子都可以进行$O(\sqrt n)$的计算,所以总的时间复杂度是$O(n)$

参考代码

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<ctime>
#include<cmath>
#include<algorithm>
using namespace std;
typedef long long ll;
#define FILE "read"
#define cmax(a,b) a=max(a,b)
#define cmin(a,b) a=min(a,b)
const int MAXN=(int)1e7,D(50),mod=20101009;
int n,m,cnt,prime[MAXN+D],check[MAXN+D],mu[MAXN+D];
ll ans,sum[MAXN+D];
namespace INIT{
char buf[1<<15],*fs,*ft;
inline char getc(){return (fs==ft&&(ft=(fs=buf)+fread(buf,1,1<<15,stdin),fs==ft))?0:*fs++;}
inline int read(){
int x=0,f=1; char ch=getc();
while(!isdigit(ch)) {if(ch=='-') f=-1; ch=getc();}
while(isdigit(ch)) {x=x*10+ch-'0'; ch=getc();}
return x*f;
}
}using namespace INIT;
void pre(){
mu[1]=1;
for(ll i=2;i<=min(n,m);++i){
if(!check[i]) prime[++cnt]=i,mu[i]=-1;
for(ll j=1;j<=cnt&&prime[j]*i<=min(n,m);++j){
check[prime[j]*i]=1;
if(i%prime[j]==0) {mu[i*prime[j]]=0; break;}
mu[i*prime[j]]=-mu[i];
}
}
for(ll i=1;i<=min(n,m);++i) sum[i]=((sum[i-1]+(mu[i]*i*i%mod)+mod)%mod+mod)%mod;
}
ll cal(ll x,ll y){return ((x*(x+1)/2%mod)*(y*(y+1)/2%mod))%mod;}
ll get(ll x,ll y){
if(x>y) swap(x,y); ll temp(0);
for(ll i=1,last;i<=x;i=last+1){
last=min(x/(x/i),y/(y/i));
temp+=((sum[last]-sum[i-1])*cal(x/i,y/i))%mod;
temp%=mod;
}
return temp;
}
int main(){
freopen(FILE".in","r",stdin);
freopen(FILE".out","w",stdout);
n=read(); m=read(); pre();
if(n>m) swap(n,m);
for(ll i=1,last;i<=n;i=last+1){
last=min(n/(n/i),m/(m/i));
ans+=(i+last)*(last-i+1)/2%mod*get(n/i,m/i)%mod;
ans%=mod;
}
printf("%lld\n",(ans+mod)%mod);
return 0;
}
文章目录
  1. 1. 题目大意
  2. 2. 题解
  3. 3. 参考代码
,