codevs 1166 矩阵取数游戏

二次联通门 : codevs 1166 矩阵取数游戏

/*

    codevs 1166 矩阵取数游戏

    SB区间dp

    dp[l][r] = max (dp[l + 1][r] + number[l], dp[l][r - 1] + number[r]) * 2;
    不过要套高精

    我用的高精是全部封装好的
    可以像平时的int等类型用

    缺点就是慢。。。
    慢差不多1/3吧。。
*/
#include <iostream>
#include <cstdio>
#include <vector>
#include <iomanip>
#include <cassert>
#include <algorithm>

#define int64 long long

using namespace std;

#define Max 105

const int B = 1000000000;
const int L = 9;

inline int intcmp(int a, int b)
{
    if (a > b)
         return 1;
    else if (a < b)
          return -1;
    else
        return 0;
}

void read (int &now)
{
    now = 0;
    register char word = getchar ();
    while (word < ‘0‘ || word > ‘9‘)
        word = getchar ();
    while (word >= ‘0‘ && word <= ‘9‘)
    {
        now = now * 10 + word - ‘0‘;
        word = getchar ();
    }
}
struct BigInt
{

    vector<int> a;
    BigInt(){}
    BigInt(int n)
    {
        while (n > 0)
            a.push_back(n % B), n /= B;
    }

    BigInt(int64 n)
    {
        while (n > 0)
            a.push_back(n % B), n /= B;
    }
    inline void clr0()
    {
        while (!a.empty() && a.back() == 0)
            a.pop_back();
    }
    inline BigInt &operator+=(const BigInt &rhs)
    {
        a.resize(max(a.size(), rhs.a.size()));
        int t = 0;
        for (int i = 0; i < (int)rhs.a.size(); i++)
        {
            a[i] += rhs.a[i] + t;
            t = a[i] >= B;
            a[i] -= B & (-t);
        }
        for (int i = (int)rhs.a.size(); t != 0 && i < (int)a.size(); i++)
        {
            a[i] += t;
            t = a[i] >= B;
            a[i] -= B & (-t);
        }
        if (t != 0)
            a.push_back(t);
        return *this;
    }
    inline BigInt &operator-=(const BigInt &rhs)
    {
        a.resize(max(a.size(), rhs.a.size()));
        int t = 0;
        for (int i = 0; i < (int)rhs.a.size(); i++)
        {
            a[i] -= rhs.a[i] + t;
            t = a[i] < 0;
            a[i] += B & (-t);
        }
        for (int i = (int)rhs.a.size(); t != 0 && i < (int)a.size(); i++)
        {
            a[i] -= t;
            t = a[i] < 0;
            a[i] += B & (-t);
        }
        clr0();
        return *this;
    }
    inline BigInt &operator*=(const BigInt &rhs)
    {
        int na = (int)a.size();
        a.resize(na + rhs.a.size());
        for (int i = na - 1; i >= 0; i--)
        {
            int ai = a[i];
            int64 t = 0;
            a[i] = 0;
            for (int j = 0; j < (int)rhs.a.size(); j++)
            {
                t += a[i + j] + (int64)ai * rhs.a[j];
                a[i + j] = t % B;
                t /= B;
            }
            for (int j = (int)rhs.a.size(); t != 0 && i + j < (int)a.size(); j++)
            {
                t += a[i + j];
                a[i + j] = t % B;
                t /= B;
            }
            assert(t == 0);
        }
        clr0();
        return *this;
    }
    inline BigInt &operator/=(const BigInt &rhs)
    {
        return *this = div(rhs);
    }
    inline BigInt &operator%=(const BigInt &rhs)
    {
        return div(rhs), *this;
    }
    inline BigInt &shlb(int l = 1)
    {
        if (a.empty())
            return *this;
        a.resize(a.size() + l);
        for (int i = (int)a.size() - 1; i >= l; i--)
            a[i] = a[i - l];
        for (int i = 0; i < l; i++)
            a[i] = 0;
        return *this;
    }
    inline BigInt &shrb(int l = 1)
    {
        for (int i = 0; i < (int)a.size() - l; i++)
            a[i] = a[i + l];
        a.resize(max((int)a.size() - l, 0));
        return *this;
    }
    inline int cmp(const BigInt &rhs) const
    {
        if (a.size() != rhs.a.size())
            return intcmp(a.size(), rhs.a.size());
        for (int i = (int)a.size() - 1; i >= 0; i--)
            if (a[i] != rhs.a[i])
                return intcmp(a[i], rhs.a[i]);
        return 0;
    }
    inline BigInt div(const BigInt &rhs)
    {
        assert(!rhs.a.empty());
        if (rhs > *this)
            return 0;
        BigInt q, r;
        q.a.resize((int)a.size() - (int)rhs.a.size() + 1);
        for (int i = (int)a.size() - 1; i > (int)a.size() - (int)rhs.a.size(); i--)
        {
            r.shlb();
            r += a[i];
        }
        for (int i = (int)a.size() - (int)rhs.a.size(); i >= 0; i--)
        {
            r.shlb();
            r += a[i];
            if (r.cmp(rhs) < 0)
                q.a[i] = 0;
            else
            {
                int le = 0, ri = B;
                while (le != ri)
                {
                    int mi = (le + ri) / 2;
                    if ((rhs * mi).cmp(r) <= 0)
                        le = mi + 1;
                    else
                        ri = mi;
                }
                q.a[i] = le - 1;
                r -= rhs * q.a[i];
            }
        }
        q.clr0();
        *this = r;
        return q;
    }
    friend inline BigInt operator+(const BigInt &lhs, const BigInt &rhs)
    {
        BigInt res = lhs;
        return res += rhs;
    }
    friend inline BigInt operator-(const BigInt &lhs, const BigInt &rhs)
    {
        BigInt res = lhs;
        return res -= rhs;
    }
    friend inline BigInt operator*(const BigInt &lhs, const BigInt &rhs)
    {
        BigInt res = lhs;
        return res *= rhs;
    }
    friend inline BigInt operator/(const BigInt &lhs, const BigInt &rhs)
    {
        BigInt res = lhs;
        return res.div(rhs);
    }
    friend inline BigInt operator%(const BigInt &lhs, const BigInt &rhs)
    {
        BigInt res = lhs;
        return res.div(rhs), res;
    }
    friend inline ostream &operator<<(ostream &out, const BigInt &rhs)
    {
        if (rhs.a.size() == 0)
            out << "0";
        else
        {
            out << rhs.a.back();
            for (int i = (int)rhs.a.size() - 2; i >= 0; i--)
                out << setfill(‘0‘) << setw(L) << rhs.a[i];
        }
        return out;
    }
    friend inline bool operator<(const BigInt &lhs, const BigInt &rhs)
    {
        return lhs.cmp(rhs) < 0;
    }
    friend inline bool operator<=(const BigInt &lhs, const BigInt &rhs)
    {
        return lhs.cmp(rhs) <= 0;
    }
    friend inline bool operator>(const BigInt &lhs, const BigInt &rhs)
    {
        return lhs.cmp(rhs) > 0;
    }
    friend inline bool operator>=(const BigInt &lhs, const BigInt &rhs)
    {
        return lhs.cmp(rhs) >= 0;
    }
    friend inline bool operator==(const BigInt &lhs, const BigInt &rhs)
    {
        return lhs.cmp(rhs) == 0;
    }
    friend inline bool operator!=(const BigInt &lhs, const BigInt &rhs)
    {
        return lhs.cmp(rhs) != 0;
    }
};

