2017.3.25校内赛（balance）（hill）（draw）（knight）

最新推荐文章于 2021-10-01 22:06:57 发布
最新推荐文章于 2021-10-01 22:06:57 发布
阅读量446 收藏 1
点赞数 2
分类专栏： ————测试题———— 校内赛 ————DP———— 文章标签： 动态规划
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Hawo11/article/details/65937974
版权

动态规划的专练，由远在上海的CKY学长命题。总体而言，题目还是有难度。

天平（balance.in/balance.out）

物理老师YJ有一个长杆天平，天平的两臂长均为15，将长杆看作x轴，则平衡点在0位置处，负数位置在左臂上，正数位置在右臂上。长杆上有n个位置有挂钩可以挂秤砣。YJ有m个秤砣，质量分别为gi，每个挂钩可以不挂也可以挂任意个秤砣。YJ想要知道，在使用所有秤砣的条件下，有多少种不同的挂秤砣的方案，可以使得天平平衡？问题太过复杂，仅凭物理知识难以解决，所以请你来帮助他。
天平的平衡条件是所有秤砣的位置质量之和为0。例如有质量为2,3,4的秤砣分别挂在-3,-2,3位置处，则2(-3) + 3*(-2) + 4*3 == 0，天平是平衡的。
【输入格式】
第一行两个数n，m。表示挂钩的数目和秤砣的数目。
第二行n个不同且递增的数，第i个数表示第i个挂钩的位置，数的范围在[-15,15]内。
第三行m个不同且递增的数，第i个数表示第i个秤砣的质量，数的范围在[1,25]内。
【输出格式】
一个整数，代表能使得天平平衡的方案数。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int M=20;
const int N=15000;
int n,m,tot;
long long c[M+5],map[M+5],f[M+5][N+5];//放了多少东西，背包大小 
int main()
{
    freopen("balance.in","r",stdin);
    freopen("balance.out","w",stdout);
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
    scanf("%I64d",&map[i]);
    for(int i=1;i<=m;i++)
    scanf("%I64d",&c[i]);
    f[0][7500]=1;
    for(int i=1;i<=m;i++)
    for(int k=1;k<=n;k++)//两个相当于一个 
    for(int j=0;j<=N;j++)//背包大小 
    {
        if(j-map[k]*c[i]<0||j-map[k]*c[i]>N) continue;//不超界 
        f[i][j]+=f[i-1][j-map[k]*c[i]];
    }
    printf("%I64d",f[m][7500]);//正好不剩 
    return 0;
}
方法比较简单。

考虑天平的状态可以用秤砣质量距离来表示，而题目情况下，天平的全部秤砣质量距离在(-7500,7500)(20*25*15 = 7500)以内。故可以以此状态来dp。
类似背包问题，把每一个秤砣挂在不同地方，对应成n*m个物品，转移方程：f[i][j]+=f[i-1][j-dis[k]*mg[i]] (i=1..n,j=0..15000,k=1..m)

山峰数（hill.in/hill.out）

山峰数是指数字排列中不存在山谷(先降后升)的数，例如0,5,13,12321都是山峰数，101，1110000111都不是山峰数。
现给出n个数，请依次判断它们是否为山峰数，如果不是，输出-1。如果是，求出比它小的数中有多少个山峰数。
【输入格式】
第一行一个数n，表示询问数目。
接下来n行，每一行一个数x，表示询问的数。
【输出格式】
输出有n行，x如果不是山峰数，输出-1。x如果是山峰数，则输出有多少个比它小的山峰数。

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<string>
#include<cmath>
#include<queue>
#include<algorithm>
#include<iostream>
#define ms(x,y) memset(x,y,sizeof(x))
#define ll long long
using namespace std;

ll f[75][10][2][2];
int len,a[75];
char str[105];

