【DP】TEST20170520

最新推荐文章于 2022-11-18 20:39:03 发布

Etta19

最新推荐文章于 2022-11-18 20:39:03 发布

阅读量398

点赞数

分类专栏： DP test

本文链接：https://blog.csdn.net/Etta19/article/details/72860301

版权

DP 同时被 2 个专栏收录

20 篇文章 0 订阅

订阅专栏

test

10 篇文章 0 订阅

订阅专栏

举火燎天(lights.pas/c/cpp Time:1s Memory:256M)

问题描述

征夷王将越南(南夷)攻下来之后,决定用 N 个火炬摆成一个圆圈,围住越南的都城,以举火燎天宣告自己的胜利。这 N 个火炬的置已经固定了,但是,每个火炬的颜色都没定,总共有 M 种颜色的火炬。为了显示我们中华文化的博大精深,征夷王决定,相邻的火炬的颜色不能相同。现在征夷王想知道,满足他要求的方案有多少个。由于火炬的位置已经固定,所以即便旋转翻转之后两种方案相同也视作不同方案。

输入

输入文件名为 lights.in。
输入一行两个正整数 N 和 M,意义如上所述。

输出

输出文件名为 lights.out。
输出一行一个正整数,代表方案总数。

输入样例

3 4

输出样例

数据范围

对于 100%的数据,1<=N,M<=100

分析

f[i][0]表示第i个火炬与第一个火炬颜色不同的方案数，f[i][1]表示第i个火炬与第1个火炬颜色相同的方案数
初始化：f[1][0]=0,f[1][1]=m;
递推：f[i][1]=f[i-1][0], f[i][0]=f[i-1][1](m-1)+f[i-1][0](m-2);
ans=f[n][0];

代码

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;

typedef long long LL;

const int MOD = 1e8;

LL n,m;

struct bign {
    LL s[100];
    int len;
    bign() {
        memset(s,0,sizeof(s));
        len=1;
    }

    bign operator = (const LL num) {
        s[1]=num;
        len=1;
        return *this;
    }

    bign operator * (const LL num) {
        bign c;
        c.len=len+1;
        for(int i=1;i<=c.len;++i) {
            c.s[i+1]=(c.s[i]+s[i]*num)/MOD;
            c.s[i]=(c.s[i]+s[i]*num)%MOD;   
        }
        while(c.s[c.len]==0&&c.len>1)--c.len;
        return c;
    }

    bign operator + (const bign & num) {
        bign c;
        c.len=max(num.len,len)+1;
        for(int i=1;i<=c.len;++i) {
            c.s[i+1]=(c.s[i]+s[i]+num.s[i])/MOD;
            c.s[i]=(c.s[i]+s[i]+num.s[i])%MOD;
        }
        while(c.s[c.len]==0&&c.len>1)--c.len;
        return c;
    }

    void print() {
        for(int i=len;i>=1;--i) {
            if(i==len)printf("%lld",s[i]);
            else printf("%08lld",s[i]);
        }
        printf("\n");
    }

}f[100+10][2];

int main() {
    freopen("lights.in","r",stdin);
    freopen("lights.out","w",stdout);
    scanf("%lld%lld",&n,&m);
    if(n==1)printf("%lld\n",m);
    else if(n==2)printf("%lld\n",m*(m-1));
    else {
        f[1][1]=m;f[1][0]=0;
        for(int i=2;i<=n;++i) {
            f[i][1]=f[i-1][0];
            f[i][0]=f[i-1][1]*(m-1)+f[i-1][0]*(m-2);
        }
        f[n][0].print();
    }
    return 0;
}

罗大师的采药(medic.pas/c/cpp Time:0.2s Memory:256M)

问题描述

罗大师掉到了一个山洞,里面有 N 株药材。但是罗大师要教济远小朋友题目,所以他必须赶快离开。但是罗大师需要采药来赚点外快。采某株药需要 Ti 个单位时间,同时这株草药价值 Vi 个金币。罗大师为了考验高二小朋友,所以就提出了 Q 个问题,每个询问是“如果我只有 Mi 个单位时间,我最多能得到多少个金币?”现在,为了高二机房的尊严,请你们回答罗大师的询问。

输入

    输入文件名为 medic.in。
    输入文件第一行为 N 和 Q,代表有 N 株草药和 Q 个询问。
    后接 N 行,每行两个整数,分别为采第 I 株草药的时间 Ti 和其价格 Vi。
    后街 Q 行,每行一个正整数 Mi,代表对于每个询问。

输出

输出文件名为 medic.out。
输出 Q 行,每行一个整数,代表对于罗大师的每个询问你得出的答案。

输入样例

输出样例

数据范围

对于 50%的数据,1≤N,Mi,Q≤1000,;
对于 100%的数据,1≤N,Mi,Q≤100000,Ti,Vi ≤10。

分析

单调队列优化多重背包

代码

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;

const int sm = 10+5, sn = 1e5+5;

int n,q,sum,maxn;
int dq[sn],dqv[sn];
int m[sn],ex[sm][sm];
int w[sn],c[sn],f[sn],ss[sn];

