UVA12105 Bigger is Better 详解

UVA12105 Bigger is Better

题目描述：紫书P293

用不超过 $n（1\le n\le 100）$ 根火柴摆出一个尽量大的、能被 $m（1\le m\le 3000）$ 整除的数。

解题思路

方法一：

来自紫书，设$dp(i,j)$表示用$i$根火柴能拼出除以$m$余数为$j$的最大的数,枚举在最右边添加数字$k$,用$dp(i,j)\ast 10+k$去更新$dp(i+c(k),(j\ast 10+k)\%m)$,其中$c(k)$表示数字$k$需要的火柴数。但是由于最多能有50位，需要高精度。

方法二：

也来自紫书，$dp(i,j)$表示拼出一个“除以$m$余数为$j$的$i$位数”至少需要多少根火柴，状态转移和上面差不多，我们也枚举最后一位$k$，用$dp(i,j)+c[k]$去更新$dp(i,(j\ast 10+k)$即可。
这时我们知道了答案是几位数，那么接下来从高位开始依次输出，因为要让这个数尽量大，所以对于每一位我们都按9，8，7…降序依次判断是否可行，例如$m=7$，并已经确定了最大位数为3位数。首先试着让最高位为9，如果能摆出形如9ab的整数，那它一定数最大的，又因为900除以7的余数为4，所以后两位（ab）除以7的余数应该为3，如果$d(2,3)+c(9)\le n$即可。
如何计算出x00…除以$m$的余数需要预处理，还要注意注意最高位不能为0。

方法三：

设$dp(i,j)$表示用不超过$i$根火柴摆成“还需要剩余部分除以$m$余数为$j$”的最大位数（当前余数是m-j）。

状态转移就是$dp(i,j)=max(dp[i-c[k]][(j\ast 10+k)\%m])，(0\le k\le 9)$（自己想下其实还是很好理解的和方法二的差不多）。为了方便输出，再定义一个$p(i,j)$表示$dp(i,j)$状态时的最大位的数是多少，边界条件$dp(j,0)=0$,如果我当前的数已经符合条件，那么我可以不再加，也就是余数为0时，长度为0，因为$dp(i,j)$中的$i$表示的是不超过$i$，表示即使我不在再添加数字了也可以有多出来的火柴，这时$dp(j,0)$的$j$表示的就是哪些多出来的火柴。（如果$dp(i,j)$的$i$表示的是刚好用$i$根火柴，那么只需要$dp(0,0)=0$）。

让$dp(i,j)$表示剩余$i$根（已经用了$n-i$根）余数为$j$的最大位数也是一样的，不过是从高位开始枚举，有兴趣的可以试试。

代码快

方法一大数没去做，方法二写的有点乱，没整理，这里只提供方法三的

#include <bits/stdc++.h>
using namespace std;
int dp[105][3005];
int p[105][3005];
int c[]={6,2,5,5,4,5,6,3,7,6};
int main(){
	int n,m,cas=0;
	while(cin>>n&&n)
	{
		cin>>m;
		memset(dp,-1,sizeof(dp));
		for(int i=0;i<=n;i++)
			dp[i][0]=0;
		int ans=-1,x,y;
		for(int i=1;i<=n;i++) 
			for(int j=0;j<m;j++)
				for(int k=9;k>=0;k--)
					if(i>=c[k])
					{
						if(dp[i-c[k]][(j*10+k)%m]>=0&&dp[i][j]<(dp[i-c[k]][(j*10+k)%m]+1))
						{
							dp[i][j]=dp[i-c[k]][(j*10+k)%m]+1;
							p[i][j]=k;
						}
					}

		ans=dp[n][0];
		printf("Case %d: ",++cas);
		if(ans<=0){
			printf("-1\n");continue;
		}
		x=n;y=0;
		while(ans--)
		{
			int e=p[x][y];
			printf("%d",e);
			x-=c[e];
			y=(y*10+e)%m;
		}
		printf("\n");
	}
	return 0;
}