LOJ P10149 凸多边形的划分 题解

Analysis

区间dp+压位高精

dp五分钟，高精两小时

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define maxn 110
#define INF 2147483647
using namespace std;
typedef long long ll;
ll dp[maxn][maxn][maxn],a[maxn];
ll s1[maxn],s2[maxn],s3[maxn];
int n;
inline void mark(ll c[])
{
    for(int i=;i<=c[];i++)
    {
        c[i+]+=c[i]/;
        c[i]%=;
    }
    while(c[c[]+])
    {
        c[]++;
        c[c[]+]+=c[c[]]/;
        c[c[]]%=;
    }
}
inline void mul(ll a1,ll a2,ll a3,ll c[])
{
    c[]=;
    c[]=;
    for(int i=;i<=c[];i++) c[i]*=a1;
    mark(c);
    for(int i=;i<=c[];i++) c[i]*=a2;
    mark(c);
    for(int i=;i<=c[];i++) c[i]*=a3;
    mark(c);
}
inline void add(ll a[],ll b[],ll c[])
{
    c[]=max(a[],b[]);
    for(int i=;i<=c[];i++) c[i]=a[i]+b[i];
    mark(c);
}
inline bool compare(ll a[],ll b[])
{
    if(a[]>b[]) return ;
    else if(a[]<b[]) return ;
    for(int i=a[];i>=;i--)
    {
        if(a[i]>b[i]) return ;
        else if(a[i]<b[i]) return ;
    }
    return ;
}
inline void print()
{
    printf("%lld",dp[][n][dp[][n][]]);
    for(int i=dp[][n][]-;i>;i--)
    {
        printf("%lld",dp[][n][i]/);
        printf("%lld",dp[][n][i]/%);
        printf("%lld",dp[][n][i]/%);
        printf("%lld",dp[][n][i]%);
    }
}
int main()
{
    scanf("%d",&n);
    for(int i=;i<=n;i++) scanf("%lld",&a[i]);
    for(int i=;i<=n;i++)
        for(int j=;j<=n;j++)
            dp[i][j][]=;
    for(int l=;l<=n-;l++)
        for(int i=;i<=n;i++)
        {
            int j=i+l;
            if(j>n) break;
            dp[i][j][]=INF;
            for(int k=i+;k<j;k++)
            {
                memset(s1,,sizeof(s1));
                memset(s2,,sizeof(s2));
                memset(s3,,sizeof(s3));
                mul(a[i],a[j],a[k],s1);
                add(dp[i][k],dp[k][j],s2);
                add(s1,s2,s3);
                if(compare(dp[i][j],s3))
                    memcpy(dp[i][j],s3,sizeof(s3));
            }
        }
    print();
    return ;
}

请各位大佬斧正（反正我不认识斧正是什么意思）