P1005 矩阵取数游戏

最新推荐文章于 2022-09-17 17:18:18 发布
最新推荐文章于 2022-09-17 17:18:18 发布
阅读量540 收藏 1
点赞数
分类专栏： DP 文章标签： 高精度 动态规划DP
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_34786702/article/details/52145066
版权
题目描述

帅帅经常跟同学玩一个矩阵取数游戏：对于一个给定的n*m的矩阵，矩阵中的每个元素aij均为非负整数。游戏规则如下：

1.每次取数时须从每行各取走一个元素，共n个。m次后取完矩阵所有元素；

2.每次取走的各个元素只能是该元素所在行的行首或行尾；

3.每次取数都有一个得分值，为每行取数的得分之和，每行取数的得分 = 被取走的元素值*2^i，其中i表示第i次取数（从1开始编号）；

4.游戏结束总得分为m次取数得分之和。

帅帅想请你帮忙写一个程序，对于任意矩阵，可以求出取数后的最大得分。

输入输出格式

输入格式：
输入文件game.in包括n+1行：

第1行为两个用空格隔开的整数n和m。

第2~n+1行为n*m矩阵，其中每行有m个用单个空格隔开的非负整数。

数据范围：

60%的数据满足：1<=n, m<=30，答案不超过10^16

100%的数据满足：1<=n, m<=80，0<=aij<=1000

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

输入输出样例

输入样例#1：
2 3
1 2 3
3 4 2

输出样例#1：
82

说明

NOIP 2007 提高第三题

题解

从DP角度来说，是一道简单题；

对于每一层，层与层之间的各种状态互不影响。所以，问题简单了，求每一层的解 其sigma为答案

f[i][j] 表示 在区间[i][j]取数的最优解 f[i][j]=max(f[i+1][j]+a[i](2^当前次幂),f[i][j-1]+a[j](2^当前次幂))

可以记忆化搜索

然而本题最终数值过大 需要高精 不用高精的能过60

60分代码

#include <iostream>
#include <cstring>
using namespace std;
unsigned long long n,m,f[1111][1111],map[1111][1111],ans;
unsigned long long dfs(int x,int l,int r,unsigned long long s)
{
    if(f[l][r]) return f[l][r];
    if(l==r)
    {
        f[l][r]=map[x][l]*s;
        return f[l][r];
    }
    f[l][r]=max(dfs(x,l,r-1,s*2)+map[x][r]*s,dfs(x,l+1,r,s*2)+map[x][l]*s);
    return f[l][r];
}
int main()
{
    cin>>n>>m;
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            cin>>map[i][j];
        }
    }
    for(int i=1;i<=n;i++)
    {    
        memset(f,0,sizeof(f));
        ans+=dfs(i,1,m,2);
    }
    cout<<ans<<endl;

}

高精的话 则有f[i][j][k] 第三维相当于高精度的每一位

并且我进行了压位（压四位）；

并在开始预处理出了2的n次幂

下面是既恶心又恶心而且恶心的高精。。。。。

