USACO Stringsobits

最新推荐文章于 2019-06-13 09:24:21 发布

jinTester

最新推荐文章于 2019-06-13 09:24:21 发布

阅读量441

点赞数

分类专栏： USACO 文章标签： USACO

本文链接：https://blog.csdn.net/jinjiahao5299/article/details/39119715

版权

USACO 专栏收录该内容

17 篇文章 0 订阅

订阅专栏

不是特别好想的一道题，首先领dp(i,j)表示第i位位1且前i位有j个1的数即个数，很明显j<=i&&j<=L，这是为了之后更新方便而预处理的，别忘了dp(0,0)=1，之后领dp(i,i+1)（这些单元一定用不到的）为第i行j<=i&&j<=L的和，那么，先判断最高位的1应该放在哪一位上，直接I-dp[i][i+1],若这个值小于等于0，则确定最高位在第i位上，否则的话，I-=dp[i][i+1]。之后问题便转化为前N位中至多有L-1个1的数字钟找第I个（I此时已更新）。值得注意的是在计算dp的时候，一定要防止溢出。

/*
ID: jinusac1
PROG: kimbits
LANG: C++
*/
#include <iostream>
#include <cstdio>
using namespace std;
bool r[40]={false};
long long dp[32][33]={0};
long long jie(long long num)
{
	if(num==0) return 1;
	long long r=1;
	for(long long i=1;i<=num;i++) r*=i;
	return r;
}
long long gcd(long long a, long long b)
{
	long long r=1;
	while(r){
		r=a%b;
		a=b;b=r;
	}
	return a;
}
long long cc(long long n,long long m)
{
	long long fenzi=1,fenmu=1;
	if(m>(n<<1)){
		//fenmu=jie(n);i
		//for(long long i=m-n+1;i<=m;i++) fenzi*=i;j
		long long i=1,j=m-n+1;
		while(i<=n||j<=m){
			if(i<=n) fenmu*=i++;
			if(j<=m) fenzi*=j++;
			long long t=gcd(fenzi,fenmu);
			fenzi/=t;fenmu/=t;
		}
	}
	else{
		//fenmu=jie(m-n);i
		//for(long long i=n+1;i<=m;i++) fenzi*=i;j
		long long i=1,j=n+1;
		while(i<=(m-n)||j<=m){
			if(i<=(m-n)) fenmu*=i++;
			if(j<=m) fenzi*=j++;
			long long t=gcd(fenzi,fenmu);
			fenzi/=t;fenmu/=t;
		}
	}
	return fenzi/fenmu;
}
int main()
{
	freopen("kimbits.in","r",stdin);
	freopen("kimbits.out","w",stdout);
	//freopen("in.txt","r",stdin);
	long long N,L,I;
	cin>>N>>L>>I;
	for(int i=1;i<=N;i++)
		for(int j=1;j<=i&&j<=L;j++)
			dp[i][j]=cc(j-1,i-1);
	dp[0][0]=dp[0][1]=1;
	while(L){
		for(int i=1;i<=N;i++){
			dp[i][i+1]=0;
			for(int j=1;j<=i&&j<=L;j++) dp[i][i+1]+=dp[i][j];
		}
		for(int i=0;i<=N;i++){
			if((I-dp[i][i+1])<=0){
				L--;
				r[i]=true;
				break;
			}
			else I-=dp[i][i+1];
		}
	}
	for(int i=N;i>=1;i--) cout<<r[i];
	cout<<endl;
	return 0;
}