【CEOI2002】【Poj 1038】Bugs Integrated, Inc.

http://poj.org/problem?id=1038

发一下中文题面(今天考试直接被改了):

生记茶餐厅由于受杀人事件的影响,生意日渐冷清,不得不暂时歇业。四喜赋闲在家,整天抱着零食看电视,在大家的提醒下才开始注意自己日益发福的形象,下定决心减肥,萌发了去工作压力大的电脑公司打工的念头。于是,她应聘到了 Bugs 公司,这是一家专门生产硬件的企业。初来乍到,四喜被分配到车间进行产品组装,工作就是把公司生产的一种 2*3 单位尺寸的芯片嵌入 N*M 单位尺寸的模板内。模板接受过严格检查,损坏的单位小方格已被标上黑色记号,如下图所示。

嵌入芯片的要求是,放置芯片的区域内不能有黑色记号,同时芯片与芯片之间不能重叠。公司为了追求利益,希望将尽量多的芯片嵌入模板,这可难坏了那批工人。四喜向你求助,希望你能帮她计算出可能嵌入的最大芯片数量。

①记忆化搜索

悲剧

怎么改都MLE 。。。只能用map动态开内存,但是常数。。。

状态很小,3^10,记录每个点往下还要扩展几格(三格的直接忽略-直接转到下一层了)

// <bugs-b.cpp> - Thu Oct 20 08:13:16 2016
// This file is made by YJinpeng，created by XuYike's black technology automatically.
// Copyright (C) 2016 ChangJun High School, Inc.
// I don't know what this program is.

#include <iostream>
#include <vector>
#include <algorithm>
#include <cstring>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <map>
#pragma GCC push_options
#pragma GCC optimize ("O2")
#define IN inline 
#define RG register 
using namespace std;
const int MAXN=160;
inline int gi() {
    register int w=0,q=0;register char ch=getchar();
    while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
    if(ch=='-')q=1,ch=getchar();
    while(ch>='0'&&ch<='9')w=w*10+ch-'0',ch=getchar();
    return q?-w:w;
}
int n,m,g[MAXN][12],mul[12],a[12];
map<int,int>f[MAXN];map<int,bool>u[MAXN];
IN int dfs(RG int d){
    if(d==n+1)return 0;
    int s=0;
    for(int i=1;i<=m;i++)s+=mul[i]*a[i];
    if(u[d][s])return f[d][s];u[d][s]=1;
    struct kai{
        int d,c,b[12];
        IN void work(int x,int o){
            for(int i=x;i<=m;i++)if(a[i])a[i]--;
            c=max(c,dfs(d+1)+o);//this
            for(int i=x;i<=m;i++)a[i]=b[i];
            if(x>m)return;
            for(int i=x;i<m;i++){
                if(d<n&&a[i]==0&&a[i+1]==0&&!g[d][i]&&!g[d][i+1]&&!g[d+1][i]&&!g[d+1][i+1]){
                    if(d<=n-2&&!g[d+2][i]&&!g[d+2][i+1]){
                        a[i]=a[i+1]=2;//this
                        work(i+2,o+1);
                        a[i]=a[i+1]=0;
                    }
                    if(i<m-1&&a[i+2]==0&&!g[d][i+2]&&!g[d+1][i+2]){
                        a[i]=a[i+1]=a[i+2]=1;//this
                        work(i+3,o+1);
                        a[i]=a[i+1]=a[i+2]=0;
                    }
                }
                if(a[i])a[i]--;
            }for(int i=x;i<=m;i++)a[i]=b[i];
        }
        void pre(int x){
            c=0;d=x;for(int i=1;i<=m;i++)b[i]=a[i];
        }
    }e;e.pre(d);e.work(1,0);
    return f[d][s]=e.c;
}
int main()
{
    int T=gi();
    while(T--){
        for(int i=1;i<MAXN;i++)
            f[i].clear(),u[i].clear();
        memset(g,0,sizeof(g));
        n=gi(),m=gi();int k=gi();
        for(int i=1;i<=k;i++)g[gi()][gi()]=1;
        mul[1]=1;for(int i=2;i<=m;i++)mul[i]=mul[i-1]*4;
        printf("%d
",dfs(1));
    }
    return 0;
}

吐槽:POJ居然卡空间,不得不用map(十倍常数)用数组0.03s就过了今天考试~~表示愤怒，我还是要把我的MLE的代码蒯上.(跑得比状压dp快)

// <bugs-s.cpp> - Thu Oct 20 08:13:16 2016
// This file is made by YJinpeng，created by XuYike's black technology automatically.
// Copyright (C) 2016 ChangJun High School, Inc.
// I don't know what this program is.

#include <iostream>
#include <vector>
#include <algorithm>
#include <cstring>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <map>
#pragma GCC push_options
#pragma GCC optimize ("O2")
#define IN inline 
#define RG register 
using namespace std;
const int MAXN=151;
const int MAXM=60010;
inline int gi() {
    register int w=0,q=0;register char ch=getchar();
    while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
    if(ch=='-')q=1,ch=getchar();
    while(ch>='0'&&ch<='9')w=w*10+ch-'0',ch=getchar();
    return q?-w:w;
}
int n,m,g[MAXN][12],mul[12],a[12];
int f[MAXN][MAXM];bool u[MAXN][MAXM];
IN int dfs(RG int d){
    if(d==n+1)return 0;
    int s=0;
    for(int i=1;i<=m;i++)s+=mul[i]*a[i];
    if(u[d][s])return f[d][s];u[d][s]=1;
    struct kai{
        int d,c,b[12];
        IN void work(int x,int o){
            for(int i=x;i<=m;i++)if(a[i])a[i]--;
            c=max(c,dfs(d+1)+o);//this
            for(int i=x;i<=m;i++)a[i]=b[i];
            if(x>m)return;
            for(int i=x;i<m;i++){
                if(d<n&&a[i]==0&&a[i+1]==0&&!g[d][i]&&!g[d][i+1]&&!g[d+1][i]&&!g[d+1][i+1]){
                    if(d<=n-2&&!g[d+2][i]&&!g[d+2][i+1]){
                        a[i]=a[i+1]=2;//this
                        work(i+2,o+1);
                        a[i]=a[i+1]=0;
                    }
                    if(i<m-1&&a[i+2]==0&&!g[d][i+2]&&!g[d+1][i+2]){
                        a[i]=a[i+1]=a[i+2]=1;//this
                        work(i+3,o+1);
                        a[i]=a[i+1]=a[i+2]=0;
                    }
                }
                if(a[i])a[i]--;
            }for(int i=x;i<=m;i++)a[i]=b[i];
        }
        void pre(int x){
            c=0;d=x;for(int i=1;i<=m;i++)b[i]=a[i];
        }
    }e;e.pre(d);e.work(1,0);
    return f[d][s]=e.c;
}
int main()
{
    freopen("bugs.in","r",stdin);
    freopen("bugs.out","w",stdout);
    int T=gi();
    while(T--){
        memset(f,0,sizeof(f));
        memset(u,0,sizeof(u));
        memset(g,0,sizeof(g));
        n=gi(),m=gi();int k=gi();
        for(int i=1;i<=k;i++)g[gi()][gi()]=1;
        mul[1]=1;for(int i=2;i<=m;i++)mul[i]=mul[i-1]*3;
        printf("%d
",dfs(1));
    }
    return 0;
}

②状压DP

YZC今天给我讲了一下思路,第一行和第二行的状态一定是一样的,关键第三行,开两个数组分别记录第一行和第三行

挖坑待填~(有时间再写)