D - Together Square

传送门

D - Together Square

题意:问有多少对(i,j)满足i<=n&&j<=n&&i*j是一个平方数

分析打表然后OEIS。考虑枚举i,再加上j可能的个数。对于一个i,如果他是平方数,那么j的取值也必须是平方数;如果i不是平方数,那么找到大于i的最小平方数,也就是把i的质因数中出现奇数次的质因子补上,我们把vx[i]定义为i的质因数中出现奇数次的质因子的乘积,也就是j最小的取值。之后只需加上平方数的个数即可。

代码

#include<bits/stdc++.h>
using namespace std;
#define int long long
const int N=2e5+5;
vector<int>v[N];
int vx[N];

void init()//O(能过1e6)
{
    for(int i=2;i<N;i++)//枚举质因数 
    {
        if(v[i].size()==0)//质数
        {
            for(int j=i;j<N;j+=i)
            {
                int temp=j;
                while(temp%i==0)
                {
                    v[j].push_back(i);
                    temp/=i;
                }
            }
        }
    }
}

signed main()
{
    init();
    int n;
    cin>>n;
    for(int i=1;i<=n;i++)
    {
        vx[i]=1;
        for(int j=0;j<v[i].size();j++)
        {
            if(v[i][j]==v[i][j+1]&&j+1<v[i].size())
            {
                j++;
            }
            else vx[i]*=v[i][j];
        }
    }
    int sum=0;
    for(int i=1;i<=n;i++)
    {
        sum+=sqrt(n/vx[i]);
    }
    cout<<sum<<endl;
}

  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

chmpy

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值