inline BigInt BigInt_max (BigInt a, BigInt b)
{
    return a > b ? a : b;
}

int N, M;

int number[Max];

BigInt dp[Max][Max];
BigInt Answer;

int main (int argc, char *argv[])
{
    ios :: sync_with_stdio (false);
    read (N);
    read (M);

    for (int i = 1; i <= N; i ++)
    {
        for (register int pos = 1; pos <= M; pos ++)
            for (register int __pos = 1; __pos <= M; __pos ++)
                dp[pos][__pos] = 0;
        for (register int pos = 1; pos <= M; pos ++)
        {
            read (number[pos]);
            dp[pos][pos] = 2 * number[pos];
        }

        for (register int size = 1; size < M; size ++)
            for (register int l = 1; l + size <= M; l ++)
            {
                register int r = l + size;
                dp[l][r] = 2 * BigInt_max (dp[l + 1][r] + number[l], dp[l][r - 1] + number[r]);
            }
        Answer += dp[1][M];
    }
    cout << Answer;
    return 0;
}
时间: 2024-11-06 09:42:27

codevs 1166 矩阵取数游戏的相关文章

1166 矩阵取数游戏[区间dp+高精度]

1166 矩阵取数游戏 2007年NOIP全国联赛提高组 时间限制: 1 s 空间限制: 128000 KB 题目等级 : 黄金 Gold 题解 题目描述 Description [问题描述]帅帅经常跟同学玩一个矩阵取数游戏:对于一个给定的n*m 的矩阵,矩阵中的每个元素aij均为非负整数.游戏规则如下:1. 每次取数时须从每行各取走一个元素,共n个.m次后取完矩阵所有元素:2. 每次取走的各个元素只能是该元素所在行的行首或行尾:3. 每次取数都有一个得分值,为每行取数的得分之和,每行取数的得分

矩阵取数游戏 NOIP 2007

2016-05-31 17:26:45 题目链接: NOIP 2007 矩阵取数游戏(Codevs) 题目大意: 给定一个矩阵,每次在每一行的行首或者行尾取一个数乘上2^次数,求取完最多获得的分数 解法: 动态规划 DP[i][j]表示当前行i位置到j位置获得的最大分数 转移方程: DP[i][j]=max(DP[i+1][j]+a[i]*2^x,DP[i][j-1]+a[i]*2^x); 需要注意的地方: 1.M和N给的范围略坑啊,2^80只有悄悄不说话了,所以直接上了100000的压位,感觉

