CQOI2016 伪光滑数

这里讲讲搜索。

二分答案ans,若要求ans~N中伪光滑数的个数,可以用1~N的个数减去1~ans-1的个数。搜索x范围内伪光滑数的个数,先枚举伪光滑数中最大的质因子,就可以知道伪光滑数因子个数的上限,以及因数大小的上限,可以设计dfs函数dfs(num,x,last)表示伪光滑数因子个数的上限,伪光滑数大小上限,上次用过的因子。直接枚举所用的因子转移状态,其中last不递增防止计算重复。

优化就是若last^num<=x,意味着剩下的因子可以在last之内随便选,利用挡板原理那么方案数就是C(last以内因子个数,last以内因子个数+num),其中不一定要选够num个因子,所以相当于多了一个可以选的因子。可以卡过。

  第二个优化就是当x比较小时记忆化一下,优化明显。

这个方法可以通过k为任意值的数据,N=10^18时,k最大为5百万左右,附上代码。

#include<cstdio>
#include<cstring>
#include<algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
typedef long long ll;
ll N;
ll p[200],top;
const ll inf=3000000000000000000;
ll lim[200][200],c[200][200];
void prepare()
{
    fo(i,2,128)
    {
        int ok=1;
        fo(j,1,top)
        if (i%p[j]==0) ok=0;
        if (ok) p[++top]=i;
    }
    c[0][0]=1;
    fo(i,1,100)
        fo(j,0,i)
        c[j][i]=(j>0?c[j-1][i-1]:0)+c[j][i-1];
}
ll f[70][3000][40];
ll dfs(ll cs,ll n,ll last)
{
    if (n<3000 && f[cs][n][last]) return f[cs][n][last];
    if (n==0) return 0;
    if (cs==0) return 1;
    if (lim[last][cs]<=n) return c[cs][last+cs];
    ll re=1;
    fd(i,last,1) re+=dfs(cs-1,n/p[i],i);
    if (n<3000) f[cs][n][last]=re;
    return re;
}
ll get(ll n)
{
    ll re=0;
    fo(i,1,top) 
    {
        lim[i][1]=p[i];
        fo(j,2,70)
        {
            lim[i][j]=lim[i][j-1]*p[i];
            if (lim[i][j]/p[i]!=lim[i][j-1]) lim[i][j]=inf;
        }
    }
    fo(i,1,top)
    {
        ll cs=0,tp=N;
        while (tp>=p[i]) cs++,tp/=p[i];
        ll num=cs-1;
        re+=dfs(num,n/p[i],i);
    }
    return re;
}
int main()
{
    prepare();
    ll k;
    scanf("%lld%lld",&N,&k);
    ll l=1,r=N,srm=get(N),ans=0;
    while (l<=r)
    {
        ll mid=(l+r)>>1;
        if (srm-get(mid-1)>=k) ans=mid,l=mid+1; else
        r=mid-1;
    }
    printf("%lld",ans);
}
相关推荐
程序员的必经之路! 【限时优惠】 现在下单,还享四重好礼: 1、教学课件免费下载 2、课程案例代码免费下载 3、专属VIP学员群免费答疑 4、下单还送800元编程大礼包 【超实用课程内容】  根据《2019-2020年中国开发者调查报告》显示,超83%的开发者都在使用MySQL据库。使用量大同时,掌握MySQL早已是运维、DBA的必备技能,甚至部分IT开发岗位也要求对据库使用和原理有深入的了解和掌握。 学习编程,你可能会犹豫选择 C++ 还是 Java;入门据科学,你可能会纠结于选择 Python 还是 R;但无论如何, MySQL 都是 IT 从业人员不可或缺的技能!   套餐中一共包含2门MySQL据库必学的核心课程(共98课时)   课程1:《MySQL据库从入门到实战应用》   课程2:《高性能MySQL实战课》   【哪些人适合学习这门课程?】  1)平时只接触了语言基础,并未学习任何据库知识的人;  2)对MySQL掌握程度薄弱的人,课程可以让你更好发挥MySQL最佳性能; 3)想修炼更好的MySQL内功,工作中遇到高并发场景可以游刃有余; 4)被面试官打破沙锅问到底的问题问到怀疑人生的应聘者。 【课程主要讲哪些内容?】 课程一:《MySQL据库从入门到实战应用》 主要从基础篇,SQL语言篇、MySQL进阶篇三个角度展开讲解,帮助大家更加高效的管理MySQL据库。 课程二:《高性能MySQL实战课》主要从高可用篇、MySQL8.0新特性篇,性能优化篇,面试篇四个角度展开讲解,帮助大家发挥MySQL的最佳性能的优化方法,掌握如何处理海量业务据和高并发请求 【你能收获到什么?】  1.基础再提高,针对MySQL核心知识点学透,用对; 2.能力再提高,日常工作中的代码换新貌,不怕问题; 3.面试再加分,巴不得面试官打破沙锅问到底,竞争力MAX。 【课程如何观看?】  1、登录CSDN学院 APP 在我的课程中进行学习; 2、移动端:CSDN 学院APP(注意不是CSDN APP哦)  本课程为录播课,课程永久有效观看时长 【资料开放】 课件、课程案例代码完全开放给你,你可以根据所学知识,自行修改、优化。  下载方式:电脑登录课程观看页面,点击右侧课件,可进行课程资料的打包下载。
©️2020 CSDN 皮肤主题: 大白 设计师:CSDN官方博客 返回首页