#include <iostream>
#include <cstring>
using namespace std;
#define pp 10000
#define ll long long
int n,m,f[111][111][111],map[1111][1111],ans[111],z[111][1111];
int gd1[111],gd2[111],gd3[111],gd4[111],gd5[111];
int bo[111][111];
int init()//初始化到 2的m次幂 
{
    z[1][1]=2;z[1][0]=1;
    for(int ii=2;ii<=m;ii++)
    {
        int tot=z[ii-1][0]+1;
        for(int i=1;i<=z[ii-1][0];i++)
        {
            z[ii][i]=z[ii-1][i]*2;
        }
        for(int i=1;i<=tot;i++)
        {
            z[ii][i+1]+=z[ii][i]/pp;
            z[ii][i]%=pp;
        }
        while(!z[ii][tot]&&tot>1) tot--;
        z[ii][0]=tot;
    }
}
int print(int * a)//输出数组a中的数   可用于debug 
{
    for(int i=a[0];i>=1;i--)
    {
        if(i==a[0])
        {
            //printf("%d",a[i]);
            cout<<a[i];
            continue;
        }
        int now=a[i];
        int xx[6];
        xx[4]=now%10;
        now/=10;
        xx[3]=now%10;
        now/=10;
        xx[2]=now%10;
        now/=10;
        xx[1]=now%10;
        now/=10;
        cout<<xx[1]<<xx[2]<<xx[3]<<xx[4];
        //printf("%d%d%d%d%d",xx[1],xx[2],xx[3],xx[4]);
    }

}
int gjj(int *a,int *b,int *c)//把b数组与c数组的数相加  赋给a数组 
{
    for(int i=1;i<=max(b[0],c[0]);i++)
    {
        a[i]=b[i]+c[i];
    }
    a[0]=max(b[0],c[0])+1;
    for(int i=1;i<=a[0];i++)
    {
        a[i+1]+=a[i]/pp;
        a[i]%=pp;
    }
    while(!a[a[0]]&&a[0]>1) a[0]--;/****要a[0]大于1   有一组数值ans==0   不大于1 会无输出。。。。*********/ 
}
int gjc(int *a,int x,int b)//把x与2的b次方的乘积赋值给a数组 
{
    //memset()

    a[0]=z[b][0]+4;
    //for(int i=1;i<=a[0];i++)
    //  a[i]=0;
    //cout<<a[2]<<"   a 2    ";
    for(int i=1;i<=z[b][0];i++)
    {
        a[i]=z[b][i]*x;
        //cout<<a[i]<<' '<<i<<' '<<z[b][i]<<' '<<x<<' '<<i<<' '<<z[b][0]<<' '<<a[2]<<endl;
    }
    for(int i=1;i<=a[0];i++)
    {
    //  cout<<a[i]<<' '<<i<<' '<<pp<<' '<<a[i+1]<<endl;
        a[i+1]+=a[i]/pp;
        a[i]%=pp;
    }
    while(!a[a[0]]&&a[0]>1) a[0]--;
    //cout<<"fdfdfdfdfd"<<endl;
    //for(int i=a[0];i>=1;i--)
    //  cout<<a[i];
    //cout<<endl<<x<<' '<<b<<' '<<"fdfdfdf"<<' '<<z[b][1]<<' '<<z[b][1]*x<<endl;
}
int max(int *a,int *b,int *c)//把b数组与c数组中数值大的赋值给a数组 
{
    int pd=0;
    if(b[0]>c[0])   pd=1;
    if(b[0]<c[0])   pd=2;
    if(b[0]==c[0])
    {
        int h=b[0];
        while(h>=1)
        {
            if(b[h]>c[h])
            {
                pd=1;break;
            }
            if(b[h]<c[h])
            {
                pd=2;
                break;
            }
            h--;
        }
    }
    if(pd==1)//DEBUG!!!!!!!!  pd=1!!!!!!!@!$@!$@!
        for(int i=0;i<=b[0];i++)
            a[i]=b[i];
    else
        for(int i=0;i<=c[0];i++)
            a[i]=c[i];
    while(!a[a[0]]&&a[0]>1) a[0]--;
    /*print(b);
    cout<<"~~~"<<' '<<b[0]<< ' '<<c[0]<<' '<<pd<<' '<<b[1]<<' '<<c[1]<<endl;
    print(a);
    cout<<"~~~~~~~~~!!"<< ' ';
    print(c);*/
}

