「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;
}