SGU167（dp，输出方案）

题目描述
在N×MN×M的矩阵中，每个格子有一个权值，要求寻找一个包含K个格子的凸连通块（连通块中间没有空缺，并且轮廓是凸的，如右图所示），使这个连通块中的格子的权值和最大。求出这个最大的权值和，并给出连通块的具体方案。

输入格式
在输入的第一行上有3个整数N，M，K（1 <= N，M <= 15,0 <= K <= N * M）。 接下来的N行包含M个整数，它们是该正方形上的石油资源数量。 这些数字中的每一个都在0到1000的范围内。

输出格式
在输出的第一行，写入字符串“Oil：X”，其中整数X是可由A国控制的最大油数。 接下来你应该输出K对数字—将被A国占据的正方形坐标。 第一个坐标是行数（从上到下，从1开始），第二个是列数（从左到右，从1开始）。

样例数据
input
2 3 4
10 20 30
40 2 3
output
Oil : 100
1 1
1 2
1 3
2 1
数据规模与约定
时间限制：0.75s
空间限制：65MB

---

首先一个关键的部分就是凸包如何处理
其实放在这个图上一个凸包就是对于每行的区间$[l_i,r_i]$
满足$l_i$先递减再递增，$r_i$先递增再递减
然后就可$dp$了
$f_{i,j,l,r,f1(1/0),f2(1/0)}$表示前$i$行放了$j$个位置，当前行放了$[l,r]$的区间
放了之后左右端点分别是增(0)或减(1)的
若$f_1=1,f_2=0$

$$f_{i,j,l,r,f1,f2}=\max_{l≤p≤q≤r}\{f_{i-1,j-l+r-1,p,q,1,0}\}+v_{l,r}$$
若 $f_1=1,f_2=1$
$$f_{i,j,l,r,f1,f2}=\max_{l≤p≤r≤q}\{ \max_{x∈\{0,1\}}\{f_{i-1,j-l+r-1,p,q,1,x}\} \}+v_{l,r}$$
若 $f_1=0,f_2=0$
$$f_{i,j,l,r,f1,f2}=\max_{p≤l≤q≤r}\{ \max_{x∈\{0,1\}}\{f_{i-1,j-l+r-1,p,q,x,0}\} \}+v_{l,r}$$
若 $f_1=0,f_2=1$
$$f_{i,j,l,r,f1,f2}=\max_{p≤l≤r≤q}\{ \max_{x∈\{0,1\},y∈\{0,1\}}\{f_{i-1,j-l+r-1,p,q,x,y}\} \}+v_{l,r}$$
转移即可
输出方案再开个数组记录哪转移过来即可

