感觉看着别人推的式子看了一万年才大概也许看懂了!?
QAQ感觉这一篇写的挺详细的就直接贴别人的好了。
然后至于后面的分块。其实我也看了半天才大概也许明白。。
式子最后推成了
(令a为较小值)
∑d=1a(μ(d)∗⌊ad⌋⌊bd⌋)
⌊ad⌋⌊bd⌋ 单调不上升
那么就会取到很多相同的值。
找一下规律。若 a=100 那么当d在[26,33]这段区间内 ⌊ad⌋ 都取3,[34,50]内都取2。
若一段的开头为L,结尾为R,那么令 ⌊aL⌋ 为x,则 R=a/x
那么每次把 a/(a/d) 和 b/(b/d) 取个 min 就好了。
【代码】
#include <cstdio>
#include <iostream>
#include <algorithm>
#include <set>
#define N 50005
using namespace std;
typedef long long ll;
int read()
{
int x=0,f=1;char ch=getchar();
while(!isdigit(ch)){if(ch=='-') f=-1;ch=getchar();}
while(isdigit(ch)){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
return x*f;
}
int T;
int Miu[N],p[5134],sum[N];
bool Not_Prime[N];
void Get_Miu()
{
Miu[1]=1;
for(int i=2;i<N;i++)
{
if(!Not_Prime[i]) p[++p[0]]=i,Miu[i]=-1;
for(int j=1;j<=p[0]&&i*p[j]<N;j++)
{
Not_Prime[i*p[j]]=1;
if(i%p[j]==0) {
Miu[i*p[j]]=0;break;
}
Miu[i*p[j]]=-Miu[i];
}
}
for(int i=1;i<N;i++) sum[i]=sum[i-1]+Miu[i];
}
void Calc(int x,int y)
{
if(x>y) swap(x,y);
ll ans=0;int pos;
for(int i=1;i<=x;i=pos+1)
{
pos=min(x/(x/i),y/(y/i));
ans+=1LL*(sum[pos]-sum[i-1])*(x/i)*(y/i);
}
printf("%lld\n",ans);
}
int main()
{
Get_Miu();
T=read();
while(T--)
{
static int a,b,d;
a=read(),b=read(),d=read();
Calc(a/d,b/d);
}
return 0;
}