HDU 3208 Integer’s Power 【容斥】

题目:点击打开链接

题意:给定一个区间,求出区间内所有数的最大的指数和。

分析:可以先求出1~n中有的区间指数和,再取端点相减。考虑到最大的指数不超过60(2^60>1e18),因此可以枚举所有指数,根据指数k对n开根,因为用a=pow(n,i)可能会产生不超过1的精度误差,所以可以每次判断a-1,a和a+1哪个值的k次方不超过n且更接近,并且a+1的k次方可能会超过long long表示的最大范围,要特殊处理。对于每个指数k,得到的a的值是预处理出来的小于等于a^k的个数,因为一个指数t1若是能被另一个指数t2整除,那么t2对应的a个数都在t1中出现过,相同的数会重复计算,所以在t1对应的a个数减去t2对应的a个数,所以指数从大到小枚举,会遇到上述情况,类似dp的思想,处理出所有的a后,乘上各自的指数即可。

#pragma comment(linker, "/STACK:1024000000,1024000000")
#include <stdio.h>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <queue>
#include <set>
#include <vector>
#include <map>
#define PR pair<int,int>
#define MP make_pair
#define fi first
#define se second
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define sqr(x) ((x)*(x))
#define ll long long
const ll INF = 1e18+300;
const int inf=0x3f3f3f3f;
const int M=100010;
const int N=100100;
const ll MOD=1000000007;
const double eps=1e-3;
const double pi=acos(-1.0);
using namespace std;
ll l,r;
ll dp[63];
ll qpow(ll x,int n)
{
    ll ans=1;
    while(n)
    {
        if(n&1)
        {
            double tmp=1.0*INF/ans;
            if(x>tmp) return -1;
            ans*=x;
        }
        n>>=1;
        if(x>((ll)1<<31)&&n) return -1;
        x=x*x;
    }
    return ans;
}
ll change(ll x,ll n)
{
    ll a=(ll)pow(x,1.0/n);
    ll qt=qpow(a,n);
    ll ql=qpow(a-1,n);
    ll qr=qpow(a+1,n);
    if(qr!=-1&&qr<=x) return a+1;
    else if(qt!=-1&&qt<=x) return a;
    else if(ql!=-1&&ql<=x) return a-1;
}
ll calc(ll n)
{
    if(n==1) return 1;
    int i,j;
    memset(dp,0,sizeof(dp));
    dp[1]=n;
    for(i=2;i<63;i++)
    {
        dp[i]=change(n,i)-1;
        if(dp[i]==0) break;
    }
    int k=i;
    for(i=k-1;i>0;i--)
        for(j=1;j<i;j++)
            if(i%j==0) dp[j]-=dp[i];
    ll ans=dp[1];
    for(int i=2;i<k;i++)
        ans+=(i*dp[i]);
    return ans;
}
int main()
{
    while(~scanf("%I64d%I64d",&l,&r))
    {
        if(l==0&&r==0) break;
        printf("%I64d\n",calc(r)-calc(l-1));
    }
}


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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值