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);
}
根据引用所述,交错序列是一个仅由0和1构成的序列,其中没有相邻的1(可以有相邻的0)。特征值定义为x^ay^b,其中x和y分别表示0和1出现的次。长度为n的交错序列可能有多个。问题要求计算所有长度为n的交错序列特征值的和除以m的余。 根据引用所述,输入文件包含一个行,该行包含三个整n、a、b和m。其中,1≤n≤10000000,0≤a、b≤45,m<100000000。 为了解决这个问题,可以使用动态规划和矩阵快速幂优化的方法,具体实现可以参考引用提到的相关算法。算法的思路是通过计算长度为n的交错序列的特征值,然后将所有特征值求和并对m取余。 具体步骤如下: 1. 使用动态规划计算长度为n的所有交错序列的特征值,将结果保存在一个矩阵中。 2. 使用矩阵快速幂优化,将动态规划的过程进行优化。 3. 对优化后的结果进行求和,并对m取余。 4. 输出结果。 参考引用给出的博客中的代码实现,可以帮助你更好地理解和实现该算法。<span class="em">1</span><span class="em">2</span><span class="em">3</span> #### 引用[.reference_title] - *1* *2* *3* [BZOJ5298 CQOI2018 交错序列 【DP+矩阵快速幂优化】*](https://blog.csdn.net/weixin_30892987/article/details/99470493)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v92^chatsearchT0_1"}}] [.reference_item style="max-width: 100%"] [ .reference_list ]
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值