Sicily 1327 Pinary （SOJ 1327） 【dp 动态规划】

原题地址：点击打开链接

这一题能不能算DP呢，其实只是一个简单的迭代，不过也可以说是DP。

初看这条题我是没有什么想法，思路比较乱，但是想到了状态转移之后就豁然开朗。

——————正文——————

这题的时间卡的比较紧，数据规模最大有9*10^7，而1s大概就允许复杂度为10^8左右，所以比较高效的算法才能够AC

假设f[i]表示当结果的位数为i时，i为所能表示的所有数字的数量，即便当i为3时，f[i]=2，因为有100和101

这样就比较容易得到递推公式：f[i]=f[i-2]+f[i-1]+...+f[2]+f[1]+1

这是因为对于i位数，其总共的情况数就等于——（左数起）第2位到第j-1位都为0，第j位为1的时候——所有的情况数之和。

同时我们发现，计算f[i-1]时，有f[i-1]=f[i-3]+f[i-4]+...+f[1]+1，计算f[i]时，有f[i]=f[i-2]+f[i-3]+...+f[1]+1，等式的右边被重复计算，于是我们使用一个变量sum[i]来保存f[1]+...+f[i]的值

这样就有以下表达式：

f[1]=f[2]=1;

sum[1]=1; sum[2]=2;

f[i]=sum[i-2]+1;

sum[i]=sum[i-1]+f[i] => sum[i]=sum[i-1]+sum[i-2]+1;

因此发现根本不需要计算f[i]，因此最终的代码里面直接使用sum就可以了。

计算得到sum之后，从sum[i]我们就知道了有i位及以下位数的数字其总共可以表达的数有多少个。

根据这个，我们对于每一个输入tmp，都不断查找刚好超过当前tmp的sum[i]，这样就表示最终的结果在i位必定是1，然后tmp减去sum[i-1]+1（思考一下为什么呢？），再重复本句的步骤即可，知道了所有1的位置，那就知道了最终的结果。

代码非常短，如下：

//f[i]=f[i-2]+f[i-3]...+f[1]+1#include<iostream>#include<stdio.h>#include<memory.h>#include<cstring>//long long f[63];long long sum[45];long long tmp=1;int b[1000];int n;void compute(){memset(b,0,sizeof(b));int i,maxi=0,bound=40;while(tmp!=0){  i=1;  for(;i<=bound && sum[i]<tmp;i++)  {}  maxi=(maxi>i)?maxi:i;  b[i]=1;  tmp-=sum[i-1]+1;  bound=i-2;}for(i=maxi;i>=1;i--)  printf("%d",b[i]);printf("\n");}int main(){//freopen("in","r",stdin);//freopen("out","w",stdout);//f[0]=0; f[1]=f[2]=1;sum[1]=1; sum[2]=2;for(int i=3;i<=40;++i)sum[i]=sum[i-1]+sum[i-2]+1;scanf("%d",&n);while(n--){scanf("%lld",&tmp);compute();}}

******<转载说明>******
转载注明：诚实的偷包贼
原文地址：http://blog.csdn.net/fanfank/article/details/8960563
******<转载说明/>******