int dfs(int x,int l,int r,int s)//记忆化搜索 
{
    if(bo[l][r]) return 0;//记忆化 
    bo[l][r]=1;
    if(l==r)
    {
        gjc(f[l][r],map[x][l],s);/****相当于下面一行的状态转移方程*******/ 
    //  f[l][r]=map[x][l]*s;
    //  print(f[l][r]);
        //cout<<f[l][r][1]<<' '<<l<<' '<<r<<' '<<map[x][l]<<' '<<s<<' '<<z[s][1]<<endl;
    //  return f[l][r];
        return 0;
    }
    memset(gd1,0,sizeof(gd1));/*******************/
    memset(gd2,0,sizeof(gd2));/*初始化很重要**/ 
    memset(gd3,0,sizeof(gd3));/**！！！！！*******/
    memset(gd4,0,sizeof(gd4));/*******************/
    dfs(x,l,r-1,s+1);
    dfs(x,l+1,r,s+1);
    memset(gd1,0,sizeof(gd1));/*******************/
    memset(gd2,0,sizeof(gd2));/*初始化很重要->AC**/ 
    memset(gd3,0,sizeof(gd3));/**！！！！！*******/
    memset(gd4,0,sizeof(gd4));//不能只在上面初始化    因为中间有dfs 和递归过程  期间 gd1 2 3 4 数组的值早已改变！！！！ 
    gjc(gd1,map[x][r],s);           ///***********************/
    gjc(gd2,map[x][l],s);           // ~~~~**********************/
    gjj(gd3,f[l][r-1],gd1);         //~~~~**********************/
    gjj(gd4,f[l+1][r],gd2);         //~~~ *********************/
    max(f[l][r],gd3,gd4);           //~~~此五行 相当于下面 无高精的状态转移方程  

    //f[l][r]=max(dfs(x,l,r-1,s*2)+map[x][r]*s,dfs(x,l+1,r,s*2)+map[x][l]*s);
    //return f[l][r];
}
int main()
{
    cin>>n>>m;
    init();//预处理 2的n次方
    //gjc(f[1][1],68,7);
    /*map[1][1]=99;
    dfs(1,1,1,80);*/
    /*for(int i=40;i>=1;i--)
    {
        int now=f[1][1][i];
        int xx[6];
        xx[4]=now%10;
        now/=10;
        xx[3]=now%10;
        now/=10;
        xx[2]=now%10;
        now/=10;
        xx[1]=now%10;
        now/=10;
        cout<<xx [1]<<xx [2]<<xx [3]<<xx [4];
    }*/
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            cin>>map[i][j];
        }
    }
    for(int i=1;i<=n;i++)
    {   
        memset(f,0,sizeof(f));// 
        memset(bo,0,sizeof(bo));//对每一层for i的循环  需要初始化 用来记忆的bo数组  以及每层的f数组（滚动） 
        dfs(i,1,m,1);
        memset(gd5,0,sizeof(gd5));
        for(int ii=0;ii<=ans[0];ii++)
            gd5[ii]=ans[ii];
        gjj(ans,gd5,f[1][m]);
        /*print(f[1][1]);
        print(f[2][2]);
        print(f[3][3]);*/
    //  print(f[1][3]);
    //  cout<<endl; 
    }
    print(ans);
}

无debug的注释的

#include <iostream>
#include <cstring>
using namespace std;
#define pp 10000
#define ll long long
int n,m,f[111][111][111],map[1111][1111],ans[111],z[111][1111];
int gd1[111],gd2[111],gd3[111],gd4[111],gd5[111];
int bo[111][111];
int init()//初始化到 2的m次幂 
{
    z[1][1]=2;z[1][0]=1;
    for(int ii=2;ii<=m;ii++)
    {
        int tot=z[ii-1][0]+1;
        for(int i=1;i<=z[ii-1][0];i++)
        {
            z[ii][i]=z[ii-1][i]*2;
        }
        for(int i=1;i<=tot;i++)
        {
            z[ii][i+1]+=z[ii][i]/pp;
            z[ii][i]%=pp;
        }
        while(!z[ii][tot]&&tot>1) tot--;
        z[ii][0]=tot;
    }
}
int print(int * a)//输出数组a中的数   可用于debug 
{
    for(int i=a[0];i>=1;i--)
    {
        if(i==a[0])
        {
            cout<<a[i];
            continue;
        }
        int now=a[i];
        int xx[6];
        xx[4]=now%10;
        now/=10;
        xx[3]=now%10;
        now/=10;
        xx[2]=now%10;
        now/=10;
        xx[1]=now%10;
        now/=10;
        cout<<xx[1]<<xx[2]<<xx[3]<<xx[4];
    }

}
int gjj(int *a,int *b,int *c)//把b数组与c数组的数相加  赋给a数组 
{
    for(int i=1;i<=max(b[0],c[0]);i++)
    {
        a[i]=b[i]+c[i];
    }
    a[0]=max(b[0],c[0])+1;
    for(int i=1;i<=a[0];i++)
    {
        a[i+1]+=a[i]/pp;
        a[i]%=pp;
    }
    while(!a[a[0]]&&a[0]>1) a[0]--;/****要a[0]大于1   有一组数值ans==0   不大于1 会无输出。。。。*********/ 
}
int gjc(int *a,int x,int b)//把x与2的b次方的乘积赋值给a数组 
{
    a[0]=z[b][0]+4;
    for(int i=1;i<=z[b][0];i++)
    {
        a[i]=z[b][i]*x;
    }
    for(int i=1;i<=a[0];i++)
    {
        a[i+1]+=a[i]/pp;
        a[i]%=pp;
    }
    while(!a[a[0]]&&a[0]>1) a[0]--;
}
int max(int *a,int *b,int *c)//把b数组与c数组中数值大的赋值给a数组 
{
    int pd=0;
    if(b[0]>c[0])   pd=1;
    if(b[0]<c[0])   pd=2;
    if(b[0]==c[0])
    {
        int h=b[0];
        while(h>=1)
        {
            if(b[h]>c[h])
            {
                pd=1;break;
            }
            if(b[h]<c[h])
            {
                pd=2;
                break;
            }
            h--;
        }
    }
    if(pd==1)//DEBUG!!!!!!!!  pd=1!!!!!!!@!$@!$@!
        for(int i=0;i<=b[0];i++)
            a[i]=b[i];
    else
        for(int i=0;i<=c[0];i++)
            a[i]=c[i];
    while(!a[a[0]]&&a[0]>1) a[0]--;
}