int main() {

    freopen("medic.in","r",stdin);
    freopen("medic.out","w",stdout);

    scanf("%d%d",&n,&q);
    for(int i=1,x,y;i<=n;++i) {
        scanf("%d%d",&x,&y);
        ex[x][y]++;
    }
    for(int i=1;i<=q;++i) {
        scanf("%d",&m[i]);
        maxn=max(maxn,m[i]);
    }

    for(int i=1;i<=10;++i)
        for(int j=1;j<=10;++j) {
            if(!ex[i][j])continue;
            w[++sum]=i;
            c[sum]=j;
            ss[sum]=ex[i][j];
        }

    for(int i=1;i<=sum;++i) 
        for(int j=0;j<w[i];++j) {
            int s=1,t=1;
            for(int k=0;k*w[i]+j<=maxn;++k) {
                while(s<t&&dq[t-1]<=f[j+k*w[i]]-k*c[i])--t;
                dq[t]=f[j+k*w[i]]-k*c[i];
                dqv[t++]=k;
                f[j+k*w[i]]=dq[s]+k*c[i];
                if(k-dqv[s]==ss[i])++s;
            }
        }

    for(int i=1;i<=q;++i)
        printf("%d\n",f[m[i]]);

    return 0;
}

罗大师的棋盘(chessboard.pas/c/cpp Time:0.05s Memory:256M)

问题描述

罗大师有一个 8*8 的棋盘。有一天罗大师有点无聊,于是想将这个切割这个棋盘。每次将原棋盘割下一块矩形棋盘并使剩下部分也是矩形,再将剩下的部分继续如此分割,这样割了(n-1)次后,连同最后剩下的矩形棋盘共有 n 块矩形棋盘。(每次切割都只能沿着棋盘格子的边进行)。原棋盘的每个格子上都有一个权值,剩下的矩形棋盘的总分为其所含各个格子的权值之和。罗大师希望分得的 N 个矩形棋盘的权值相差比较小,也就是均方差最小。罗
大师想知道,最小的可能的均方差为多少?

输入

输入文件名为 chessboard.in。
输入一个数 N,代表罗大师要切割成的块数。后接一个 8*8 的矩阵,表示每个格子上的。

输出

输出文件名为 chessboard.out。
输出一个实数,为最小的可能的均方差。(四舍五入到小数点后 3 位)

输入样例

3
1 1 1 1 1 1 1 3
1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 0
1 1 1 1 1 1 0 3

输出样例

1.633

数据范围

测试点	N的值
1	3
2	4
3	5
4	6
5	7
6	8
7	9
8	10
9	11
10	12

分析

将均方差公式变形，发现平均值是一定的，只需要让每个矩形总分的平方和尽量小，记忆化搜索实现。

代码

#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;

const int inf = 0x3f3f3f3f;

int sum,n;
int a[9][9],f[9][9][9][9][13],val[9][9][9][9];

int dfs(int i,int j,int u,int v,int k) {
    if(f[i][j][u][v][k]!=inf)return f[i][j][u][v][k];
    if(k==1){ 
        return f[i][j][u][v][k]=val[i][j][u][v]*val[i][j][u][v];
    }
    if(i==u&&j==v) return f[i][j][u][v][k]=inf+1;
    if(i<=u&&j<=v&&k>=2) {
        int tmp=inf;
        for(int h=i;h<u;++h) {
            tmp=min(tmp,dfs(i,j,h,v,k-1)+dfs(h+1,j,u,v,1));
            tmp=min(tmp,dfs(i,j,h,v,1)+dfs(h+1,j,u,v,k-1));
        }
        for(int h=j;h<v;++h) {
            tmp=min(tmp,dfs(i,j,u,h,k-1)+dfs(i,h+1,u,v,1));
            tmp=min(tmp,dfs(i,j,u,h,1)+dfs(i,h+1,u,v,k-1));
        }
        if(tmp==inf)++tmp;
        return f[i][j][u][v][k]=tmp;
    }   
    else return f[i][j][u][v][k]=inf+1;
}

double calc(int x) {
    return (double) sqrt(x*1.0/n-(sum*1.0/n)*(sum*1.0/n));
}

int main() {
    freopen("chessboard.in","r",stdin);
    freopen("chessboard.out","w",stdout);
    memset(f,0x3f,sizeof(f));

    scanf("%d",&n);
    for(int i=1;i<=8;++i) 
        for(int j=1;j<=8;++j) {
            scanf("%d",&a[i][j]);
            sum+=a[i][j];   
        }

    for(int i=1;i<=8;++i)
        for(int j=1;j<=8;++j)
            for(int u=i;u<=8;++u)
                for(int v=j;v<=8;++v) {
                    if(i==u&&j==v)val[i][j][u][v]=a[i][j];
                    else val[i][j][u][v]=val[i][j][u][v-1]+val[i][j][u-1][v]-val[i][j][u-1][v-1]+a[u][v];
                }


    dfs(1,1,8,8,n);
    printf("%.3lf\n",calc(f[1][1][8][8][n]));
    return 0;
}