BZOJ2440(全然平方数)二分+莫比乌斯容斥

BZOJ2440(全然平方数)二分+莫比乌斯容斥

大家好,又见面了,我是全栈君,祝每个程序员都可以多学几门语言。

题意:全然平方数是指含有平方数因子的数。求第ki个非全然平方数。


解法:比較明显的二分,getsum(int middle)求1-middle有多少个非全然平方数,然后二分。求1-middle的非全然平方数个数能够用总数减掉全然平方数个数。计算全然平方数的个数用容斥:

    首先加上n/(2*2)+n/(3*3)+n/(5*5)+n/(7*7)…+…然后减掉出现两次的,然后加上三次的…奇加偶减。这就是mou的原型,用mou数组计算非常easy;

       BZOJ2440(全然平方数)二分+莫比乌斯容斥

代码:

/******************************************************
* @author:xiefubao
*******************************************************/
#pragma comment(linker, "/STACK:102400000,102400000")
#include <iostream>
#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <queue>
#include <vector>
#include <algorithm>
#include <cmath>
#include <map>
#include <set>
#include <stack>
#include <string.h>
//freopen ("in.txt" , "r" , stdin);
using namespace std;

#define eps 1e-8
#define zero(_) (abs(_)<=eps)
const double pi=acos(-1.0);
typedef unsigned long long LL;
const int Max=100010;
const LL INF=2e16+7;

LL mou[Max];
void init()
{
    for(LL i=2; i<Max; i++)
    {
        if(!mou[i])
        {
            mou[i]=i;
            for(LL j=i*i; j<Max; j+=i)
                mou[j]=i;
        }
    }
    mou[1]=1;
    for(int i=2; i<Max; i++)
    {
        if(i/mou[i]%mou[i]==0) mou[i]=0;
        else mou[i]=-mou[i/mou[i]];
    }
}
int k;
LL getnum(LL middle)
{
    LL ans=0;
    for(LL i=1; i*i<=middle; i++)
    {
        ans+=mou[i]*(middle/(i*i));
    }
    return ans;
}
int main()
{
    init();
    int t;
    cin>>t;
    while(t--)
    {
        scanf("%d",&k);
        LL left=1,right=INF;
        while(left<=right)
        {
            int middle=(left+right)/2;
            if(getnum(middle)<k)
                left=middle+1;
            else
                right=middle-1;
        }
        cout<<left<<'\n';
    }
    return 0;
}

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请联系我们举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/118488.html原文链接:https://javaforall.net

(0)
上一篇 2021年12月5日 下午2:00
下一篇 2021年12月5日 下午2:00


相关推荐

  • QT之Android下获取手机传感器数据学习笔记

    QT之Android下获取手机传感器数据学习笔记QT+=coreguisensorspositioning其中sensors是获取手机上传感器数据的组件,positioning是获取位置信息的组件1、获取陀螺仪传感器数据#include&lt;QGyroscope&gt;QGyroscope*gyroscope;QGyroscopeReading*reader;gyroscope=newQGyro…

    2022年7月13日
    18
  • 5G基站:宏基站&&微基站&&皮基站&&飞基站

    5G基站:宏基站&&微基站&&皮基站&&飞基站四月来临 春暖花开 终于慢慢滴都开始恢复正常了 5G 基站分为宏基站和微小基站两种 宏基站主要用于室外覆盖 微小基站发射功率较小 主要用于室内场景 根据 3GPP 组织的规则 无线基站分为 4 类 分别是宏基站 微基站 皮基站和飞基站 划分基站主要依据是功率和容量 宏基站的功率在 10W 以上 覆盖能力 覆盖半径 在 200 米以上 可同时接入用户数视基站规模而定 一般在 1000 个以上 微基站功率为 500mW

    2026年3月19日
    2
  • 运行怎么进入文件路径_cmd命令怎么进入某个文件夹

    运行怎么进入文件路径_cmd命令怎么进入某个文件夹1.通过Windows+R进入命令调出运行2.输入cmd进入命令窗口(默认的一般是c:\Users下的某个文件夹,例如我的是c:\Users\LML)3.若想进入c盘的其他文件路径下,可以通过在目录下输入cd..进入上一层目录,直到进入c盘根目录;通过命令行输入c:\cd+文件或文件夹路径 进入目标文件夹4.若想进入其他盘下的文件路径,通过在命令行默认路径后输入想进入的盘名加上冒号,例如:c:…

    2022年10月15日
    4
  • NSGA2 Python实现

    NSGA2 Python实现#importingnecessarymodulesimportmathimportrandomimportmatplotlib.pyplotasPlt#FirstFunctiontooptimizedeffunction1(x1,x2):value=-x1*2+x2returnvalue#SecondFunctiontooptimizedeffunction2(x1,x2):value=-x1*5*x2

    2022年5月12日
    46
  • JS 实现替换字符串中所有指定字符总结

    JS 实现替换字符串中所有指定字符总结最近在写前端需要把字符串中的空格全部替换掉 对 js 不是很了解 现在对此进行总结 letstr console log str replace 如果按照上面的写法将会打印 而如果按照如下写法 则或实现此功能 letstr console log str replace g 将会打印 所以 可以总结一下 str str replace 你想要替

    2026年3月19日
    2
  • poj 2689 巧妙地运用素数筛选

    poj 2689 巧妙地运用素数筛选

    2022年1月8日
    46

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号