bzoj-2693 jzptab

本文探讨了LCM矩阵的求和问题,通过将LCM转换为GCD,利用线性筛技术优化算法效率,最终实现O(n)的时间复杂度。通过预处理和积性函数的性质,进一步简化计算过程。

题意:

有一个n*m的数表,(i,j)位置上的数是LCM(i,j);
求这个数表的和;
n,m<=10^7,多组数据

Crash的数字表格加强版


题解:

考虑把lcm转化成gcd
那答案就是

然后神奇的设:

就有:

一样可以枚举


 的取值,这是O(√n)的;

然后求f(x,y);


大概证明了一下= =

线性筛之后也可以O(√n)求出f(x,y)
总复杂度O(n),常数略大。。

这题显然是卡O(n)过不了呗
那就还得优化


预处理这玩意


然后O(√n)就搞出来啦!



“积性函数的约数和也是积性函数”  ->好像比较显然?
所以g(D)是积性函数
线性筛裸上就好


代码:


#include<stdio.h>
#include<string.h>
#include<algorithm>
#define N 10010000
#define mod 100000009ll
using namespace std;
typedef long long ll;
ll pri[N>>3],tot;
ll g[N],sum[N];
bool vis[N];
void init()
{
	g[1]=sum[1]=1;
	for(ll i=2;i<N;i++)
	{
		if(!vis[i])
		{
			pri[++tot]=i;
			g[i]=((-(ll)(i-1)*i)%mod+mod)%mod;
		}
		for(ll j=1;j<=tot&&i*pri[j]<N;j++)
		{
			vis[i*pri[j]]=1;
			if(i%pri[j]==0)
			{
				g[i*pri[j]]=g[i]*pri[j]%mod;
				break;
			}
			g[i*pri[j]]=g[i]*g[pri[j]]%mod;
		}
		sum[i]=sum[i-1]+g[i];
	}		
}
ll F(ll x,ll y)
{
	return ((ll)(x+1)*x/2)%mod*(((ll)(y+1)*y/2)%mod)%mod;
}
int main()
{
	init();
	ll c,T,n,m,i,j,k;
	ll ans;
	scanf("%lld",&T);
	for(c=1;c<=T;c++)
	{
		scanf("%lld%lld",&n,&m);
		if(n>m)	swap(n,m);
		for(i=1,ans=0;i<=n;i=j+1)
		{
			j=min(n/(n/i),m/(m/i));
			ans+=F(n/i,m/i)*(sum[j]-sum[i-1]+mod)%mod;
			ans%=mod;
		}
		printf("%lld\n",ans);
	}
	return 0;
}



评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值