ll dfs(int p,int nu,int isdown,int lim) {
    if(p==len+1)return 1;
    if(f[p][nu][isdown][lim]!=-1)return f[p][nu][isdown][lim];
    ll tmpans=0;
    int r=lim?a[p]:9;
    for(int i=0;i<=r;i++) {
        if(!isdown){
            if(i>=nu) tmpans+=dfs(p+1,i,0,lim && i==r);
            else tmpans+=dfs(p+1,i,1,lim && i==r);
        }
        else if (i<=nu) tmpans+=dfs(p+1,i,1,lim && i==r);
    }
    return f[p][nu][isdown][lim]=tmpans;
}
int main()
{
    freopen("hill.in","r",stdin);
    freopen("hill.out","w",stdout);
    int tc;
    scanf("%d\n",&tc);
    while(tc--){
        scanf("%s",str);
        len=strlen(str);
        for(int i=0;i<len;i++)
            a[i+1]=str[i]-'0';

        bool isdown=false;
        bool ishill=true;
        for(int i=2;i<=len;i++){
            if(a[i]<a[i-1])
                isdown=true;
            if(isdown && a[i]>a[i-1]){
                printf("-1");
                ishill=false;
                break;
            }
        }
        if(ishill) {
            ms(f,-1);
            cout<<dfs(1,0,0,1)-1;
        }
        if(tc!=0)
            printf("\n");
    }
    return 0;
}
数位dp

模板基础上，加入参数nu和isdown，nu表示上一个数填的多少，isdown表示当前是在上升期还是在下降期，然后使用模板记忆化搜索即可。

粉刷匠2（draw.in/draw.out）

有一个4*N的木板需要粉刷，第i行j列的颜色记为A(i, j)。
有256种颜色，记为0..255，为了使得粉刷比较好看，粉刷需要满足如下要求:
1. A(x,y) >= A(x,y-1)
2. 有一些指定的(x1, y1)和(x2, y2)，要求A(x1, y1) = A(x2, y2)
请问有多少种满足要求的粉刷方式？输出答案的最后5位即可。
【输入格式】
第一行两个数n, m，表示木板的长度，和指定规则的条目个数。
接下来m行，每行四个数x1,y1,x2,y1，此规则表示A(x1, y1) 需要等于 A(x2, y2)
。
【输出格式】
一个整数，表示答案的末5位。

#include<cstdio>
#include<cstring>
#define fi first 
#define se second
#define ll long long
#define ms(x,y) memset(x,y,sizeof(x))
using namespace std;


int f[16][16][16][16];
bool vis[16][16][16][16];
int r1x[105],r2x[105],r1y[105],r2y[105];

int main(){ 
    freopen("draw.in","r",stdin);
    freopen("draw.out","w",stdout);
    int n,m,tc;
    int c[4];
    tc=1;
    for(int tt=1;tt<=tc;tt++)
    {
        ms(vis,0);
        scanf("%d%d",&n,&m);
        for(int i=1;i<=m;i++){
            scanf("%d%d%d%d",&r1x[i],&r1y[i],&r2x[i],&r2y[i]);
            r1x[i]--;r2x[i]--;
        }

        for(int i=1;i<=m;i++){
            for(c[0]=0;c[0]<=n;c[0]++)
            for(c[1]=0;c[1]<=n;c[1]++)
            for(c[2]=0;c[2]<=n;c[2]++)
            for(c[3]=0;c[3]<=n;c[3]++){
                if(c[r1x[i]]>=r1y[i] ^ c[r2x[i]]>=r2y[i])
                    vis[c[0]][c[1]][c[2]][c[3]]=1;
            }
        }


        ms(f,0);
        f[0][0][0][0]=1;
        for(int color=0;color<=255;color++){
            for(int col=0;col<=3;col++)
            for(c[0]=0;c[0]<=n;c[0]++)
            for(c[1]=0;c[1]<=n;c[1]++)
            for(c[2]=0;c[2]<=n;c[2]++)
            for(c[3]=0;c[3]<=n;c[3]++)
                if(c[col]<n){
                    int tmp=f[c[0]][c[1]][c[2]][c[3]];
                    c[col]++;
                    f[c[0]][c[1]][c[2]][c[3]] = (f[c[0]][c[1]][c[2]][c[3]] + tmp) % 100000;
                    c[col]--;
                }

            for(c[0]=0;c[0]<=n;c[0]++)
            for(c[1]=0;c[1]<=n;c[1]++)
            for(c[2]=0;c[2]<=n;c[2]++)
            for(c[3]=0;c[3]<=n;c[3]++)
                if(vis[c[0]][c[1]][c[2]][c[3]])
                    f[c[0]][c[1]][c[2]][c[3]]=0;
        }
        printf("%05d\n",f[n][n][n][n]);
    }
    return 0;
}
30分的做法：

