math

题目描述

简单题

我们观察d(i*j)什么是奇数
i*j是完全平方数时！
我们现在要找√(i*m)中有多少至少包含i每个因子一半的。
为了得到这个i的“一半”，线筛一发就好了。

#include<cstdio>
#include<algorithm>
#include<cmath>
#define fo(i,a,b) for(i=a;i<=b;i++)
using namespace std;
typedef long long ll;
typedef double db;
const int maxn=10000000+10;
int f[maxn],g[maxn],pri[maxn];
bool bz[maxn];
int i,j,k,l,n,top,ans;
ll r,t,m;
void prepare(){
    f[1]=1;
    fo(i,2,maxn-10){
        if (!bz[i]){
            pri[++top]=i;
            g[i]=1;
            f[i]=i;
        }
        fo(j,1,top){
            if (i>(maxn-10)/pri[j]) break;
            bz[i*pri[j]]=1;
            if (i%pri[j]==0){
                f[i*pri[j]]=f[i];
                g[i*pri[j]]=g[i]+1;
                if (g[i*pri[j]]%2) f[i*pri[j]]*=pri[j];
                break;
            }
            f[i*pri[j]]=f[i]*pri[j];
            g[i*pri[j]]=1;
        }
    }
}
ll read(){
    ll x=0;
    char ch=getchar();
    while (ch<'0'||ch>'9') ch=getchar();
    while (ch>='0'&&ch<='9'){
        x=x*10+ch-'0';
        ch=getchar();
    }
    return x;
}
int main(){
    //freopen("t1.in","r",stdin);freopen("t1.out","w",stdout);
    prepare();
    n=read();m=read();
    fo(i,1,n){
        r=floor(sqrt((db)i*m));
        t=r/(ll)f[i];
        if (t%2) ans--;else ans++;
    }
    printf("%d\n",ans);
}