USACO-Section3.2 Stringsobits【递归】

题目描述：

考虑排好序的N(N<=31)位二进制数。
你会发现，这很有趣。因为他们是排列好的，而且包含所有长度为N且这个二进制数中1的位数的个数小于等于L(L<=N)的数。
你的任务是输出第i（1<=i<=长度为N的二进制数的个数）小的（注：题目这里表述不清，实际是，从最小的往大的数，数到第i个符合条件的，这个意思），长度为N，且1的位数的个数小于等于L的那个二进制数。
（例：100101中，N=6，含有位数为1的个数为3）。

INPUT FORMAT:

共一行，用空格分开的三个整数N，L，i。

OUTPUT FORMAT:

共一行，输出满足条件的第i小的二进制数。.

---

SAMPLE INPUT

5 3 19

---

SAMPLE OUTPUT

10011

---

解题思路：

这道题是关于01串的，由于上限为2^32，加上判断数字0、1个数的时间，暴力枚举肯定是不可取的。那就需要研究其他的方案。这道题的意思是从N个位置中最多选出L个，并且按顺序排列，输出第i个，所以我们可以通过算出C(0,n)到C(m,n)个数来判断第n位是否用到，具体解法见代码，需要稍微思考一下。

#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#include<math.h>
#include<algorithm>
using namespace std;
FILE *fout;
int N,L,s[32];
long long int P;
int ans=0;
long long int C(int a,int b){//计算C(a,b)
long long int q=1;
    for(int i=0;i<a;i++){
        q=q*(b-i);
        q=q/(i+1);
    }
    return q;
}
long long int fun(int a,int b){//计算C(0,b)到C(a,b)的和 
    long long int temp=0;
    for(int i=0;i<=a;i++){
        temp=temp+C(i,b);
    }
    return temp;
}
void solve(int m,int n,long long int t){
    if(t==0)return ;//如果t为0返回 
    if(m==0||n==0){//如果m，n没了返回 
        return ;
    }
    long long int temp;
    if(m>n)temp=fun(n,n);   
    else temp=fun(m,n); 
    if(temp<=t){//计算fun(m,n)是否大于t 
        if(temp==t){//如果相等说明取从n-1开始的前m位为1，在这里就一位一位的取了，写起来方便 
            ans+=(long long int)pow(2,n-1); 
            solve(m-1,n-1,t-fun(m,n-1));    
        }
        else if(fun(m,n+1)>t){//如果fun(m,n)小于t且fun(m,n+1)大于t说明第n+1位必取 
            ans+=(long long int)pow(2,n);
            solve(m-1,n,t-temp);
        }
        else{//如果fun(m,n+1)也没t大就继续 
            solve(m,n+1,t);
        }
    }
    else{
        solve(m,n-1,t);
    }
}
void print(){//输出结果 
    int ccount=0;
    while(ans){
        if(ans&1)s[ccount++]=1;
        else s[ccount++]=0;
        ans>>=1;
    }
    for(int i=N-1;i>=0;i--)
    fprintf(fout,"%d",s[i]);
    fprintf(fout,"\n"); 
} 
int main(){
    FILE *fin  = fopen ("kimbits.in", "r");
    fout = fopen ("kimbits.out", "w");
    fscanf(fin,"%d%d%lld",&N,&L,&P);
    solve(L,L,P);
    print();
    exit(0);
}