dp[i][c1][c2][c3]表示第i行，第一列填到c1，第二列填到c2,，第三列填到c3的方法总数。转移时如果暴力枚举转移仍然会TLE，需要使用完全背包的降低维度思想，列出转移方程：dp[i][c1][c2][c3]= dp[i-1][c1][c2][c3] + dp[i] [c1-1][c2][c3] + dp[i] [c1][c2-1][c3] + dp[i][c1][c2][c3-1]。这样才不会TLE。

100分的做法：

换一种dp思路，我们不按照常规思维枚举行列，而是从小到大枚举颜色。若当前枚举到的颜色为i，我们使用dp[l1][l2][l3][l4]表示每一行使用当前颜色涂到哪一个位置。
假设有一行现在是123344556，现在涂7，显然涂7只能涂序列的后缀，比如涂成123777777，或者123344577才能符合条件。所以考虑枚举当前颜色涂到哪个后缀来转移，转移时同样需要使用完全背包的降维思想优化。
对于限制条件，我们需要把不符合限制条件的去掉，什么样的方案符合限制条件呢？对于限制条件A(x1, y1) = A(x2, y2)，和当前枚举到的颜色i，要么x1行和x2行，当前颜色都涂到了y1,y2位置，要么都没有涂到y1,y2位置。除此之外的都是不合法的方案，我们事先处理好一个vis数组，vis[l1][l2][l3][l4]表示四行分别涂到l1,l2,l3,l4，是否可行，在dp时如果遇到标记不可行的vis，这dp值设为0即可。

棋盘（knight.in/knight.out）

有一个N*M的棋盘，要在上面摆上knight，每个格子可以放至多一个knight。
knight的攻击范围如下图：
所有knight不能互相攻击，请问总共有多少可行的摆法？答案对1000000007取模。
【输入格式】
第一行个数t，表示测试的组数。
接下来t组，每组两个整数，n和m。
【输出格式】
一共t行，第i行表示第i组的答案。

#include<cstdio>
#include<cstring>
#include<iostream>
#define ms(x,y) memset(x,y,sizeof(x))
#define ll long long
using namespace std;

const ll mod=(1e9)+7;
int n,m;
int fi[8]={-2,-1,1,2,2,1,-1,-2};
int fj[8]={1,2,2,1,-1,-2,-2,-1};
int lim;
typedef ll arr[1<<8][1<<8];
arr A,B,s,tmparr;
int p[3][4];
arr storea[5];
arr stores[5];
bool isstore[4];

