luoguP1018 乘积最大 dfs，dp，高精

链接：https://www.luogu.com.cn/problem/P1018

题意，给定一串数字，长度最大为40，在数字中插入k个乘号，问怎样插入使得乘积最大。

解法：最开始用dfs写，还是很不错的，过了8个样例，但是__int128也不支持超过40位的数字，所以最后两个样例总过了。发现题解也没有dfs的，就采取了他们的dp做法。

dp:dp[k][i]表示前i个数字插入k个乘号的最大值，每次枚举最后一个乘号的位置，然后dp[k][i]=max(dp[k][i],dp[k-1][j]*xx(j+1,i)).xx(j+1,i)表示第k+1串字符串组成的数字emmmm.懂的都懂。

那么只需要在每一次解dp[k][i]的时候已经解出了dp[k-1][i]（这是一个学习dp很优的思维）。首先初始化dp[0][i]。然后循环初始化，注意高精的写法，也是值得学习的。

代码：

#include <bits/stdc++.h>

#define ll int
#define pi acos(-1)
#define pb push_back
#define mst(a, i) memset(a, i, sizeof(a))
#define pll pair<ll, ll>
#define fi first
#define se second
#define dbg(x) cout << #x << "===" << x << endl
#define dbgg(i, x) cout << #x << "[" << i << "]===" << x[i] << endl
using namespace std;
template<class T>void read(T &x){T res=0,f=1;char c=getchar();while(!isdigit(c)){if(c=='-')f=-1;c=getchar();}while(isdigit(c)){res=(res<<3)+(res<<1)+c-'0';c=getchar();}x=res*f;}
void print(ll x){if(x<0){putchar('-');x=-x;}if(x>9)print(x/10);putchar(x%10+'0');}
const ll maxn = 50 + 10;
const ll mod = 1e9 + 7;


struct node{
    ll len,ans[maxn];
}dp[maxn][maxn];//?str[k][i]表示第数组第i个数前放k个乘号得到的最大值。
ll n,m,a[maxn];
char s[maxn];
node cmp(node x,node y){
    ll lx=x.len,ly=y.len;
    if(lx>ly) return x;
    if(lx<ly) return y;
    for(ll i=lx;i>=1;i--){
        if(x.ans[i]>y.ans[i]) return x;
        if(x.ans[i]<y.ans[i]) return y;
    }
    return x;
}
node mul(node x,ll l,ll r){
    node y,z;
    mst(y.ans,0);
    mst(z.ans,0);
    for(ll i=r;i>=l;i--) y.ans[r-i+1]=a[i];
    ll lx=x.len,ly=r-l+1;
    for(ll i=1;i<=lx;i++){
        for(ll j=1;j<=ly;j++){
            z.ans[i+j-1]+=x.ans[i]*y.ans[j];
        }
    }
    ll L=lx+ly-1;
    for(ll i=1;i<=L;i++){
        z.ans[i+1]+=z.ans[i]/10;
        z.ans[i]%=10;//?这里不用判断z.ans[i]>10，也不能判断！！！
    }
    if(z.ans[L+1]) ++L;
    z.len=L;
    return z;
}
ll gcd(ll a,ll b){return (b==0)?a:gcd(b,a%b);}
ll qpow(ll a,ll p,ll mod){ll ans=1;a=a%mod;while(p){if(p&1)ans=(ans*a)%mod;p>>=1;a=(a*a)%mod;}return ans;}
int main() {
    ll _s = 1;
    //read(_s);
    for (ll _=1;_<=_s;_++) {
        read(n),read(m);
        scanf("%s",s+1);
        for(ll i=1;i<=n;i++) a[i]=s[i]-'0';
        for(ll i=1;i<=n;i++){
            //dp[0][i].len=0;
            for(ll j=i;j>=1;j--){
                dp[0][i].ans[++dp[0][i].len]=a[j];
            }
        }


        for(ll i=2;i<=n;i++){//?第k个乘号在i-th之前
            for(ll k=1;k<=min(m,i-1);k++){
                for(ll j=k;j<i;j++){//?枚举乘号的最后一个位置
                    dp[k][i]=cmp(dp[k][i],mul(dp[k-1][j],j+1,i));
                }
            }
        }
        //dbg(dp[m][n].len);
        ll p=dp[m][n].len;
        for(ll i=p;i>=2;i--){
            if(dp[m][n].ans[i]){
                p=i;break;
            }
        }//?有还是无妨的
        for(ll i=p;i>=1;i--)
            printf("%d",dp[m][n].ans[i]);
       // print(ans);
    }
    return 0;
}
/*
input:::
40 3
9380408868346641446585285380108133256352
output:::
6051462042301381677936607451948047334400
*/