BUPT Summer Journey #test1 E

E. 数的关系 2014新生暑假个人排位赛01

时间限制 5000 ms     内存限制 65536 KB    

题目描述

用关系“＜”和“=”将3个数A、B和C依序排列时有13种不同的序关系：
A＝B＝C，A＝B＜C，A＜B＝C，A＜B＜C，A＜C＜B，A＝C＜B，B＜A＝C， 

 

B＜A＜C，B＜C＜A，B＝C＜A，C＜A＝B，C＜A＜B，C＜B＜A 

 

。

现在输入数字的个数，要求你给出上述关系的数目。

数的个数不大于100
 

输入格式

多组数据，EOF结束

每行一个输入

输出格式

对于每个输入，输出一行，即对应答案

输入样例

3

输出样例

13

 

思路：一开始。我想枚举=号与《号的个数，然后后来发现A=B<C=D<E与A=B=C<D<E时会出现错误。后来思路就一直卡死在那儿。后来经过May的指点，其实相等的数地位是相同的，不应该去考虑等号，这样会衍生很多种情况。递推的话比较难。本来递推就是在递推之间略去了原来的计算因素从而简化计算量，故所以我们可以枚举一共需要多少个值m.把n个数放入m个值里面。此时递推关系就引刃而解了。设dp[n][m]为把n个数取m个值的方法。首先，它可以相当于n-1个数放入m个值里面，则第n个数在m个值里面任取一个值有m*dp[n-1][m]，其次，还可以相当于n-1个数放入m-1个值后，剩下的一个数可以取m也可以取1到m-1所以它的种数为m*dp[n-1][m-1]。别急。看有没有遗漏的。首先，N-2个数的状态是N-1个数的状态的子状态，包含在N-1的情况里面。其余，新增加的那个元素最多取一个新加入的值，故m也只有m和m-1两种前驱状态。不重不漏。递推方程为：

dp[n][m]=m*(dp[n-1][m-1]+dp[n-1][m]).

初始值为：dp[1][1]=1;dp[i][1]=0(i!=1)(可以自己推几个数确定未定义出的值比如2dp[2][1],dp[2][2]的值)

代码：

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
using namespace std;
int n;
int dp[101][101][1000],ans[101][1000];
void print(int p[])
{
    printf("%d",p[p[0]]);
    for(int i=p[0]-1;i>=1;i--)
    {
        printf("%d",p[i]/1000%10);
        printf("%d",p[i]/100%10);
        printf("%d",p[i]/10%10);
        printf("%d",p[i]/1%10);
    }
    printf("\n");
}
void HighAdd(int c[],int a[],int b[])
{
    c[0]=a[0];
    if(a[0]<b[0])c[0]=b[0];
    for(int i=1;i<=c[0];i++)c[i]=a[i]+b[i];
    for(int i=1;i<=c[0];i++){c[i+1]+=c[i]/10000;c[i]%=10000;}
    if(c[c[0]+1]>0)c[0]++;
}
void SingleMultiply(int a[],int b)
{
    if(b==0){for(int i=0;i<2000;i++)a[i]=0;return;}
    for(int i=1;i<=a[0];i++)a[i]*=b;
    for(int i=1;i<=a[0];i++){a[i+1]+=a[i]/10000;a[i]%=10000;}
    while(a[a[0]+1]>0){a[0]++;a[a[0]+1]=a[a[0]]/10000;a[a[0]]%=10000;}
}
int main()
{
    memset(dp,0,sizeof(dp));
    memset(ans,0,sizeof(ans));
    dp[1][1][0]=1;dp[1][1][1]=1;
    for(int i=2;i<=100;i++){dp[1][i][0]=1;dp[1][i][1]=0;}
    //for(int i=1;i<=100;i++)print(dp[1][i]);
    for(int i=2;i<=100;i++)
        for(int j=100;j>=1;j--)
        {
            HighAdd(dp[i][j],dp[i-1][j],dp[i-1][j-1]);
            //printf("dp[%d][%d]=",i,j);print(dp[i][j]);
            SingleMultiply(dp[i][j],j);
        }
    //for(int i=1;i<=2;i++)print(dp[2][i]);
    for(int i=1;i<=100;i++)
        for(int j=1;j<=i;j++)HighAdd(ans[i],ans[i],dp[i][j]);
    while(scanf("%d",&n)==1)
    {
        print(ans[n]);
    }
    return 0;
}