BZOJ1101: [POI2007]Zap

BZOJ1101

感觉看着别人推的式子看了一万年才大概也许看懂了！？
QAQ感觉这一篇写的挺详细的就直接贴别人的好了。
然后至于后面的分块。其实我也看了半天才大概也许明白。。
式子最后推成了
（令a为较小值）

$$\sum_{d=1}^a{(μ(d)*\lfloor{\frac{a}{d}}\rfloor\lfloor{\frac{b}{d}}\rfloor})$$
$\lfloor{\frac{a}{d}}\rfloor \lfloor{\frac{b}{d}}\rfloor$单调不上升
那么就会取到很多相同的值。
找一下规律。若 $a=100$那么当d在[26,33]这段区间内 $\lfloor{\frac{a}{d}}\rfloor$都取3，[34,50]内都取2。
若一段的开头为L，结尾为R，那么令 $\lfloor{\frac{a}{L}}\rfloor$为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;
}