「SDOI2014」数表

「SDOI2014」数表
题解参看dalao
个人认为这一道题的公式不是很难推(注意有两项相乘的将乘积作为一个变量T来枚举，再枚举T的约数d,这样就有可能能够预处理)
我认为该题的最大的难点(也是坑死我的地方)就是小于等于a这一个限制。 看了dalao们的题解和代码后明白了其实只需要先对f排序,再将查询离线做(按a从小到大排序)用树状数组维护即可。
对于每一个查询的a来说对f[d]小于等于a的修改.
修改要改的是a,2*a,3*a…….;

代码

#include <iostream>
#include <cstdio>
#include <algorithm>
using namespace std;
typedef long long ll;
const int maxn=100000+5;
int tot,prime[maxn],mu[maxn];
ll f[maxn],minv[maxn],Mod;
int t,p[maxn],pos;
int tree[maxn],ans[maxn];
bool used[maxn];
struct node
{
    int n,m,id; ll a;
    bool operator < (const node &o) const
    {
        return a<o.a;
    }
}q[maxn];
bool cmp(int x,int y)
{
    return f[x]<f[y];
}
void get()
{
    f[1]=mu[1]=1; ll y;
    for(int i=2; i<=100000; i++)
    {
        if(!used[i])
        {
            prime[++tot]=i; mu[i]=-1; minv[i]=i;
            f[i]=i+1;
        }
        for(int j=1; j<=tot && prime[j]*i<=100000; j++)
        {
            y=prime[j]*i; used[y]=true;
            if(i%prime[j]==0)
            {
                mu[y]=0; f[y]=f[i]/((1-minv[i]*prime[j])/(1-prime[j]));
                minv[y]=minv[i]*prime[j]; f[y]*=((1-minv[y]*prime[j])/(1-prime[j])); break;
            }
            else
            {
                minv[y]=prime[j];
                f[y]=f[i]*f[prime[j]]; mu[y]=-mu[i];
            }
        }
    }
}
void modify(int k,ll x)
{
    while(k<=100000)
    {
        tree[k]+=x; k+=(k&(-k));
    }
}
ll query(int k)
{
    int ans1=0;
    while(k)
    {
        ans1+=tree[k];k-=(k&(-k));
    }
    return ans1;
}
int main()
{
    freopen("test.in","r",stdin);
    freopen("test.out","w",stdout);
    Mod=1LL<<(31);  get();
    scanf("%d",&t);
    for(int i=1; i<=100000; i++)
    {
        p[i]=i; //cout<<f[i]<<endl;
    }
    for(int i=1; i<=t; i++)
    {
        int n,m,a;scanf("%d %d %d",&n,&m,&a);
        if(n>m) swap(n,m);
        q[i].n=n; q[i].m=m; q[i].a=a; q[i].id=i;
    }
    sort(q+1,q+t+1);
    sort(p+1,p+100001,cmp);
    pos=1;
    for(int i=1; i<=t; i++)
    {
        while(f[p[pos]]<=q[i].a)
        {
            for(int j=p[pos]; j<=100000; j+=p[pos])
            {
                modify(j,f[p[pos]]*mu[j/p[pos]]);   
            }
            pos++;
        }
        int last=1;
        while(last<=q[i].n)
        {
            int p=min(q[i].n/(q[i].n/last),q[i].m/(q[i].m/last));
            ans[q[i].id]+=(q[i].n/last)*(q[i].m/last)*((query(p)-query(last-1)));
            last=p+1;
        }
    }
    for(int i=1; i<=t; i++)
    {
        printf("%d\n",(ans[i]+Mod)%Mod);
    }
    return 0;
}