Sicily 1327 Pinary (SOJ 1327) 【dp 动态规划】
2013-05-22 15:36
477 查看
原题地址:点击打开链接
这一题能不能算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的位置,那就知道了最终的结果。
代码非常短,如下:
******<转载说明>******
转载注明:诚实的偷包贼
原文地址:/article/1805765.html
******<转载说明/>******
这一题能不能算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(); } }
******<转载说明>******
转载注明:诚实的偷包贼
原文地址:/article/1805765.html
******<转载说明/>******
相关文章推荐
- Sicily 1346 金明的预算方案 (SOJ 1346) 【DP 动态规划-背包问题】
- Sicily 1033 City Road (SOJ 1033)【dp动态规划】
- sicily 1327. Pinary
- soj 1327. Pinary
- Sicily 1091 Maximum Sum(SOJ 1091)【dp动态规划】
- sicily 1327. Pinary
- Sicily 1019 Apple Tree (SOJ 1019)【树形dp 树形动态规划】
- POJ 1163 DP 动态规划
- 【LeetCode】132. Palindrome Partitioning II 基于动态规划DP、C++、Java的分析及解法
- DP动态规划
- USACO 1.1 Broken Necklace (DP动态规划)
- poj1661 Help Jimmy-dp动态规划
- zoj 1163 The Staircases 动态规划(dp)
- [zz] POJ 动态规划DP题目列表
- 最长递增子序列-动态规划dp-51node
- 01背包问题(动态规划DP)
- soj 1700 ping_简单dp
- LeetCode 139. Word Break 动态规划DP Python解法
- 【HDU6146】【2017百度之星复赛-C】Pokémon GO(动态规划)(dp)