bool check(int state) {
    ms(p,0);
    int tmp=state;
    for(int i=1;i>=0;i--)
        for(int j=m-1;j>=0;j--,tmp>>=1)
            p[i][j]=tmp&1;
    for(int i=0;i<2;i++)
        for(int j=0;j<m;j++)
            if(p[i][j])
                for(int d=0;d<8;d++) {
                    int ti=i+fi[d],tj=j+fj[d];
                    if(ti>=0 && tj>=0 && ti<=2 && tj<m){
                        if(ti<2 && p[ti][tj]) return false;
                        if(ti==2) p[ti][tj]=1;
                    }
                }
    int line3=0;
    for(int j=0;j<m;j++)
        line3=(line3<<1)+p[2][j];
    int limm=(1<<m)-1;
    line3=(~line3)&limm;   
    for(int i=0;i<=limm;i++)
        if((i|line3)==line3)
            A[state][i+((state&limm)<<m)]=1;
    return true;
}
void Multi(arr A,arr B,arr C) {
    ms(tmparr,0);
    for(int i=0;i<=lim;i++)
        for(int j=0;j<=lim;j++)
            for(int k=0;k<=lim;k++)
                tmparr[i][j]=(tmparr[i][j]+A[i][k]*B[k][j])%mod;
    for(int i=0;i<=lim;i++)
        for(int j=0;j<=lim;j++)
            C[i][j]=tmparr[i][j];
}
int main()
{
    freopen("knight.in","r",stdin);
    freopen("knight.out","w",stdout);
    int tc;
    ms(isstore,0);
    scanf("%d\n",&tc);
    while(tc--){
        scanf("%d%d",&m,&n);
        if(n==1){
            printf("%d\n",1<<m);
            continue;
        }
        ms(A,0);ms(s,0);ms(B,0);
        n-=2;
        lim=(1<<2*m)-1;
        if(!isstore[m]) {
            for(int i=0;i<=lim;i++) {
                if(check(i))s[0][i]=1;
                else s[0][i]=0;
            }
            for(int i=0;i<=lim;i++)
                for(int j=0;j<=lim;j++)
                {
                    stores[m][i][j]=s[i][j];
                    storea[m][i][j]=A[i][j];
                }
            isstore[m]=true;
        }
        else{
            for(int i=0;i<=lim;i++)
                for(int j=0;j<=lim;j++)
                {
                    s[i][j]=stores[m][i][j];
                    A[i][j]=storea[m][i][j];
                }
        }
        for(int i=0;i<=lim;i++)
            B[i][i]=1;

        while(n) {
            if(n&1)Multi(B,A,B);
            Multi(A,A,A);
            n>>=1;
        }
        Multi(s,B,s);
        ll ans=0;
        for(int i=0;i<=lim;i++)
            ans=(ans+s[0][i])%mod;
        cout<<ans<<endl;
    }
    return 0;
}
状态压缩动态规划

每一行使用一个二级制串表示当前行的摆knight情况。用f[i][j][k]表示第i行，前一行为j，前一行的前一行为k时的方法总数，转移时枚举三行的状态，判断合法情况转移即可。
如果枚举每一行转移，可以得70分。如果将状态转移矩阵化，使用矩阵快速幂，可以通过全部测试数据。

以上：2017.3.25

