CodeForces-431D Random Task
题目描述
求一个\(n\),使得\(n+1\)到\(2n\)这些数的二进制中恰好有\(k\)个\(1\)的数有\(m\)个。
Input
输入包含两个正整数\(m,k\)。$(0<=m<=1e18, 1<=k<=64) $
Output
输出\(n\)。
Sample Input
1 1
3 2
Sample Output
1
5
首先我们要知道这个答案是具有单调性的。
即若满足\(m>n\),则区间\([m,2m]\)内满足条件的数的数量\(>\)区间\([n,2n]\)内满足条件的数的数量。
证明如下:
将区间\([n+1,2n]\)分割成两个区间\([n+1,n+1]\)和\([n+2,2n]\)。
而区间\([n+2,2*(n+1)]\)可以分割成\([n+2,2n]\)和\([2n+1,2n+2]\)。
而\(2n+2\)二进制中\(1\)的数量和\(n+1\)中\(1\)的数量是相同的。
而多了一个\(2n+1\)只有可能让数量增加,所以答案满足单调性。
即对于\(m>n\),则区间\([m,2m]\)内满足条件的数的数量\(>\)区间\([n,2n]\)内满足条件的数的数量。
证毕。
于是我们就可以二分求解了。
现在我们就只用判断\([n+1,2n]\)中满足条件的数的数量即可。
我们发现对于区间\([n+1,2n]\)满足条件的数的数量可以转换为\([1,2n]\)-\([1,n]\)的数量。
而一个前缀的数量肯定会更好求解。
问题进一步转换为求解\([1,n]\)内满足条件数的数量。
这个问题我们可以利用两种方法求解:
方法一:计数DP
对于当前\(n\),我们从二进制位高的开始考虑,首先我们只有在遇到一个\(1\)时才能更新答案。
对于当前遇到的\(1\),我们发现前面的位上的数是固定的,一共有\(tot\)个\(1\),而我们再后面的\(i\)个位上可以选择\(k-tot\)个位置放置\(1\),所以答案增加\(C(i,n-tot)\),然后\(tot++\)即可。
代码如下
#include <bits/stdc++.h>
using namespace std;
#define int long long
#define reg register
#define Raed Read
#define clr(a,b) memset(a,b,sizeof a)
#define debug(x) cerr<<#x<<" = "<<x<<endl;
#define rep(a,b,c) for(reg int a=(b),a##_end_=(c); a<=a##_end_; ++a)
#define ret(a,b,c) for(reg int a=(b),a##_end_=(c); a<a##_end_; ++a)
#define drep(a,b,c) for(reg int a=(b),a##_end_=(c); a>=a##_end_; --a)
#define erep(i,G,x) for(int i=(G).Head[x]; i; i=(G).Nxt[i])
#pragma GCC target("avx,avx2,sse4.2")
#pragma GCC optimize(3)
inline int Read(void) {
int res=0,f=1;
char c;
while(c=getchar(),c<48||c>57)if(c=='-')f=0;
do res=(res<<3)+(res<<1)+(c^48);
while(c=getchar(),c>=48&&c<=57);
return f?res:-res;
}
template<class T>inline bool Min(T &a, T const&b) {
return a>b?a=b,1:0;
}
template<class T>inline bool Max(T &a, T const&b) {
return a<b?a=b,1:0;
}
const int N=1e5+5,M=2e5+5,mod=1e9+7;
bool MOP1;
int m,k,C[70][70];
int solve(int x) {
int Ans=0,tot=0;
drep(i,63,0)if(x&(1ll<<i)) {
if(tot<=k)Ans+=C[i][k-tot];
tot++;
}
return Ans;
}
bool MOP2;
inline void _main() {
m=Read(),k=Read();
int L=1,R=1e18,Ans=0;
C[0][0]=1;
rep(i,1,64) {
C[i][0]=1;
rep(j,1,i)C[i][j]=C[i-1][j-1]+C[i-1][j];
}
while(L<=R) {
int mid=(L+R)>>1;
if(solve(mid*2)-solve(mid)>=m)Ans=mid,R=mid-1;
else L=mid+1;
}
printf("%lld\n",Ans);
}
signed main() {
_main();
return 0;
}
方法二:数位DP
我们可以直接将\(n\)转换为二进制数,那么就可以直接\(dp\)了。
\(dp[i][j]\)代表考虑到第\(i\)位,\(1\)的数量为\(j\)的数的个数。
代码如下
#include <bits/stdc++.h>
using namespace std;
#define int long long
#define reg register
#define Raed Read
#define clr(a,b) memset(a,b,sizeof a)
#define debug(x) cerr<<#x<<" = "<<x<<endl;
#define rep(a,b,c) for(reg int a=(b),a##_end_=(c); a<=a##_end_; ++a)
#define ret(a,b,c) for(reg int a=(b),a##_end_=(c); a<a##_end_; ++a)
#define drep(a,b,c) for(reg int a=(b),a##_end_=(c); a>=a##_end_; --a)
#define erep(i,G,x) for(int i=(G).Head[x]; i; i=(G).Nxt[i])
#pragma GCC target("avx,avx2,sse4.2")
#pragma GCC optimize(3)
inline int Read(void) {
int res=0,f=1;
char c;
while(c=getchar(),c<48||c>57)if(c=='-')f=0;
do res=(res<<3)+(res<<1)+(c^48);
while(c=getchar(),c>=48&&c<=57);
return f?res:-res;
}
template<class T>inline bool Min(T &a, T const&b) {
return a>b?a=b,1:0;
}
template<class T>inline bool Max(T &a, T const&b) {
return a<b?a=b,1:0;
}
const int N=1e5+5,M=2e5+5,mod=1e9+7;
bool MOP1;
int m,k,dp[70][70],A[70];
int dfs(int pos,int tot,int limit) {
if(!pos)return tot==k;
if(!limit&&dp[pos][tot]!=-1)return dp[pos][tot];
int up=limit?A[pos]:1,res=0;
rep(i,0,up)res+=dfs(pos-1,tot+i,limit&(i==up));
if(!limit)dp[pos][tot]=res;
return res;
}
int solve(int x) {
int len=0;
while(x)A[++len]=x&1,x>>=1;
return dfs(len,0,1);
}
bool MOP2;
inline void _main() {
// cerr<<"M="<<(&MOP2-&MOP1)/1024.0/1024.0<<endl;
m=Read(),k=Read();
int L=1,R=1e18,Ans=0;
memset(dp,-1,sizeof dp);
while(L<=R) {
int mid=(L+R)>>1;
if(solve(mid*2)-solve(mid)>=m)Ans=mid,R=mid-1;
else L=mid+1;
}
printf("%lld\n",Ans);
}
signed main() {
_main();
return 0;
}
CodeForces-431D Random Task的更多相关文章
- Codeforces Round #247 (Div. 2) D. Random Task
D. Random Task time limit per test 1 second memory limit per test 256 megabytes input standard input ...
- codeforces 478B Random Teams
codeforces 478B Random Teams 解题报告 题目链接:cm.hust.edu.cn/vjudge/contest/view.action?cid=88890#probl ...
- codeforces 431 D. Random Task 组合数学
题意: 给定m,k 0 <= m <= 10^18 ,1 <= k <= 64 求一个数n,满足n+1,n+2,...n+n这n个数中,刚好有m个数的2进制表示法刚好有k个1 ...
- codeforces 478B Random Teams 解题报告
题目链接:http://codeforces.com/problemset/problem/478/B 题目意思:有 n 个人,需要将这班人分成 m 个 组,每个组至少含有一个人,同一个组里的人两两可 ...
- Codeforces 653F Paper task SA
Paper task 如果不要求本质不同直接st表二分找出最右端, 然后计数就好了. 要求本质不同, 先求个sa, 然后用lcp求本质不同就好啦. #include<bits/stdc++.h& ...
- CF431D Random Task 二分+数位dp
One day, after a difficult lecture a diligent student Sasha saw a graffitied desk in the classroom. ...
- Codeforces 846F - Random Query
原题链接:http://codeforces.com/contest/846/problem/F 题意:给一个数列,任意取区间[l, r],问区间内不同数字的个数的期望是多少. 思路: 对于第i个数a ...
- Codeforces 1067E - Random Forest Rank(找性质+树形 dp)
Codeforces 题面传送门 & 洛谷题面传送门 一道不知道能不能算上自己 AC 的 D1E(?) 挺有意思的结论题,结论倒是自己猜出来了,可根本不会证( 开始搬运题解 ing: 碰到这样 ...
- Codeforces 788E - New task(线段树)
Codeforces 题目传送门 & 洛谷题目传送门 这是一道 *2900 的 D1E,而且被!我!自!己!搞!出!来!了! 虽然我承认它难度及摆放的位置异常异常虚高,并且就算我到了现场也不可 ...
- CodeForces - 846F Random Query(期望)
You are given an array a consisting of n positive integers. You pick two integer numbers l and r fro ...
随机推荐
- Angular 文档中链接的修改路径
在 Angular 文档程序中的左侧链接的修改路径在哪里? 如下图所示的修改路径. 左侧链接的修改路径在 angular-cn\aio\content\navigation.json 这个文件中. 你 ...
- 4.JSP内置对象
JSP内置对象,JSP提供了由容器实现和管理的内置对象,也可以称之为隐含对象,这些内置对象不需要通过 JSP页面编写来实例化,在所有的JSP页面中都可以直接使用,它起到了简化页面的作用. 在JSP中一 ...
- hive on tez 任务失败
最近再hue 集群查询任务经常失败,经过几天的观察,终于找到原因,报错如下 Status: FailedVertex failed, vertexName=Map 1, vertexId=vertex ...
- conda install -c anaconda
有些包在conda默认的channels中不包含,比如cudatoolkit-8.0,cudnn等,这时只需要在conda install指令后加上-c anaconda即可.比如要下载cudatoo ...
- redis学习之路
cp /home/zzq/桌面/M+.txt /home/zzq/mySercet#将桌面文件下的M+.txt 复制到mySercetrm -rf M+.txt#强制递归删除M+.txt1.网络适配器 ...
- vscode + php+ftp
首先,php网站的文件都整理到一个文件夹中: 然后,用vscode的File.Open Folder打开刚才的文件夹: 3,Ctrl+Shift+P,输入SFTP:Config,会打开一个配置文件,编 ...
- Android 面试汇总<三>
1.3 计算机网络 基础 Q:五层协议的体系结构分别是什么?每一层都有哪些协议? 技术点:网络模型.协议 思路:分条解释每层名字以及协议 参考回答: 物理层 数据链路层:逻辑链路控制LLC.媒体接入控 ...
- 利用Python处理向不受信任主机发送请求
,HTTP vs HTTPS 超文本传输协议HTTP协议被用于在Web浏览器和网站服务器之间传递信息.HTTP协议以明文方式发送内容,不提供任何方式的数据加密,如果攻击者截取了Web浏览器和网站服务器 ...
- apache用户认证
创建一个目录abc:mkdir abc在此目录下建一个文件:12.txt正常情况下可以访问. 建立用户认证,从而使用户访问特定目录文件需要认证 在虚拟主机配置文件中即vim /usr/local/ap ...
- Redis 常用命令 大全
Redis 常用命令 发现几个很好的 Redis 常用命令汇总大全网页,分享给小伙伴们~ 1.Redis 命令参考 http://redisdoc.com/string/index.html 2.W3 ...