int dfs(int x,int l,int r,int s)//记忆化搜索 
{
    if(bo[l][r]) return 0;//记忆化 
    bo[l][r]=1;
    if(l==r)
    {
        gjc(f[l][r],map[x][l],s);/****相当于状态转移方程*******/ 
        return 0;
    }
    memset(gd1,0,sizeof(gd1));/*******************/
    memset(gd2,0,sizeof(gd2));/*初始化很重要**/ 
    memset(gd3,0,sizeof(gd3));/**！！！！！*******/
    memset(gd4,0,sizeof(gd4));/*******************/
    dfs(x,l,r-1,s+1);
    dfs(x,l+1,r,s+1);
    memset(gd1,0,sizeof(gd1));/*******************/
    memset(gd2,0,sizeof(gd2));/*初始化很重要->AC**/ 
    memset(gd3,0,sizeof(gd3));/**！！！！！*******/
    memset(gd4,0,sizeof(gd4));//不能只在上面初始化    因为中间有dfs 和递归过程  期间 gd1 2 3 4 数组的值早已改变！！！！ 
    gjc(gd1,map[x][r],s);           ///***********************/
    gjc(gd2,map[x][l],s);           // ~~~~**********************/
    gjj(gd3,f[l][r-1],gd1);         //~~~~**********************/
    gjj(gd4,f[l+1][r],gd2);         //~~~ *********************/
    max(f[l][r],gd3,gd4);           //~~~此五行 相当于无高精的状态转移方程  
}
int main()
{
    cin>>n>>m;
    init();//预处理 2的n次方
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            cin>>map[i][j];
        }
    }
    for(int i=1;i<=n;i++)
    {   
        memset(f,0,sizeof(f));
        memset(bo,0,sizeof(bo));//对每一层for i的循环  需要初始化 用来记忆的bo数组  以及每层的f数组（滚动） 
        dfs(i,1,m,1);
        memset(gd5,0,sizeof(gd5));
        for(int ii=0;ii<=ans[0];ii++)
            gd5[ii]=ans[ii];
        gjj(ans,gd5,f[1][m]);
    }
    print(ans);
}
olers
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
洛谷 P1005 矩阵取数游戏
我要吃熊猫
05-28 1148
高精度+线性dp
矩阵取数游戏——典型的动态规划问题，值得一看
最新发布
CylMK的博客
07-10 1062
帅帅经常跟同学玩一个矩阵取数游戏：对于一个给定的n×m的矩阵，矩阵中的每个元素aij​均为非负整数。nm×2iii1m帅帅想请你帮忙写一个程序，对于任意矩阵，可以求出取数后的最大得分。
NOIP 2007 矩阵取数游戏 源代码
03-23
NOIP 2007 的题目矩阵取数游戏的pascal源代码，测试数据全部通过。
P1005专属代码
07-25
P1005专属代码，欢迎大家下载，不知道能不能用
P1005 最大公约数.cpp
11-29
P1005 最大公约数.cpp
c#调试程序一例（p1005）
03-20
根据提供的标题、描述以及部分日志内容，我们可以深入解析C#中的程序调试技巧与实践案例。本案例聚焦于一个具体的C#程序调试过程，通过跟踪程序执行的日志记录来理解程序的行为并解决潜在的问题。...
TYVJ题库P1005题 滑雪问题源代码
11-21
描述中提到的“我采用了向四个方向搜索的办法”，暗示了解决方案可能涉及到一种基于网格或矩阵的搜索算法，可能是广度优先搜索（BFS）或深度优先搜索（DFS）。在这种情况下，搜索通常用于遍历可能的路径或者状态空间...
luogu P1005 矩阵取数游戏
zsyz_ZZY的博客
04-03 285
题目传送门：https://www.luogu.org/problemnew/show/P1005题意：有一个n*m的矩阵，求取所有数的最大值。注意：取数有且仅能从当前行两边取且对于每一次取数，都要从每一行取一个，取数的值=被取走的元素值*2^i，i表示当前取到第i个数。思路：区间dp，我们可以知道每一行与其他行没有联系，将每一行的最大值加起来，就可以推方程了。代码：#include&lt;cst...
洛谷OJ: P1005 矩阵取数游戏
308宿舍的鶸
04-16 371
今天因为高数模拟考没怎么做题，就睡前做了这么一题，舍友也开始一起在洛谷刷题了，很开心，但是高数只考了63，很不开心，心情复杂。思路：每一行都是独立的选择，不要因为有多行而被带到错误的思考方向，所以我们开始分析独立的一行数据，那么只有两种选择，1.选择最开头的数字 2.选择最末尾的数字， 而每次增加的分数都是2^i*选择的数字，那么我们稍加思考就可以得到状态转移方程啦。首先，我们定义dp[i][j]...
洛谷P1005 [NOIP2007 提高组] 矩阵取数游戏 题解
q779的博客
07-09 771
区间dp+高精度
P1005 [NOIP2007 提高组] 矩阵取数游戏
qqqingyi的博客
09-17 758
NOIP2007提高组
P1005 矩阵取数游戏（区间dp+java）
Endeavor_G的博客
05-18 542
题目描述 帅帅经常跟同学玩一个矩阵取数游戏：对于一个给定的n \times mn×m的矩阵，矩阵中的每个元素a_{i,j}ai,j​均为非负整数。游戏规则如下： 每次取数时须从每行各取走一个元素，共nn个。经过mm次后取完矩阵内所有元素； 每次取走的各个元素只能是该元素所在行的行首或行尾； 每次取数都有一个得分值，为每行取数的得分之和，每行取数的得分 = 被取走的元素值\times 2^i...
洛谷P1005 矩阵取数游戏
Myriad Dreamin 愿逐将属于我的奇迹
12-30 250
也是从这里开始更会写高精度运算了。 贪心算法和dp很简单，只是加上了高精度而已。#include <iostream> #include <algorithm> #include <cstring> #define mem(_x) memset(_x,0,sizeof(_x)) #define fs(_x,_y) for(int _x=1;i<=_y;_x++) #define fb(_x,_y)
p1005最大公约数c语言
08-25
### 回答1： 下面是一个使用辗转相除法求最大公约数的 C 语言函数： ```c int gcd(int a, int b) { if (b == 0) { return a; } else { return gcd(b, a % b); } } ``` 在此函数中，如果 `b` 为 0，则 `a` 即为最大公约数；否则递归调用 `gcd(b, a % b)`。 ### 回答2： 题目为求解P1005题中的最大公约数，下面是使用C语言编写的解答： ```c #include <stdio.h> // 计算两个数的最大公约数 int gcd(int a, int b) { if (b == 0) { return a; } return gcd(b, a % b); } int main() { int n, m; // 两个数 scanf("%d %d", &n, &m); // 输入两个数 int result = gcd(n, m); // 调用最大公约数函数 printf("%d\n", result); // 输出结果 return 0; } ``` 该程序首先定义了一个`gcd`函数，用于计算两个数的最大公约数。在函数中，使用递归的方式进行求解，如果其中一个数为0，则另一个数即为最大公约数；否则，递归调用求取两个数的余数和较小数的最大公约数。 在`main`函数中，首先定义了两个整数变量`n`和`m`，用于存储输入的两个数。然后使用`scanf`函数从标准输入中读取两个数的值。接下来，调用`gcd`函数求取最大公约数，并将结果保存在`result`变量中。最后，使用`printf`函数将结果输出到标准输出中。 该程序能够求解P1005题中给定的两个数的最大公约数，并输出结果。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值