【Luogu】 P3312 [SDOI2014]数表

最新推荐文章于 2024-08-08 15:12:59 发布

Farmer_D

最新推荐文章于 2024-08-08 15:12:59 发布

阅读量44

点赞数

分类专栏： Luogu 文章标签：算法

本文链接：https://blog.csdn.net/djx2021/article/details/129098534

版权

Luogu 专栏收录该内容

45 篇文章 1 订阅

订阅专栏

题目描述

点击打开链接

题目解法

我们令 $g (x)$ 为 $x$ 的约数和
答案为：
$\sum_{i=1}^{n}\sum_{j=1}^{m}[g[gcd(i,j)]<=a]g[gcd(i,j)]$
按照套路枚举 $g c d (i, j)$ ：
$\sum_{d=1}^{min(n,m)}[g(d)<=a]\sum_{i'=1}^{\left \lfloor \frac{n}{d} \right \rfloor}\sum_{j'=1}^{\left \lfloor \frac{m}{d} \right \rfloor}[gcd(i',j')=1]g(d)$
如果 $x > 1$ ，那么 $\sum_{d|x} \mu(d)$ 为 $0$ ：
$\sum_{d=1}^{min(n,m)}[g(d)<=a]g[d]*\sum_{i'=1}^{\left \lfloor \frac{n}{d} \right \rfloor}\sum_{j'=1}^{\left \lfloor \frac{m}{d} \right \rfloor}\sum_{d'|i'且d'|j'}\mu(d')$
交换 $i^{'}, j^{'}, d^{'}$ 的循环顺序：
$\sum_{d=1}^{min(n,m)}[g(d)<=a]g[d]*\sum_{d'=1}^{min(\left \lfloor \frac{n}{d} \right \rfloor,\left \lfloor \frac{m}{d} \right \rfloor)}\mu(d')\left \lfloor \frac{n}{dd'} \right \rfloor\left \lfloor \frac{m}{dd'} \right \rfloor$
令 $T = d d^{'}$ ，且先枚举 $T$ 来交换循环顺序：
$\sum_{T=1}^{min(N,M)}\left \lfloor \frac{n}{T} \right \rfloor\left \lfloor \frac{m}{T} \right \rfloor\sum_{d|T}[g(d)<=a]g(d)*\mu(\frac{T}{d})$
我么令 $h (T)$ 为 $\sum_{d|T}[g(d)<=a]g(d)*\mu(\frac{T}{d})$
我们可以把询问离线下来，按照 $a$ 从小到大排序
然后用每个 $g (d)$ 来更新 $h (T)$ ，这样全部更新一遍的复杂度大概是 $O (n l o g n)$
每个答案可以用数论分块 $O(\sqrt n)$ 求出
这里的 $h (T)$ 因为需要满足单点修改和区间求和的操作，所以可以用常数较小的树状数组来实现
时间复杂度 $O(nlogn\sqrt n)$

#include <bits/stdc++.h>
#define lowbit(x) x&-x
#define int long long
using namespace std;
typedef pair<int,int> pii;
const int N(100100),Q(20100);
struct Query{
	int n,m,a,id;
}query[Q];
int g[N],h[N],ANS[Q];
int pr[N],cnt,v[N],mu[N];
pii sortg[N];
inline int read(){
	int FF=0,RR=1;
	char ch=getchar();
	for(;!isdigit(ch);ch=getchar()) if(ch=='-') RR=-1;
	for(;isdigit(ch);ch=getchar()) FF=(FF<<1)+(FF<<3)+ch-48;
	return FF*RR;
}
bool cmp(const Query &x,const Query &y){
	return x.a<y.a;
} 
void init(){
	mu[1]=1;
	for(int i=2;i<N;i++){
		if(!v[i]) v[i]=pr[++cnt]=i,mu[i]=-1;
		for(int j=1;i<=(N-1)/pr[j];j++){
			v[pr[j]*i]=pr[j];
			if(v[i]==pr[j]) break;
			mu[pr[j]*i]=-mu[i];
		}
	}
	for(int i=1;i<N;i++)
		for(int j=i;j<N;j+=i)
			g[j]+=i;
	for(int i=1;i<N;i++)
		sortg[i]=make_pair(g[i],i);
	sort(sortg+1,sortg+N);
}
void add(int x,int c){
	for(;x<N;x+=lowbit(x))
		h[x]+=c;
}
int ask(int x){
	int res=0;
	for(;x;x-=lowbit(x))
		res+=h[x];
	return res;
}
int get_next(int n,int k){
	return n/(n/k);
}
signed main(){
	init();
	int q=read();
	for(int i=1;i<=q;i++)
		query[i].n=read(),query[i].m=read(),query[i].a=read(),query[i].id=i;
	sort(query+1,query+q+1,cmp);
	int p=1;
	for(int i=1;i<=q;i++){
		while(p<N&&sortg[p].first<=query[i].a){
			int d=sortg[p].second;
			for(int j=d;j<N;j+=d)
				add(j,mu[j/d]*g[d]);
			p++;
		}
		int Newn=min(query[i].n,query[i].m);
		int ans=0;
		for(int l=1;l<=Newn;){
			int r=min(Newn,min(get_next(query[i].n,l),get_next(query[i].m,l)));
			ans+=(ask(r)-ask(l-1))*(query[i].n/l)*(query[i].m/l);
			l=r+1;
		}
		ANS[query[i].id]=ans;
	}
	for(int i=1;i<=q;i++)
		printf("%lld\n",ANS[i]%(1ll<<31));
	return 0;
}