Hawo11
关注
专栏目录
Vivado 2017.2 安装教程（含多版本各类安装包）
09-04 5291
今天给大侠带来FPGA Xilinx Vivado 2017.2 安装教程，话不多说，上货。 各类软件安装包获取方式 Vivado 2017.2版本安装包获取，可在公众号内部回复“Vivado 2017.2安装包”。 这里也为各位大侠附上其他厂商及版本的安装包获取方式，如下： Intel altera Quartus 系列： Quartus版本 获取安装包回复内容 QuartusII9.1 ...
ubuntu16.04安装petalinux 2017.4 教程
Logic+
12-29 6572
前面介绍了vivado 2017.4的安装和启动方法，这篇主要介绍petalinux 2017.4的安装方法， 主要参考资料： UG1144 - PetaLinux Tools Documentation: Reference Guide ( ver2017.4, 2446 KB ) [PDF] 下载petalinux 2017.4： PetaLinux 2017.4 Inst
DP专题（题目总结）
qq_45997091的博客
04-02 194
DP做题总结 一、背包问题（nnn件物品，背包总体积是vvv,如何选择，使得获得的总价值最大？） 背包问题的状态表示都是f[i][j]，只从前i个物品中选，且总体积不超过j的所有选法的集合，每一个f[i][j]都有一个值，这个值就是集合的属性，一般有Max，Min，Count。状态的转换对应的就是集合划分的过程。 1）0-1背包 easy 2）完全背包 与0-1背包的区别就是：0-1背包每个物品只能选择一次，而完全背包可以选择无限次。（前提是不超过背包的体积容量） 以第i个物品选的个数进行划分 f[i,j]
dp专题总结
wyxxzsy的博客
08-19 1132
1
DP专题
lemonoil的博客
03-25 1142
220分，加上long long 成320分。。。无语，果真需要更加细致地分析数据啊！！！天平（balance.in/balance.out） 物理老师YJ有一个长杆天平，天平的两臂长均为15，将长杆看作x轴，则平衡点在0位置处，负数位置在左臂上，正数位置在右臂上。长杆上有n个位置有挂钩可以挂秤砣。YJ有m个秤砣，质量分别为gi，每个挂钩可以不挂也可以挂任意个秤砣。YJ想要知道，在使用所有秤砣的条
dp 专题
meixiuxiudd的博客
08-09 434
一大波的dp题目，算是总结一下 基础dp 1.hdu 1284 背包问题 有1，2，3三种硬币，给你数字求它的兑换方法 硬币可以无限取，背包问题 /******************************************** Author :Crystal Created Time : File Name : *****************
dp专题
薛定谔的Hello Kitty
10-12 189
对“dp[i-1][j-1] 表示替换操作，dp[i-1][j] 表示删除操作，dp[i][j-1] 表示插入操作。”的补充理解： 以 word1 为 “horse”，word2 为 “ros”，且 dp[5][3] 为例，即要将 word1的前 5 个字符转换为 word2的前 3 个字符，也就是将 horse 转换为 ros，因此有： (1) dp[i-1][j-1]，即先将 word1 的前 4 个字符 hors 转换为 word2 的前 2 个字符 ro，然后将第五个字符 word1[4]（因为..
sqlite3.exe绿色版 下载
10-01
SQLite是一款轻型的数据库，它的设计目标是嵌入式的，而且目前已经在很多嵌入式产品中使用了它，它占用资源非常的低，在嵌入式设备中，可能只需要几百K的内存就够了。它能够支持Windows/Linux/Unix等等主流的操作...
drawio.io二次开发采坑记录
maoguan121的博客
10-01 4327
drawio.io二次开发教程 二次开发 基础资料 导入ide 增加本地和oss存储 纯前端功能 各种修改位置介绍 修改代码 准备环境 git clone https://github.com/jgraph/drawio.git cd drawid-dev\src\main\webapp python -m http.server 8081 项目默认是使用打包好的 app.min.js,如果需要修改代码，则需要在请求参数增加 dev=1 同时需要修改 index.html line:245 注
FastReports.NET4-2014.1.3，最好的.net报表控件（破解版）
03-23
FastReport.Net是专为Microsoft Visual Studio 2012和2013或任何其他开发工具功能齐全的报表工程而设计的。使用FastReport.Net，您将能够创建应用程序的独立报告。换句话说，FastReport.Net可以作为一个独立的报告...
dp专题一
oldz
04-24 328
dp专题一 1.数塔 dp[i][j] = max(dp[i-1][j] + dp[i-1][j-1]) + a[i][j] #include &lt;iostream&gt; #include &lt;cstring&gt; #include &lt;cstdio&gt; using namespace std; int a[200][200]; int f[200][20...
Windows PyCharm 2018.3.2 安装图文指南
本篇教程详细介绍了如何在Windows环境下安装PyCharm 2018.3.2，一款流行的Python集成开发环境。首先，用户需确认自己的操作系统位数，无论是32位还是64位，以便下载相应的安装包，可以从JetBrains官方网站获取：...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值