USACO 2012 February Silver

最新推荐文章于 2022-09-16 10:47:41 发布

Caristra

最新推荐文章于 2022-09-16 10:47:41 发布

阅读量229

点赞数

分类专栏： USACO 文章标签：组合数高精度

本文链接：https://blog.csdn.net/Caristra/article/details/83855439

版权

USACO 专栏收录该内容

3 篇文章 0 订阅

订阅专栏

Description:
求第 $n$ 个刚好有 $k$ 个 $1$ 的二进制正整数。
$n\le10^{30},2\le k\le20$ ，答案长度不超过 $1000$

Solution:

比较常规的组合数，对于每一位 $i$ ，我们判断 $C_{n-1}^{i-1}$ 是否大于当前的值即可。
再套一个高精度。

Code:

#include<bits/stdc++.h>
using namespace std;
#define REP(i,f,t) for(int i=(f),i##_end_=(t);i<=i##_end_;++i)
#define SREP(i,f,t) for(int i=(f),i##_end_=(t);i<i##_end_;++i)
#define DREP(i,f,t) for(int i=(f),i##_end_=(t);i>=i##_end_;--i)
#define ll long long
template<class T>inline bool chkmin(T&x,T y){return x>y?x=y,1:0;}
template<class T>inline bool chkmax(T&x,T y){return x<y?x=y,1:0;}
template<class T>inline void rd(T&x){
	x=0;char c;
	while((c=getchar())<48);
	do x=(x<<1)+(x<<3)+(c^48);
	while((c=getchar())>47);
}

const int N=302;

int K;

#define L 255
#define P 10000
 
struct Big{
    int num[L],len;
    Big(){
        memset(num,0,sizeof(num));
        len=1;
    }
    void Rd(){
        char A[N];
        scanf("%s",A);
        int SL=strlen(A);
        len=0;
        for(int i=SL-1,res;res=0,i>=0;num[len++]=res,i-=4){
            if(i>=3)for(int j=i-3;j<=i;j++)res=(res<<1)+(res<<3)+(A[j]^48);
            else for(int j=0;j<=i;j++)res=(res<<1)+(res<<3)+(A[j]^48);
        }
    }
    void Pr(){
      	if(!num[len-1]){putchar("0");return;}
        printf("%d",num[len-1]);
        DREP(i,len-2,0)printf("%04d",num[i]);       
    }
    Big operator +(const Big &a)const{
        Big b;
        b.len=max(len,a.len);
        REP(i,0,b.len-1){
            int &B=b.num[i];
            B+=num[i]+a.num[i];
            if(B>=P)B-=P,b.num[i+1]++;
        } 
        if(b.num[b.len])b.len++;
        return b;
    }
    Big operator+(int a){
        Big b;
        b.len=len;
        REP(i,0,len-1){
            int &B=b.num[i];
            B+=num[i]+a%P;
            if(B>=P)b.num[i+1]++,B-=P;
            a/=P;
        }
        if(b.num[b.len])b.len++;
        return b;
    }
    bool operator<=(const Big &a)const{
        if(len!=a.len)return len<a.len;
        DREP(i,len-1,0) if(num[i]!=a.num[i]) return num[i]<a.num[i];
        return 1;
    }
}n,now,tmp,C[N][22];
int pos[30];

int get(int k){
	for(int i=k;;i++){
		tmp=C[i-1][k-1];
		if(n<=now+tmp) return i;
		else now=now+tmp;
	}
}

int main(){
//	freopen("num.in","r",stdin);
//	freopen("num.out","w",stdout);
	n.Rd();
	scanf("%d",&K);
	
	C[0][0]=C[0][0]+1;
	REP(i,1,300){
		C[i][0]=C[i][0]+1;
		REP(j,1,min(20,i)) C[i][j]=C[i-1][j-1]+C[i-1][j];
	}
	
	SREP(i,0,K) pos[i]=get(K-i); 
    pos[K]=0;
    SREP(i,0,K){
		putchar('1');
        while(pos[i]-1>pos[i+1])pos[i]--,putchar('0');
    }
    putchar('\n');
	
	return 0;
}