BZOJ 1111: [POI2007]四进制的天平Wag DP+高精度
2016-11-03 22:43
225 查看
时空隧道
还是旁边的ysq童鞋提醒我才发现这个题号如此…汪汪汪的悲哀QAQ
分析:
看了看题目再看看砝码都是4的整次幂…就知道是要转化成4进制再干别的事情…
转化完之后观察发现每一位只有4种可能0、1、2、3…
0:什么都不用做
1:如果只考虑这一位的话显然放砝码到右盘更优
2:放左盘右盘都一样
3:放在左盘更优(就是在左盘方1个4^x的砝码,右盘放1个4^(x+1)的砝码)
然后要求方案数?大概只能DP了…
f[i][0]代表第i位放在右盘
f[i][1]代表第i位放在左盘
f[i][0]=min(f[i+1][0],f[i+1][1])+lala[i]
f[i][1]=min(f[i+1][0]+1,f[i+1][1]-1)+4-lala[i]
在f[i][1]转移的时候f[i+1][0]+1是因为要从上一个拿走4的次幂拿走一个放在这里….(自己理解一下吧…本人语文水平有限QAQ)
代码如下:
by >_< NeighThorn
还是旁边的ysq童鞋提醒我才发现这个题号如此…汪汪汪的悲哀QAQ
分析:
看了看题目再看看砝码都是4的整次幂…就知道是要转化成4进制再干别的事情…
转化完之后观察发现每一位只有4种可能0、1、2、3…
0:什么都不用做
1:如果只考虑这一位的话显然放砝码到右盘更优
2:放左盘右盘都一样
3:放在左盘更优(就是在左盘方1个4^x的砝码,右盘放1个4^(x+1)的砝码)
然后要求方案数?大概只能DP了…
f[i][0]代表第i位放在右盘
f[i][1]代表第i位放在左盘
f[i][0]=min(f[i+1][0],f[i+1][1])+lala[i]
f[i][1]=min(f[i+1][0]+1,f[i+1][1]-1)+4-lala[i]
在f[i][1]转移的时候f[i+1][0]+1是因为要从上一个拿走4的次幂拿走一个放在这里….(自己理解一下吧…本人语文水平有限QAQ)
代码如下:
#include<algorithm> #include<iostream> #include<cstring> #include<cstdio> //by NeighThorn using namespace std; const int maxn=2000+5,MOD=1e9;; char m[maxn]; int a[maxn],b[maxn],len,num,lala[maxn],cnt,f[maxn][2],ans[maxn][2]; inline int chu(int a[],int b[]){ int x=0; for(int i=1;i<=len;i++) b[i]=(x*10+a[i])/4,x=(x*10+a[i])%4; num=1; while(b[num]==0&&num<len) num++; for(int i=num;i<=len;i++) b[i-num+1]=b[i]; len-=num-1; return x; } //f[i][0]代表第i位放在右盘 //f[i][1]代表第i位放在左盘 //f[i][0]=min(f[i+1][0],f[i+1][1])+lala[i] //f[i][1]=min(f[i+1][0]+1,f[i+1][1]-1)+4-lala[i] signed main(void){ cnt=0; scanf("%s",m);len=strlen(m); for(int i=1;i<=len;i++) a[i]=m[i-1]-'0'; while(a[1]!=0){ lala[++cnt]=chu(a,b); for(int i=1;i<=len;i++) a[i]=b[i]; }lala[++cnt]=0;lala[0]=0; /* for(int i=1;i<=cnt;i++) cout<<lala[i]; cout<<endl; */ memset(f,0x3f3f3f3f,sizeof(f)); f[cnt][0]=0;ans[cnt][0]=1; for(int i=cnt-1;i>=0;i--){ if(f[i+1][0]==f[i+1][1]) f[i][0]=f[i+1][0]+lala[i],ans[i][0]=(ans[i+1][0]+ans[i+1][1])%MOD; else if(f[i+1][0]>f[i+1][1]) f[i][0]=f[i+1][1]+lala[i],ans[i][0]=ans[i+1][1]; else f[i][0]=f[i+1][0]+lala[i],ans[i][0]=ans[i+1][0]; if(f[i+1][0]+1==f[i+1][1]-1) f[i][1]=f[i+1][0]+1+4-lala[i],ans[i][1]=(ans[i+1][0]+ans[i+1][1])%MOD; else if(f[i+1][0]+1>f[i+1][1]-1) f[i][1]=f[i+1][1]-1+4-lala[i],ans[i][1]=ans[i+1][1]; else f[i][1]=f[i+1][0]+1+4-lala[i],ans[i][1]=ans[i+1][0]; } cout<<ans[0][0]<<endl; return 0; }
by >_< NeighThorn
相关文章推荐
- bzoj1111[POI2007]四进制的天平Wag 高精度+dp
- BZOJ 1111: [POI2007]四进制的天平Wag 进制拆分,DP求方案数
- BZOJ 1111 POI2007 四进制的天平Wag 高精度+动态规划
- [二幂拆分] hihoCoder Challenge 29 B. 快速乘法 & BZOJ 1111[POI2007]四进制的天平Wag
- BZOJ1111 : [POI2007]四进制的天平Wag
- BZOJ1111: [POI2007]四进制的天平Wag
- BZOJ 1111: [POI2007]四进制的天平Wag
- BZOJ 1111 POI2007 四进制的天平
- bzoj1111 [POI2007]四进制的天平Wag
- BZOJ1111 [POI2007]四进制的天平Wag
- bzoj1111: [POI2007]四进制的天平Wag
- bzoj 1111 - 四进制的天平
- BZOJ 1097: [POI2007]旅游景点atr( 最短路 + 状压dp )
- bzoj 1089: [SCOI2003]严格n元树 (dp+高精度)
- [数位DP 高精度 拓扑排序 bitset] BZOJ 2913 [Poi1997]XOR Gates
- [BZOJ]1089: [SCOI2003]严格n元树 DP+高精度
- BZOJ 1655: [Usaco2006 Jan] Dollar Dayz 奶牛商店 dp+高精度
- bzoj 1097 [POI2007]旅游景点atr(最短路,状压DP)
- [BZOJ]1210: [HNOI2004]邮递员 插头DP+高精度
- bzoj2817[ZJOI2012]波浪 DP+高精度