#include<bits/stdc++.h>
using namespace std;
#define rep(i,j,k) for(int i = j;i <= k;++i)
#define repp(i,j,k) for(int i = j;i >= k;--i)
#define rept(i,x) for(int i = linkk[x];i;i = e[i].n)
#define P pair<int,int>
#define Pil pair<int,ll>
#define Pli pair<ll,int>
#define Pll pair<ll,ll>
#define pb push_back 
#define pc putchar
#define mp make_pair
#define file(k) memset(k,0,sizeof(k))
#define fr first
#define se second
#define ll long long
namespace fastIO{
    #define BUF_SIZE 100000
    #define OUT_SIZE 100000
    bool IOerror = 0;
    inline char nc(){
        static char buf[BUF_SIZE],*p1 = buf+BUF_SIZE, *pend = buf+BUF_SIZE;
        if(p1 == pend){
            p1 = buf; pend = buf+fread(buf, 1, BUF_SIZE, stdin);
            if(pend == p1){ IOerror = 1; return -1;}
        }
        return *p1++;
    }
    inline bool blank(char ch){return ch==' '||ch=='\n'||ch=='\r'||ch=='\t';}
    inline void read(int &x){
        bool sign = 0; char ch = nc(); x = 0;
        for(; blank(ch); ch = nc());
        if(IOerror)return;
        if(ch == '-') sign = 1, ch = nc();
        for(; ch >= '0' && ch <= '9'; ch = nc()) x = x*10+ch-'0';
        if(sign) x = -x;
    }
    inline void read(ll &x){
        bool sign = 0; char ch = nc(); x = 0;
        for(; blank(ch); ch = nc());
        if(IOerror) return;
        if(ch == '-') sign = 1, ch = nc();
        for(; ch >= '0' && ch <= '9'; ch = nc()) x = x*10+ch-'0';
        if(sign) x = -x;
    }
    #undef OUT_SIZE
    #undef BUF_SIZE
};
using namespace fastIO;
int n , m , k;
int v[16][16];
struct node{  
    int l,r,f1,f2;
};
int f[16][226][16][16][2][2];
node la[16][226][16][16][2][2];
node final;
P ans;
int Get[20][20] , ii[20];
void umax(int i,int j,node now,node last)
{
    int l = now.l,r = now.r,f1 = now.f1,f2 = now.f2;
    int ff1 = last.f1,ff2 = last.f2;
    int p = last.l,q = last.r;
    int no = r-l+1;
    if(f[i-1][j-no][p][q][ff1][ff2]+v[i][r]-v[i][l-1] > f[i][j][l][r][f1][f2])
        f[i][j][l][r][f1][f2] = f[i-1][j-no][p][q][ff1][ff2]+v[i][r]-v[i][l-1],
        la[i][j][l][r][f1][f2].l = p,la[i][j][l][r][f1][f2].r = q,
        la[i][j][l][r][f1][f2].f1 = ff1,la[i][j][l][r][f1][f2].f2 = ff2;
    return;
}
void New(P &a,int i,int j,node now)
{
    if(f[i][j][now.l][now.r][now.f1][now.f2] > a.fr)
        a.fr = f[i][j][now.l][now.r][now.f1][now.f2],
        a.se = i,
        final = now;
}
char s[1010],c[1010];
int main()
{
    read(n);read(m);read(k);
    rep(i,1,n) rep(j,1,m) read(v[i][j]) , v[i][j] += v[i][j-1];
    rep(i,1,n) rep(l,1,m) rep(r,l,m)
    {
        if(r-l+1<=k)f[i][r-l+1][l][r][1][0] = v[i][r] - v[i][l-1];
        rep(j,r-l+2,k)
        {
            rep(p,l,r) rep(q,l,r) //l<=p<=q<=r
            umax(i,j,(node){l,r,1,0},(node){p,q,1,0});
            rep(p,l,r) rep(q,r,m)//l<=p<=r<=q
            umax(i,j,(node){l,r,1,1},(node){p,q,1,0}),
            umax(i,j,(node){l,r,1,1},(node){p,q,1,1}); 
            rep(p,1,l) rep(q,l,r)//p<=l<=q<=r
            umax(i,j,(node){l,r,0,0},(node){p,q,0,0}),
            umax(i,j,(node){l,r,0,0},(node){p,q,1,0});
            rep(p,1,l) rep(q,r,m)
            umax(i,j,(node){l,r,0,1},(node){p,q,0,0}),
            umax(i,j,(node){l,r,0,1},(node){p,q,0,1}),
            umax(i,j,(node){l,r,0,1},(node){p,q,1,0}),
            umax(i,j,(node){l,r,0,1},(node){p,q,1,1});
            if(j == k) New(ans,i,j,(node){l,r,0,0}),
                       New(ans,i,j,(node){l,r,0,1}),
                       New(ans,i,j,(node){l,r,1,0}),
                       New(ans,i,j,(node){l,r,1,1});
        }
    }
    printf("Oil : %d\n",ans.fr);
    node now = final;
    int i = ans.se;
    int j = k;
    int o = 0;
    while(now.l <= now.r && now.l != 0 && now.r != 0)
    {
        ++o;ii[o] = i;
        rep(k,now.l,now.r) Get[o][++Get[o][0]] = k;
        int tmp = now.r - now.l + 1;
        now = la[i][j][now.l][now.r][now.f1][now.f2];
        i--;
        j -= tmp;
    }
    repp(i,o,1) rep(j,1,Get[i][0])
    printf("%d %d\n",ii[i],Get[i][j]);
    return 0;
}