UVA - 12105 Bigger is Better（DP）

传送门

---

首先，此类问题——模某个数的动态规划，一般都是从其他余数转移而来

然后本题用到了一个比较重要的大数求余的性质，也就是说大数可以一位一位得操作从而达到取余的目的，可以参考这里

方法一

状态转移

设$d(i,j)$表示用$i$根火柴能拼出的除以$m$余数为$j$的最大数，使用字符串string表示大数。因为状态的转移是互相独立的，那么考虑刷表法，设$d(i,j)$拼上当前个位数$k$后为$cur$，由当前状态更新未知状态：

$d(i+c(k),(j\ast 10+k)\%m)=max(d(i+c(k),(j\ast 10+k)\%m),cur);$

需要注意的是，这里不能直接用max，我一开始也想错了。因为字符串是按照字典序比较的，也就是说"111"是没有"9"大的，因此需要写一个字符串大小比较函数：

int cmp(string &s1,string &s2){
    if(s1.size()>s2.size()) return 1;
    if(s1.size()<s2.size()) return -1;
    if(s1==s2) return 0;
    return s1>s2?1:-1;
}

边界处理和最终答案

因为字符串是从右向左一位一位填充的，也就是说，一开始个位的$0-9$是所有其他状态转移的初态，那么我们就顺序枚举，首先设置$d(0,0)=0$，下面的拼接时需要注意，因为不能有前导$0$，那么如果第一位为$0$，就重新从个位开始

最后求答案时，就是$max(d(i,0)),0\le i\le n$

当然一定会注意到，最后求得的答案只有$0$而无法判断什么时候输出$-1$。实际上之所以初始化为$0$是因为$0$是一个特殊的状态——$0$对任何数取模都为$0$，也就是说只要有其他的任何状态，0就不会是最大的，那么我们添加一个特判，当$ans=0$时，当且仅当$n\ge 6$时输出$0$，否则为$-1$

时间复杂度

$O(10\ast n\ast m)$，STL果然不让我失望——1.8s!

#include <set>
#include <map>
#include <stack>
#include <queue>
#include <math.h>
#include <cstdio>
#include <string>
#include <bitset>
#include <cstring>
#include <sstream>
#include <iostream>
#include <algorithm>
#include <unordered_map>
using namespace std;
#define fi first
#define se second
#define pb push_back
#define ins insert
#define lowbit(x) (x&(-x))
#define mkp(x,y) make_pair(x,y)
#define mem(a,x) memset(a,x,sizeof a);
typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
typedef pair<int,int> P;
const double eps=1e-8;
const double pi=acos(-1.0);
const int inf=0x3f3f3f3f;
const ll INF=1e18;
const int Mod=1e9+7;
const int maxn=2e5+10;

int num[]={6, 2, 5, 5, 4, 5, 6, 3, 7, 6};
string d[110][3005];
int n,m;

int cmp(string &s1,string &s2){
    if(s1.size()>s2.size()) return 1;
    if(s1.size()<s2.size()) return -1;
    if(s1==s2) return 0;
    return s1>s2?1:-1;
}

int main(){
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    ios::sync_with_stdio(0),cin.tie(0),cout.tie(0);
    int kase=0;
    while(cin>>n && n){
        cin>>m;
        mem(d,0);
        d[0][0]="0";
        for(int i=0;i<=n;i++)
            for(int j=0;j<=m;j++) if(d[i][j].size())
                for(int k=0;k<10;k++){
                    char ch=k+'0';
                    string res=d[i][j];
                    if(res[0]=='0') res=ch;
                    else res+=ch;
                    if(cmp(d[i+num[k]][(j*10+k)%m],res)<0)
                        d[i+num[k]][(j*10+k)%m]=res;
                }
        string ans="0";
        for(int i=0;i<=n;i++)
            if(cmp(d[i][0],ans)>0)
                ans=d[i][0];

        cout<<"Case "<<++kase<<": ";
        if(n>=6 && ans[0]=='0'){
            cout<<"0"<<endl;
        }else if(ans[0]=='0') cout<<"-1"<<endl;
        else cout<<ans<<endl;
    }
    return 0;
}

方法二

实际上这个是比较难想和难以实现的，参考了LRJ老师的代码后不得不敬佩一番

理解这个方法我参考了这个博客，讲得非常详细了！建议一看

#include<cstdio>
#include<cstring>

using namespace std;

const int maxn = 100 + 5;
const int maxm = 3000 + 5;

int n, m, dp[maxn][maxm], p[maxn][maxm];

int needs[] = { 6, 2, 5, 5, 4, 5, 6, 3, 7, 6 };

int main(){
    int kase = 0;
    while(scanf("%d%d", &n, &m) == 2) {
        printf("Case %d: ", ++kase);

        for(int i = 0; i <= n; i++)
            for(int j = 0; j < m; j++){
                int& ans = dp[i][j];
                ans = p[i][j] = -1;
                if (j == 0) ans = 0;
                for(int d = 9; d >= 0; d--)
                    if (i >= needs[d]){
                        int t = dp[i - needs[d]][(j * 10 + d) % m];
                        if (t >= 0 && t + 1 > ans){
                            ans = t + 1;
                            p[i][j] = d;
                        }
                }
            }

        if (p[n][0] < 0) printf("-1");
        else {
            int i = n, j = 0;
            for(int d = p[i][j]; d >= 0; d = p[i][j]){
                printf("%d", d);
                i -= needs[d];
                j = (j * 10 + d) % m;
            }
        }
        printf("\n");
    }
    return 0;
}