HDOJ - 3555 不会按位DP的做法

本文详细解析了一道关于数字49出现频率的算法题目,通过构建数学模型和使用动态规划的思想逐步推导出解决方案。文章还介绍了如何避免重复计算,并提供了一个具体的C++实现示例。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

今天做的一套去年WHU出的多校联合题...这道题纠结了好久才挤出来....去网上搜了下...大牛们似乎都是用的高端的按位DP...我等小菜就只能做屎的想赖推了....

首先将每个位的表打出来...例如sum[3][1]代表100以下的有多少个包含49的数...sum[5][5]代表50000以下有多少包含49的数..首先要看10^k的sum会是多少...也就是sum[k][1]都是多少...其中K从0到10^19 ( 比2^64大的最小位置)....

稍微推一下...譬如知道sum[3][1]=20..也就是1000以内的是20..那么1000 ~ 2000也是20...2000~3000也是20...3000~4000也是20..以此类推..则有20*10个了..然后在4000到5000中还有49XX的没有记入..49XX的一共就是10*10一百个..但是4949会重复计算一次..所以是10*10-1 ( 往大的推会发现这里的 “ 1 " 是sum[i-2][1])...

所以10的阶乘的 sum[i][1]=sum[i-1][1]*10+10^(i-1)-sum[i-2][1]....

然后在1~4时...sum[i][j]=j*sum[i][1]...这个很好想...每一段和 0 ~ 10^i是一样的...有j段就乘以j......

但在5~9时...思想要转过来...因为在4XXXXX ~ 5XXXXX我们又有更多的包含"49"数没有记入...要加上去.则: sum[i][j]=j*sum[i][1]+10^i-sum[i-1][1];

做好了这个表后....下面的工作就轻松多了...但也要注意...

如果简单的 5983 = sum[4][5]+sum[3][9]+sum[2][8]+sum[1][3] 大部分是会对...但交上去会WA..原因是没有考虑本身这个数就含了 " 49" 这种情况..

就拿最浅显的来说 ..49 = sum[2][4] + sum[1][9]...的出来就会是0...再来..493 = sum[3][4]+sum[2][9]+sum[1][3] 得出来也会少...稍微思考一下这种情况也不难处理

如果我扫到了一个49..那么就做完9的sum累加后break...ans+=后面的数+1....如49=sum[2][4]+sum[1][9]+(0+1)..493=sum[3][4]+sum[2][9]+(3+1);这样就能保证万无一失了...

这样做比较费脑细胞...也很费时间...这道题比较正规的方法应该是按位DP...本菜刚才看了好几个文章都没来感觉...慢慢来吧..

Program:


#include<iostream> #include<stdio.h> #include<string.h> using namespace std; long long sum[20][12],_10jie[20]; int t,i,n,x; char s[100]; void Prework() { int i,j; long long k=1; memset(sum,0,sizeof(sum)); for (j=5;j<=9;j++) sum[2][j]=1; for (i=3;i<=19;i++) { sum[i][1]=sum[i-1][1]*10+k-sum[i-2][1]; k*=10; for (j=2;j<=4;j++) sum[i][j]=j*sum[i][1]; for (j=5;j<=9;j++) sum[i][j]=j*sum[i][1]+k-sum[i-1][1]; } _10jie[1]=1; for (i=2;i<=19;i++) _10jie[i]=_10jie[i-1]*10; } long long GetAnswer() { int k=1; long long ans=0,y=1; bool f=false; for (i=1;i<=n;i++) { x=s[i]-'0'; ans+=sum[n-k+1][x]; if (x==9) { if (s[i-1]=='4') { i+=1; f=true; break; } } k++; } if (f) { y=0; for (;i<=n;i++) y=y*10+s[i]-'0'; ans+=y+1; } return ans; } int main() { // freopen("B.in","r",stdin); freopen("B.out","w",stdout); Prework(); scanf("%d",&t); while (t--) { scanf("%s",s+1); n=strlen(s+1); i=1; printf("%I64d\n",GetAnswer()); } return 0; }

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值