HDU 1502 Regular Words DP+高精度

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1502

题目大意:找出总的满足条件的字符串数,num(a)=num(b)=num(c)且任何前缀均满足num(a)>=num(b)>=num(c)

解题思路:用dp[i][j][k]表示a取i个,b取j个,c取k个的状态下最多有多少种满足条件的情况,容易推得状态转移方程如下:

dp[i][j][k]=dp[i-1][j][k](i>j时)+dp[i][j-1][k](j>k时)+dp[i][j][k-1](k>0时)

根据该状态转移方程即可输出最后的结果dp[n][n][n]

因为本题的数据结果比较大,所以还需要用到高精度,对不会用java的人就只能手敲一个大数相加了。。。

 #include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define scan(a) scanf("%d",&a)
using namespace std;
#define N 61
char dp[N][N][N][],a[];
void init()
{
strcpy(dp[][][],"1\0");
strcpy(dp[][][],"1\0");
strcpy(dp[][][],"0\0");
strcpy(dp[][][],"1\0");
} void add(char *p)
{
int x,l,i,jin=;
l=strlen(a);
int now=;
for(i=;p[i]!='\0';i++)
//从加数开始一位一位访问
{
if(i>=l)
x=p[i]-''+jin;//i超过了a的长度
else
x=p[i]-''+a[i]-''+jin;
if(x>)
{
jin=x/;
x%=;
}
else
jin=;
a[now++]=x+'';
}
while(jin)
{
if(l<=now)
{
a[now++]=(jin%)+'';
jin/=;
}
else
{
x=a[now]-'';
x+=jin;
if(x>)
{
jin=x/;
x%=;
}
else
jin/=;
a[now++]=x+'';
}
}
if(now>l)
a[now]='\0';
} void put(int x)
{
int l=strlen(dp[x][x][x]);
for(int i=l-;i>=;i--)
cout<<dp[x][x][x][i];
cout<<endl<<endl;
} int main()
{
int i,j,k;
int n;
init();
for(i=;i<=;i++)
{
for(j=;j<=i;j++)
{
for(k=;k<=j;k++)
{
a[]='';
a[]='\0';
// cout<<i<<' '<<j<<' '<<k<<endl;
// cout<<dp[i-1][j][k]<<' '<<dp[i][j-1][k]<<' '<<dp[i][j][k-1]<<endl;
if(i>j&&j>=k)
add(dp[i-][j][k]);
if(j>k)
add(dp[i][j-][k]);
if(k>)
add(dp[i][j][k-]);
strcpy(dp[i][j][k],a);
}
}
}
while(~scanf("%d",&n))
{
put(n);
}
return ;
}
上一篇:SRM 627 D1L2GraphInversionsDFS查找指定长度的所有路径 Binary indexed tree (BIT)


下一篇:Java字节流:InputStream OutputStream