BZOJ 2693 jzptab 莫比乌斯反演

题目大意：同2154 多组数据

后面那坨东西 由于积性函数的约数和仍是积性函数 因此只需要线性筛一下就行

i%prime[j]==0那部分由于多出来的因数都不是无平方因子数因此μ值都为0 增加的只有原先的D/i

#include 
#include 
#include 
#include 
#define M 10001000
#define MOD 100000009
using namespace std;
typedef long long ll;
int prime[1001001],tot;
bool not_prime[M];
ll h[M],sum[M];
void Linear_Shaker()
{
	int i,j;
	h[1]=1;
	for(i=2;i>1)%MOD;
	ll re2=(y*(y+1)>>1)%MOD;
	return re1*re2%MOD;
}
int Query(int n,int m)
{
	int i,last,re=0;
	if(n>m) swap(n,m);
	for(i=1;i<=n;i=last+1)
	{
		last=min(n/(n/i),m/(m/i));
		re+=Sum(n/i,m/i)*(sum[last]-sum[i-1])%MOD;
		re%=MOD;
	}
	return (re+MOD)%MOD;
}
int main()
{
	int T,n,m;
	Linear_Shaker();
	for(cin>>T;T;T--)
	{
		scanf("%d%d",&n,&m);
		printf("%d\n", Query(n,m) );
	}
}