矩阵取数游戏

题目描述 Description [问题描述]帅帅经常跟同学玩一个矩阵取数游戏:对于一个给定的n*m 的矩阵,矩阵中的每个元素aij均为非负整数.游戏规则如下:1. 每次取数时须从每行各取走一个元素,共n个.m次后取完矩阵所有元素:2. 每次取走的各个元素只能是该元素所在行的行首或行尾:3. 每次取数都有一个得分值,为每行取数的得分之和,每行取数的得分= 被取走的元素值*2i,其中i 表示第i 次取数(从1 开始编号):4. 游戏结束总得分为m次取数得分之和.帅帅想请你帮忙写一个程序,对于任意矩

P1005矩阵取数游戏

题目描述 帅帅经常跟同学玩一个矩阵取数游戏:对于一个给定的n*m的矩阵,矩阵中的每个元素aij均为非负整数.游戏规则如下: 1.每次取数时须从每行各取走一个元素,共n个.m次后取完矩阵所有元素: 2.每次取走的各个元素只能是该元素所在行的行首或行尾: 3.每次取数都有一个得分值,为每行取数的得分之和,每行取数的得分 = 被取走的元素值*2^i,其中i表示第i次取数(从1开始编号): 4.游戏结束总得分为m次取数得分之和. 帅帅想请你帮忙写一个程序,对于任意矩阵,可以求出取数后的最大得分. 输入输

NOIP2007 矩阵取数游戏

题目描述 Description [问题描述] 帅帅经常跟同学玩一个矩阵取数游戏:对于一个给定的n*m 的矩阵,矩阵中的每个元素aij均 为非负整数.游戏规则如下: 1. 每次取数时须从每行各取走一个元素,共n个.m次后取完矩阵所有元素: 2. 每次取走的各个元素只能是该元素所在行的行首或行尾: 3. 每次取数都有一个得分值,为每行取数的得分之和,每行取数的得分= 被取走的元素值*2i, 其中i 表示第i 次取数(从1 开始编号): 4. 游戏结束总得分为m次取数得分之和. 帅帅想请你帮忙写一个

矩阵取数游戏洛谷p1005

题目描述 帅帅经常跟同学玩一个矩阵取数游戏:对于一个给定的n*m的矩阵,矩阵中的每个元素aij均为非负整数.游戏规则如下: 1.每次取数时须从每行各取走一个元素,共n个.m次后取完矩阵所有元素: 2.每次取走的各个元素只能是该元素所在行的行首或行尾: 3.每次取数都有一个得分值,为每行取数的得分之和,每行取数的得分 = 被取走的元素值*2^i,其中i表示第i次取数(从1开始编号): 4.游戏结束总得分为m次取数得分之和. 帅帅想请你帮忙写一个程序,对于任意矩阵,可以求出取数后的最大得分. 输入输

洛谷 P1005 矩阵取数游戏 题解

此文为博主原创题解,转载时请通知博主,并把原文链接放在正文醒目位置. 题目链接:https://www.luogu.org/problem/show?pid=1005 题目描述 帅帅经常跟同学玩一个矩阵取数游戏:对于一个给定的n*m的矩阵,矩阵中的每个元素aij均为非负整数.游戏规则如下: 1.每次取数时须从每行各取走一个元素,共n个.m次后取完矩阵所有元素: 2.每次取走的各个元素只能是该元素所在行的行首或行尾: 3.每次取数都有一个得分值,为每行取数的得分之和,每行取数的得分 = 被取走的元

P1005 矩阵取数游戏

P1005 矩阵取数游戏 题目描述 帅帅经常跟同学玩一个矩阵取数游戏:对于一个给定的n*m的矩阵,矩阵中的每个元素aij均为非负整数.游戏规则如下: 1.每次取数时须从每行各取走一个元素,共n个.m次后取完矩阵所有元素: 2.每次取走的各个元素只能是该元素所在行的行首或行尾: 3.每次取数都有一个得分值,为每行取数的得分之和,每行取数的得分 = 被取走的元素值*2^i,其中i表示第i次取数(从1开始编号): 4.游戏结束总得分为m次取数得分之和. 帅帅想请你帮忙写一个程序,对于任意矩阵,可以求出

洛谷 P1005 矩阵取数游戏

题目描述 帅帅经常跟同学玩一个矩阵取数游戏:对于一个给定的n*m的矩阵,矩阵中的每个元素aij均为非负整数.游戏规则如下: 1.每次取数时须从每行各取走一个元素,共n个.m次后取完矩阵所有元素: 2.每次取走的各个元素只能是该元素所在行的行首或行尾: 3.每次取数都有一个得分值,为每行取数的得分之和,每行取数的得分 = 被取走的元素值*2^i,其中i表示第i次取数(从1开始编号): 4.游戏结束总得分为m次取数得分之和. 帅帅想请你帮忙写一个程序,对于任意矩阵,可以求出取数后的最大得分. 输入输