codevs 1157 2k进制数

1157 2k进制数

 

2006年NOIP全国联赛提高组

 时间限制: 1 s

 空间限制: 128000 KB

 题目等级 : 黄金 Gold

 

题目描述 Description

设r是个2k进制数，并满足以下条件：

（1）r至少是个2位的2k进制数。

（2）作为2k进制数，除最后一位外，r的每一位严格小于它右边相邻的那一位。

（3）将r转换为2进制数q后，则q的总位数不超过w。

在这里，正整数k（1≤k≤9）和w（k<W< span>≤30000）是事先给定的。

问：满足上述条件的不同的r共有多少个？

我们再从另一角度作些解释：设S是长度为w 的01字符串（即字符串S由w个“0”或“1”组成），S对应于上述条件（3）中的q。将S从右起划分为若干个长度为k 的段，每段对应一位2k进制的数，如果S至少可分成2段，则S所对应的二进制数又可以转换为上述的2k进制数r。

例：设k=3，w=7。则r是个八进制数（23=8）。由于w=7，长度为7的01字符串按3位一段分，可分为3段（即1，3，3，左边第一段只有一个二进制位），则满足条件的八进制数有：

2位数：高位为1：6个（即12，13，14，15，16，17），高位为2：5个，…，高位为6：1个（即67）。共6+5+…+1=21个。

3位数：高位只能是1，第2位为2：5个（即123，124，125，126，127），第2位为3：4个，…，第2位为6：1个（即167）。共5+4+…+1=15个。

所以，满足要求的r共有36个。

输入描述 Input Description

只有1行，为两个正整数，用一个空格隔开：

k W

输出描述 Output Description

共1行，是一个正整数，为所求的计算结果，即满足条件的不同的r的个数（用十进制数表示），要求最高位不得为0，各数字之间不得插入数字以外的其他字符（例如空格、换行符、逗号等）。

（提示：作为结果的正整数可能很大，但不会超过200位）

样例输入 Sample Input

3 7

样例输出 Sample Output

36

数据范围及提示 Data Size & Hint

 

/*
F[i][j]表示以i开头的长度为j的有几个
F[i][j]=f[i+1][j-1]+f[i+2][j-1]+...+f[maxx][j-1]
但这样做会超时
观察到：
f[i+1][j]=f[i+2][j-1]+...+f[maxx][j-1]
F[i][j]=f[i+1][j-1]+f[i+1][j]
*/

/*90分超时代码*/
#include<iostream>
#include<cstdio>
#include<cstring>
#define LL unsigned long long
#define maxn 30010
using namespace std;
int p,w,r,l,len=1;
int a[maxn],f[15];
struct node
{
    int l,a[210]; 
}g[520][610],ans;
void add(node &x,node &y)
{
    int l1=x.l,l2=y.l,l3=1,i,j,k;
    LL c[210]={0};
    while(l3<=l1||l3<=l2)
    {
        c[l3]+=(x.a[l3]+y.a[l3]);
        c[l3+1]+=c[l3]/1000000000;
        c[l3]%=1000000000;
        l3++;
    }
    while(l3>1&&c[l3]==0)l3--;
    x.l=l3;
    for(i=1;i<=l3;i++)
    x.a[i]=c[i];
}
void prepare()
{
    LL i,j,k;
    f[0]=1;
    for(LL i=1;i<=9;i++)
    f[i]=f[i-1]<<1;
    for(i=1;i<f[p];i++)
    {
        int q=f[p]-i-1;
        while(q>0)
        {
            g[i][2].l++;
            g[i][2].a[g[i][2].l]=q%1000000000;
            q=q/1000000000;
        }
    }
    for(j=3;j<=l+1;j++)
      for(i=1;i<f[p];i++)
        for(k=i+1;k<f[p];k++)
        add(g[i][j],g[k][j-1]);
}
int main()
{
    LL i,j,k;
    scanf("%d%d",&p,&w);
    r=w%p;
    l=w/p;
    prepare();
    for(i=1;i<=l;i++)
    a[i]=f[p];
    a[l+1]=f[r];
    for(j=2;j<=l+1;j++)
      for(i=1;i<a[j];i++)
        add(ans,g[i][j]);
    printf("%d",ans.a[ans.l]);
    for(i=ans.l-1;i>=1;i--)
    printf("%09d",ans.a[i]);
    cout<<endl;
    return 0;
}

/*AC代码*/
#include<iostream>
#include<cstdio>
#include<cstring>
#define LL unsigned long long
#define maxn 30010
using namespace std;
int p,w,r,l,len=1;
int a[maxn],f[15];
struct node
{
    int l,a[210]; 
}g[520][610],ans;
void add(node &x,node &y)
{
    int l1=x.l,l2=y.l,l3=1,i,j,k;
    LL c[210]={0};
    while(l3<=l1||l3<=l2)
    {
        c[l3]+=(x.a[l3]+y.a[l3]);
        c[l3+1]+=c[l3]/1000000000;
        c[l3]%=1000000000;
        l3++;
    }
    while(l3>1&&c[l3]==0)l3--;
    x.l=l3;
    for(i=1;i<=l3;i++)
    x.a[i]=c[i];
}
void prepare()
{
    LL i,j,k;
    f[0]=1;
    for(LL i=1;i<=9;i++)
    f[i]=f[i-1]<<1;
    for(i=1;i<f[p];i++)
    {
        int q=f[p]-i-1;
        while(q>0)
        {
            g[i][2].l++;
            g[i][2].a[g[i][2].l]=q%1000000000;
            q=q/1000000000;
        }
    }
    for(j=3;j<=l+1;j++)
      for(i=f[p]-1;i>=1;i--)
      {
          add(g[i][j],g[i+1][j]);
          add(g[i][j],g[i+1][j-1]);
      }
}
int main()
{
    LL i,j,k;
    scanf("%d%d",&p,&w);
    r=w%p;
    l=w/p;
    prepare();
    for(i=1;i<=l;i++)
    a[i]=f[p];
    a[l+1]=f[r];
    for(j=2;j<=l+1;j++)
      for(i=1;i<a[j];i++)
        add(ans,g[i][j]);
    printf("%d",ans.a[ans.l]);
    for(i=ans.l-1;i>=1;i--)
    printf("%09d",ans.a[i]);
    cout<<endl;
    return 0;
}