1057. Amount of Degrees
Time limit: 1.0 second
Memory limit: 64 MB
Memory limit: 64 MB
Create a code to determine the amount of integers, lying in the set [X;Y] and being a sum of exactlyK different integer degrees of B.
Example. Let X=15, Y=20, K=2, B=2. By this example 3 numbers are the sum of exactly two integer degrees of number 2:
17 = 24+20,
18 = 24+21,
20 = 24+22.
18 = 24+21,
20 = 24+22.
Input
The first line of input contains integers X and Y, separated with a space (1 ≤ X ≤ Y ≤ 231−1). The next two lines contain integers K and B (1 ≤ K ≤ 20; 2 ≤ B ≤ 10).
Output
Output should contain a single integer — the amount of integers, lying between X and Y, being a sum of exactly K different integer degrees of B.
Sample
input | output |
---|---|
15 20
2
2
| 3
|
题目大意:
求给定区间[X,Y]中满足下列条件的整数个数:
这个数恰好等于K个互不相等的B的整数次幂之和。例如,设X=15,Y=20,K=2,B=2,则有且仅有下列三个数满足题意:
17 = 2^4+2^0,
18 = 2^4+2^1,
20 = 2^4+2^2.
1 ≤ X ≤ Y ≤ 2^31−1,1 ≤ K ≤ 20, 2 ≤ B ≤ 10。
实质上:求[X,Y] 区间转化为 B 进制 1 的个数为K 的数的出现次数。
题解:数位DP。
注意:只有系数为1的情况,所以每位只能是0,1,至于不相同的只要是取1 的时候不不相同,至于0是不变的。
大体思路:
所求的数为互不相等的幂之和,亦即其B进制表示的各位数字都只能是0和1。
因此,我们只需讨论二进制的情况,其他进制都可以转化为二进制求解。
本题区间满足区间减法,因此可以进一步简化问题:令count[i..j]表示[i..j]区间内合法数的个数,则
count[i..j]=count[0..j]-count[0..i-1]。
换句话说,给定n,我们只需求出从0到n有多少个符合条件的数。
步骤:
首先预处理f。
f[i,j]代表i位二进制数中恰好有j个1的数的个数。
f[i,j]=f[i-1,j]+f[i-1,j-1]
计算count[0..n]
像前几题一样,一位一位枚举,只需要多记录后面需要的1的个数即可。
if digit[i] = 1 then ans = ans + f[i,need]
need就是后面需要的1的个数。
#include <stdio.h>
#include <cstring>
#include <vector>
#include <algorithm>
#include <string.h>
using namespace std;
int digit[33];
typedef long long LL;
const int N = 50;
int f[N][N]; //f[i][j]表示前i个中选 j 个 1的个数
void init()
{
f[0][0] = 1;
for(int i=1;i<33;i++)
{
f[i][0] = f[i-1][0];
for(int j=1;j<=i;j++)
{
f[i][j] = f[i-1][j] + f[i-1][j-1];
}
}
}
int solve(int x,int k,int B)
{
int len=-1;
//vector<int> v;
while(x)
{
digit[++len]=x%B;
//v.push_back(x%B);
x/=B;
}
int cnt = 0,ans = 0;
for(int i=len;i>=0;i--)
{
if(digit[i]==1) // 如果为 1,则依次求解
{
ans+=f[i][k-cnt]; //need
cnt++;
if(cnt==k)
break;
}
else if(digit[i]>1) //假如大于1的话,相当于所有的位可以为 1,所以直接求解跳出
{
ans += f[i+1][k-cnt];
break;
}
}
if(cnt==k)
ans++;
return ans;
}
int main()
{
init();
int x,y,k,B;
while(~scanf("%d%d%d%d",&x,&y,&k,&B))
{
printf("%d\n",solve(y,k,B)-solve(x-1,k,B));
}
return 0;
}