【NOIP2007】矩阵取数题解

题面

【问题描述】

帅帅经常更同学玩一个矩阵取数游戏：对于一个给定的n*m的矩阵，矩阵中的每个元素$a_{ij}$据为非负整数。游戏规则如下：
1. 每次取数时须从每行各取走一个元素，共n个。m次后取完矩阵所有的元素；
2. 每次取走的各个元素只能是该元素所在行的行首或行尾；
3. 每次取数都有一个得分值，为每行取数的得分之和；每行取数的得分 = 被取走的元素值*2i，其中i表示第i次取数（从1开始编号）；
4. 游戏结束总得分为m次取数得分之和。
帅帅想请你帮忙写一个程序，对于任意矩阵，可以求出取数后的最大得分。

【输入】

输入文件game.in包括n+1行；
第一行为两个用空格隔开的整数n和m。
第2~n+1行为n*m矩阵，其中每行有m个用单个空格隔开

【输出】

输出文件game.out仅包含1行，为一个整数，即输入矩阵取数后的最大的分。

【输入输出样例1】

game.in
2 3
1 2 4
3 4 2
game.out
82

【输入输出样例1解释】

第1次：第一行取行首元素，第二行取行尾元素，本次的氛围1*21+2*21=6
第2次：两行均取行首元素，本次得分为2*22+3*22=20
第3次：得分为3*23+4*23=56。总得分为6+20+56=82

【输入输出样例2】

game.in
1 4
4 5 0 5
game.out
122

【输入输出样例3】

game.in
2 10
96 56 54 46 86 12 23 88 80 43
16 95 18 29 30 53 88 83 64 67
game.out
316994

【限制】

60%的数据满足：$1<=n, m<=30$，答案不超过$10^{16}$
100%的数据满足：$1<=n, m<=80$，$0<=a_{ij}<=1000$

题解

分析

由于不能把数放回去，因此可以用动态规划解决；因为行与行之间没有关联，于是我们逐行处理。
用$f[i][j]$表示当前行i~j个数的得分最大值，a[i]表示第i位上的数。
则：

$$f[i][j] = max\{f[i+1][j]+a[i]*2^x, f[i][j-1]+a[j]*2^x\}$$
注意高精度

代码

#include<stdio.h>
#include<stdlib.h>
#include<string.h>

typedef struct node{
    int a[40];
}hugeint;

hugeint f[81][81],k1,k2,a[81],max;

int n, m;

hugeint add(hugeint k,hugeint l)
{
    hugeint j;
    memset(j.a, 0, sizeof(j.a));

    for(int i = 1; i <= k.a[0] || i<=l.a[0]; i++)
    {
        j.a[i] = k.a[i] + l.a[i];
    }

    j.a[0] = k.a[0] > l.a[0] ? k.a[0] : l.a[0];

    for(int i = 1; i <= j.a[0]; i++)
    {
        j.a[i+1] += j.a[i]/10;
        j.a[i] %= 10;
    }

    if(j.a[j.a[0]+1]>0)
        j.a[0]++;

    return j;
}

void init()
{
    memset(a,0,sizeof(a));
    for(int i=1;i<=m;i++)
    {
        scanf("%d",&a[i].a[1]);
        a[i].a[0]=1;
        while(a[i].a[a[i].a[0]]>=10)
        {
            a[i].a[a[i].a[0]+1]=a[i].a[a[i].a[0]]/10;
            a[i].a[a[i].a[0]]%=10;
            a[i].a[0]++;
        }
    }
}

hugeint cmp(hugeint k,hugeint l)
{
    if(k.a[0]>l.a[0])
        return k;

    if(k.a[0]<l.a[0])
        return l;


    if(k.a[0]==l.a[0])
    {
        for(int i=k.a[0];i>=1;i--)
        {
            if(k.a[i]>l.a[i])
            {
                return k;
            }
            if(k.a[i]<l.a[i])
            {
                return l;
            }
        }
    }



    return l;
}

void solve()
{
    memset(f,0,sizeof(f));
    for(int i=1;i<=m;i++)
    {
        f[i][i]=add(a[i],a[i]);
    }

    for(int i = m; i >= 1; i--)
    {
        for(int j = i+1; j<=m; j++)
        {
            k1 = add(f[i+1][j],a[i]);
            k2 = add(f[i][j-1],a[j]);
            f[i][j] = cmp(k1,k2);
            f[i][j] = add(f[i][j],f[i][j]);
        }
    }
    max=add(max,f[1][m]);
}

void output()
{
    for (int i = max.a[0]; i>=1; i--)
    {
        printf("%d",max.a[i]);
    }
}

int main()
{
    freopen("game.in","r",stdin);
    freopen("game.out","w",stdout);
    scanf("%d%d",&n,&m);
    max.a[0]=1;
    for(int i = 1; i <= n; i++)
    {
        init();
        solve();
    }
    output();
    return 0;
}

---

是我说谎，第六篇还是来了，累死了