DPトピック-cfの711C

CF 711 C
、今着色後のクレームには、今、着色、n個の木の着色の色、価格を着色ツリーだけでなく、木のn個、初期の色を与える：イタリア小片に木が同じ色のセクションに分割された各連続として定義され、連続切片をK、必要kの最小コスト分割を見つけます。
私たちは、[j]は[[i]のDPを作る k]はi番目の時、色jを表す行う、kは最小コストの範囲を分割しました。
その結果、ツリーの現在の色は、正面のみから転送された場合、（precolor == nowcolor）DP [あたかも I] [J] [K] =分（DP [I-1]〜[J] [K 】）
他DP [I] [J] [K] =分（DP [-I 1] [precolor] [-K 1]）。
この1つは値を持たないためので、それはそれは、現在の色を列挙するために必要となる現在の色jが決定されていないです。
nowcolor用：
DP [I] [nowcolor] [K] =分（DP [I-1] [nowcolor] [K]）または分（DP [precolor] [K-1] +コスト[Iは、[I-1] ] [nowcolor]

    #include<iostream>
    #include<cstring>
    #include<algorithm>
    #include<cmath>
    #include<cstdlib>
    #include<climits>
    #include<stack>
    #include<vector>
    #include<queue>
    #include<set>
    #include<map>
    //#include<regex>
    #include<cstdio>
    #define up(i,a,b)  for(int i=a;i<b;i++)
    #define dw(i,a,b)  for(int i=a;i>b;i--)
    #define upd(i,a,b) for(int i=a;i<=b;i++)
    #define dwd(i,a,b) for(int i=a;i>=b;i--)
    //#define local
    typedef long long ll;
    const double esp = 1e-6;
    const double pi = acos(-1.0);
    const int INF = 0x3f3f3f3f;
    const int inf = 1e9;
    using namespace std;
    int read()
    {
        char ch = getchar(); int x = 0, f = 1;
        while (ch<'0' || ch>'9') { if (ch == '-')f = -1; ch = getchar(); }
        while (ch >= '0' && ch <= '9') { x = x * 10 + ch - '0'; ch = getchar(); }
        return x * f;
    }
    typedef pair<int, int> pir;
    #define lson l,mid,root<<1
    #define rson mid+1,r,root<<1|1
    #define lrt root<<1
    #define rrt root<<1|1
    ll dp[105][105][105];
    ll c[105];
    ll p[105][105];
    int n, m, k;
    int main()
    {
        n = read(), m = read(), k = read();
        upd(i, 1, n)
        {
            c[i] = read();
        }
        upd(i, 1, n)
        {
            upd(j, 1, m)
            {
                p[i][j] = read();
            }
        }
        upd(i, 0, 101)upd(j, 0, 101)upd(t, 0, 101)dp[i][j][t] = 1e18;
        if (c[1])
        {
            dp[1][c[1]][1] = 0;
        }
        else
        {
            upd(i, 1, m)dp[1][i][1] = p[1][i];
        }
        upd(i, 2, n)
        {
            if (c[i])
            {
                upd(t, 1, k)
                {
                    upd(j, 1, m)
                    {
                        if (t > i)break;
                        if (c[i] == j)dp[i][j][t] = min(dp[i - 1][j][t], dp[i][j][t]);
                        else dp[i][c[i]][t] = min(dp[i][c[i]][t], dp[i - 1][j][t - 1]);
                    }
                }
            }
            else
            {
                upd(j, 1, m)
                {
                    upd(t, 1, k)
                    {
                        upd(o, 1, m)
                        {
                            if (t > i)break;
                            if (o == j)dp[i][j][t] = min(dp[i - 1][o][t], dp[i][j][t]);
                            else dp[i][j][t] = min(dp[i][j][t], dp[i - 1][o][t - 1]);
                        }
                        dp[i][j][t] += p[i][j];
                    }
                }
            }
        }
        ll ans = 1e18;
        upd(i, 1, m)
        {
            ans = min(ans, dp[n][i][k]);
        }
        if (ans == 1e18)cout << "-1" << endl;
        else cout << ans << endl;
        return 0;
    }