Codefoces 955C 题解

本文详细解析了Codeforces 955C题目,通过题意简述和数据介绍,阐述了解题思路。主要讨论了当p>=3和p=2时的不同处理方法,利用数学性质优化了暴力求解的过程,最终给出代码实现。

题意简述

多组数据。每次给定lll,rrr(l,rl,rl,r&lt;=1e18&lt;=1e18<=1e18),求lll~rrr中有多少xxx满足x=apx=a^px=ap,其中aaappp是两个整数满足a&gt;0,p&gt;1a&gt;0,p&gt;1a>0,p>1

数据

输入:
6
1 4
9 9
5 7
12 29
137 591
1 1000000

输出:
2
1
0
3
17
1111

思路

会发现,当p&gt;=3p&gt;=3p>=3时,使得0&lt;=ap&lt;=1e180&lt;=a^p&lt;=1e180<=ap<=1e18aaa并不多,p=3p=3p=3大概有1e61e61e6,然后后面更少。经计算,总共的个数(去重后)是100333210033321003332,就比1e61e61e6多一点。考虑暴力存下,排序去重后存着。p&gt;=3p&gt;=3p>=3的情况只要upperbound−lowerboundupperbound-lowerboundupperboundlowerbound即可(lower_bound和upper_bound返回的都是指针,如果相减就珂以得到中间长度了。不用+1+1+1的原因是upper_bound找到的那个不珂以取)。记为a1a_1a1

那如何处理p=2p=2p=2的情况?

相当于求有多少aaa满足l&lt;=a2&lt;=rl&lt;=a^2&lt;=rl<=a2<=r

我们会发现,如果设S(n)S(n)S(n)为满足1&lt;=a2&lt;=n1&lt;=a^2&lt;=n1<=a2<=naaa个数,则我们要求的就是S(r)−S(l−1)S(r)-S(l-1)S(r)S(l1)

由数学常识得,S(n)=nS(n)=\sqrt{n}S(n)=n

故我们要求的就是r−l−1\sqrt{r}-\sqrt{l-1}rl1。记为a2a_2a2

然后只要把a1a_1a1a2a_2a2加一起即可。

代码:

#include<bits/stdc++.h>
#define int long long
using namespace std;
const int Max=1000000ll*1000000ll*1000000ll;//1e18

int Root(int x)
{
    return (int)sqrt(x*1.0+0.5);
}
vector<int>tmp;
void Build()
{
    for(int i=2;i<=1000000;i++)//枚举底数
    {
        int s=i*i*i;//至少三次方
        while(s<=Max)
        {
            int rs=Root(s);
            if (rs*rs<s) tmp.push_back(s);
            //这个就是判s是否是平方数,为了避免和后面考虑p=2的情况重合

            if (s>Max/i) break;//如果爆了赶快退出
            else s*=i;//没爆继续乘
        }
    }
    sort(tmp.begin(),tmp.end());
    tmp.erase(unique(tmp.begin(),tmp.end()),tmp.end());//排序去重,方便lower_bound和upper_bound
    //printf("%d\n",tmp.size());
    //=1003332,如果要自己写(即不完全Copy我的),珂以先对一下这个答案
}

int l,r;
void Solve()
{
    int a1=upper_bound(tmp.begin(),tmp.end(),r)-lower_bound(tmp.begin(),tmp.end(),l);
    int a2=Root(r)-Root(l-1);
    printf("%I64d\n",a1+a2);//上面说的a1和a2
}

main()
{
    Build();
    int T;scanf("%I64d",&T);
    while(T--)
    {
        scanf("%I64d%I64d",&l,&r);
        Solve();
    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值