URAL 1519 Formula 1 【插头DP模板题】

最新推荐文章于 2022-09-25 09:26:08 发布
最新推荐文章于 2022-09-25 09:26:08 发布
阅读量331 收藏
点赞数 1
分类专栏： DP
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Summer_via/article/details/52700759
版权


题意：

题意就是给出一个带障碍的n*m的格子，求有多少条回路恰好经过每个非障碍格子一次？

思路：

插头DP学自 插头DP
写这道题的时候几个地方处理不当，调了许久的bug。

代码：(最小表示法)

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#include<map>
#include<string>
#include<queue>
#include<vector>
#include<list>
#include<bitset>
//#pragma comment(linker,"/STACK:1024000000,1024000000")
using namespace std;
typedef long long ll;
#define INF 0x3f3f3f3f
const int HASH=30007;
const int STATE=1000010;
struct HASHMAP
{
    int head[HASH],next[STATE],size;
    long long state[STATE];
    long long f[STATE];
    void init()
    {
        size=0;
        memset(head,-1,sizeof(head));
    }
    void push(long long st,long long ans)
    {
        int i;
        int h=st%HASH;
        for(i=head[h];i!=-1;i=next[i])//这里要注意是next
          if(state[i]==st)
          {
              f[i]+=ans;
              return;
          }
        state[size]=st;
        f[size]=ans;
        next[size]=head[h];
        head[h]=size++;
    }
}hm[2];
char Mp[20][20];
int n,m;
int ex,ey;
void encode(ll &code,int a[])
{
    code=0;
    for(int i=m; i>=0; i--)
    {
        code<<=3;
        code|=(ll)a[i];
    }
}
void uncode(ll code,int a[])
{
    for(int i=0; i<=m; i++)
    {
        a[i]=code&7;
        code>>=3;
    }
}
void add(int a[],int c)
{
    int val=1;
    for(int i=0; i<c; i++) val=max(val,a[i]+1);
    for(int i=c+1; i<=m; i++) if(a[i]>=val) a[i]++;
    a[c]=a[c+1]=val;
}
void Merge(int a[],int c)
{
    int val=min(a[c],a[c+1]);
    int Max=max(a[c],a[c+1]);
    a[c]=a[c+1]=0;
    for(int i=0; i<=m; i++) {
        if(a[i]==Max) a[i]=val;
            if(a[i]>Max) a[i]--;
    }
}
void show(int a[])
{
    for(int i=0;i<=m;i++) printf("%d%c",a[i],i==m?'\n':' ');
}
void dpblank(int x,int y,int cur)
{
    int a[20];
    ll code;
    HASHMAP &now=hm[cur];
    HASHMAP &to=hm[cur^1];
    for(int i=0;i<now.size;i++)
    {
        uncode(now.state[i],a);
        if(!a[y]&&!a[y+1])
        {
            if(x==n-1||y==m-1||Mp[x+1][y]=='*'||Mp[x][y+1]=='*') continue;
            add(a,y);
            encode(code,a);
            if(y==m-1) code<<=3;
            to.push(code,now.f[i]);
        }
        else if(!a[y]&&a[y+1]||a[y]&&!a[y+1])
        {
            int val=max(a[y],a[y+1]);
            if(x!=n-1&&Mp[x+1][y]=='.')
            {
                a[y]=val,a[y+1]=0;
                encode(code,a);
                if(y==m-1) code<<=3;
                to.push(code,now.f[i]);
            }
            if(y!=m-1&&Mp[x][y+1]=='.')
            {
                a[y]=0,a[y+1]=val;
                encode(code,a);
                to.push(code,now.f[i]);
            }
        }
        else
        {
            if(a[y]==a[y+1]&&!(x==ex&&y==ey)) continue;
            Merge(a,y);
            encode(code,a);
            if(y==m-1) code<<=3;
            to.push(code,now.f[i]);
        }
    }
}
void dpblock(int x,int y,int cur)
{
    int a[20];
    ll code;
    HASHMAP &to=hm[cur^1];
    HASHMAP &now=hm[cur];
    for(int i=0;i<now.size;i++)
    {
        code=now.state[i];
        if(y==m-1) code<<=3;
        to.push(code,now.f[i]);
    }
}
void solve()
{
    int cur=0;
    ll ans=0;
    hm[cur].init();
    hm[cur].push(0,1);
    for(int i=0;i<n;i++)
        for(int j=0;j<m;j++)
    {
        hm[cur^1].init();
        if(Mp[i][j]=='.') dpblank(i,j,cur);
        else dpblock(i,j,cur);
        cur^=1;
    }
    int a[20];
    if(hm[cur].size) ans=hm[cur].f[0];
    printf("%I64d\n",ans);
}
int main()
{
    while(~scanf("%d%d",&n,&m))
    {
        ex=ey=-1;
        getchar();
        for(int i=0; i<n; i++) gets(Mp[i]);
        for(int i=0;i<n;i++)
            for(int j=0;j<m;j++)
        {
            if(Mp[i][j]=='.') ex=i,ey=j;
        }
        if(ex<0){puts("0");continue;}
        solve();
    }
}

代码2：（括号表示法）

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#include<map>
#include<string>
#include<queue>
#include<vector>
#include<list>
#include<bitset>
//#pragma comment(linker,"/STACK:1024000000,1024000000")
using namespace std;
typedef long long ll;
#define INF 0x3f3f3f3f
const int HASH=30007;
const int STATE=1000010;
struct HASHMAP
{
    int head[HASH],next[STATE],size;
    long long state[STATE];
    long long f[STATE];
    void init()
    {
        size=0;
        memset(head,-1,sizeof(head));
    }
    void push(long long st,long long ans)
    {
        int i;
        int h=st%HASH;
        for(i=head[h];i!=-1;i=next[i])//这里要注意是next
          if(state[i]==st)
          {
              f[i]+=ans;
              return;
          }
        state[size]=st;
        f[size]=ans;
        next[size]=head[h];
        head[h]=size++;
    }
}hm[2];
char Mp[20][20];
int n,m;
int ex,ey;
void encode(ll &code,int a[])
{
    code=0;
    for(int i=m; i>=0; i--)
    {
        code<<=2;
        code|=(ll)a[i];
    }
}
void uncode(ll code,int a[])
{
    for(int i=0; i<=m; i++)
    {
        a[i]=code&3;
        code>>=2;
    }
}
void add(int a[],int c)
{
    a[c]=1;
    a[c+1]=2;
}
void Merge(int a[],int c)
{
    int num=0;
    if(a[c]==1&&a[c+1]==1)
    {
        for(int i=c+1;i<=m;i++)
        {
            if(a[i]==1) num++;
            else if(a[i]==2) num--;
            if(num==0&&a[i]==2)
            {
                a[i]=1;break;
            }
        }
    }
    else if(a[c]==2&&a[c+1]==2)
    {
        for(int i=c;i>=0;i--)
        {
            if(a[i]==2) num++;
            else if(a[i]==1) num--;
            if(num==0&&a[i]==1)
            {
                a[i]=2;break;
            }
        }
    }
    a[c]=a[c+1]=0;
}
void show(int a[])
{
    for(int i=0;i<=m;i++) printf("%d%c",a[i],i==m?'\n':' ');
}
void dpblank(int x,int y,int cur)
{
    int a[20];
    ll code;
    HASHMAP &now=hm[cur];
    HASHMAP &to=hm[cur^1];
    for(int i=0;i<now.size;i++)
    {
        uncode(now.state[i],a);
        if(!a[y]&&!a[y+1])
        {
            if(x==n-1||y==m-1||Mp[x+1][y]=='*'||Mp[x][y+1]=='*') continue;
            a[y]=1,a[y+1]=2;
            encode(code,a);
            if(y==m-1) code<<=2;
            to.push(code,now.f[i]);
        }
        else if(!a[y]&&a[y+1]||a[y]&&!a[y+1])
        {
            int val=max(a[y],a[y+1]);
            if(x!=n-1&&Mp[x+1][y]=='.')
            {
                a[y]=val,a[y+1]=0;
                encode(code,a);
                if(y==m-1) code<<=2;
                to.push(code,now.f[i]);
            }
            if(y!=m-1&&Mp[x][y+1]=='.')
            {
                a[y]=0,a[y+1]=val;
                encode(code,a);
                to.push(code,now.f[i]);
            }
        }
        else
        {
            if(a[y]==1&&a[y+1]==2&&!(x==ex&&y==ey)) continue;
//            cout<<x<<' '<<y<<endl;
//            show(a);
            Merge(a,y);
//            show(a);
//            puts("---------------");
            encode(code,a);
            if(y==m-1) code<<=2;
            to.push(code,now.f[i]);
        }
    }
}
void dpblock(int x,int y,int cur)
{
    int a[20];
    ll code;
    HASHMAP &to=hm[cur^1];
    HASHMAP &now=hm[cur];
    for(int i=0;i<now.size;i++)
    {
        code=now.state[i];
        if(y==m-1) code<<=2;
        to.push(code,now.f[i]);
    }
}
void solve()
{
    int cur=0;
    ll ans=0;
    hm[cur].init();
    hm[cur].push(0,1);
    for(int i=0;i<n;i++)
        for(int j=0;j<m;j++)
    {
        hm[cur^1].init();
        if(Mp[i][j]=='.') dpblank(i,j,cur);
        else dpblock(i,j,cur);
        cur^=1;
    }
    int a[20];
    if(hm[cur].size) ans=hm[cur].f[0];
    printf("%I64d\n",ans);
}
int main()
{
    while(~scanf("%d%d",&n,&m))
    {
        ex=ey=-1;
        getchar();
        for(int i=0; i<n; i++) gets(Mp[i]);
        for(int i=0;i<n;i++)
            for(int j=0;j<m;j++)
        {
            if(Mp[i][j]=='.') ex=i,ey=j;
        }
        if(ex<0){puts("0");continue;}
        solve();
    }
}


summer_via
关注
  • 1
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
简单插头DP的模板
11-26
传说中插头DP的模板，参考陈丹琦大牛的论文写的，这只是前面的简单部分的模板，后面的高深部分就木有了。。。
Betsy_Tour一类DP的解法
11-27
例如，Ural1519——Formula1问题提供了类似的框架，但在数据结构和算法设计上更为复杂，因为它涉及到了有障碍物的网格，并要求寻找穿过所有好方格的回路。 #### DP解法详解 在标准的状态压缩DP中，我们通常会遇到...
dp-插头dp
w75779的博客
09-25 962
插头dp 真的烦
插头dp入门题五道（caioj 1489~1493）
fyc的博客
01-16 605
插头dp劣质题解
poj 1739 Tony's Tour 插头dp模板题
sepNINE的专栏
07-23 1279
题意： 给一个迷宫，求左下角到右下角的路径数。 分析： 插头dp的模板题，建议先看cdq的论文再看代码，这份代码在模板基础上略微有改动。论文地址http://wenku.baidu.com/view/ed2b3e23482fb4daa58d4b74.html 代码： #include using namespace std; const int maxD=16; const int H
插头与轮廓线与基于连通性状态压缩的dp 学习指南
Azarath Metrion Zinthos
07-27 1860
什么是插头dp
Ural
10-25
"Ural"是一款字体，可能是指一种特定的中文字体或者西文字体设计。在IT领域，字体扮演着至关重要的角色，特别是在图形设计、网页设计、出版印刷以及各种视觉传达中。字体的选择能够极大地影响信息的呈现效果和用户...
2010FZUACM_高级数据结构习题讲解1
08-04
动态规划的状态转移方程是DP[i] = min(DP[a[j].t1-1], ..., DP[i-1]) + a[j].s，当a[j].t2=i时。线段树或胜者树可以提供O(logn)的查询和更新效率。 5. **并查集**： 在POJ2985 The k-th Largest Group问题中，涉及...
Ural URAL 解题思路
01-11
在编程竞赛和算法挑战的世界中，Ural University的在线判题系统，简称Ural或URAL，是许多程序员磨炼技能的重要平台。它提供了丰富的算法题目，涵盖数据结构、图论、动态规划、数学等多个领域。本篇文章将深入探讨在...
Ural-开源
04-27
这个库的核心特性是利用模板类来表示不同等级的张量，包括标量、向量（等级1）、矩阵（等级2）以及更高阶的张量（如等级3）。这样的设计使得Ural能够处理各种数学对象，同时保持代码的简洁性和可读性。 1. **张量类...
URAL1519 Formula 1（插头dp的基础题 ——详细解释）
qq_36782366的博客
07-28 471
题目可以到vjudge上找来做，原网站注册比较麻烦。 题意： 给你一个n*m的矩阵,矩阵中有一些格子是坏的,现在要你在矩阵中画一个环,这个环要经过矩阵中所有好的格子一次,问你有多少种画法. 分析:典型的插头dp问题,详见cdq的论文>,这里我用的一般解法做的,一般解法速度没有特殊括号表示法快,但是通用性很好. 思路主要注意好状态的转移即可; 首先判断这个这个点是好是坏，坏的话对上
URAL 1519 Formula 1 dp(插头)
便纵有千种风情
08-20 973
题意：给出n*m的方格，有些格子不能铺线，其它格子必须铺，形成一个闭合回路。问有多少种铺法？解:4进制状压表示，用m+1位表示轮廓线上的m个位置的情况。0表示没有插头与线相交，1表示插头’左’,2表示插头’右’。类似括号匹配。重要的一点是轮廓线上方的线路不能相交。转移：1.当前格上左无插头，则该格下右建插头，新的连通分量。2.上或左有插头，必连下或右，仅一个。3.上左均有插头，考虑连通分量的消除。轮
ural 1519 Formula 1(轮廓线|插头DP|括号配对)
GODSPEED
09-09 1164
题意：给你一个m * n的棋盘，有的格子是障碍，问共有多少条回路使得经过每个非障碍格子恰好一次．m, n ≤ 12。 如图，m = n = 4，(1, 1), (1, 2)是障碍，共有2条满足要求的回路． Formula 1" alt="【动态规划】Ural1519 Formula 1" src="http://s4.sinaimg.cn/bmiddle/51cea404x7abee2b588
ural1519之插头DP
Stephen
04-24 2560
1519. Formula 1 Time limit: 1.0 second Memory limit: 64 MB Background Regardless of the fact, that Vologda could not get rights to hold the Winter Olympic games of 20**, it is well-known,
mpu6050陀螺仪使用方法开发教程文档.docx
07-30
mpu6050陀螺仪使用方法
【SCI顶级】龙格库塔算法RUN-CNN-LSTM-Multihead-Attention温度预测【含源码 5747期】.zip
07-30
CSDN海神之光上传的全部代码均可运行，亲测可用，直接替换数据即可，适合小白； 1、代码压缩包内容 主函数：Main .m； 数据； 调用函数：其他m文件；无需运行 运行结果效果图； 2、代码运行版本 Matlab 2019b；若运行有误，根据提示修改；若不会，可私信博主； 3、运行操作步骤 步骤一：将所有文件放到Matlab的当前文件夹中； 步骤二：双击打开除Main.m的其他m文件； 步骤三：点击运行，等程序运行完得到结果； 4、仿真咨询 如需其他服务，可私信博主或扫描博主博客文章底部QQ名片； 4.1 CSDN博客或资源的完整代码提供 4.2 期刊或参考文献复现 4.3 Matlab程序定制 4.4 科研合作 智能优化算法-CNN-LSTM-Mutilhead-Attention回归预测系列程序定制或科研合作方向： 4.4.1 遗传算法GA/蚁群算法ACO-CNN-LSTM-Mutilhead-Attention回归预测 4.4.2 粒子群算法PSO/蛙跳算法SFLA-CNN-LSTM-Mutilhead-Attention回归预测 4.4.3 灰狼算法GWO/狼群算法WPA-CNN-LSTM-Mutilhead-Attention回归预测 4.4.4 鲸鱼算法WOA/麻雀算法SSA-CNN-LSTM-Mutilhead-Attention回归预测 4.4.5 萤火虫算法FA/差分算法DE-CNN-LSTM-Mutilhead-Attention回归预测
大模型的数学之魂：统计学习理论精要
07-30
大型机器学习模型，通常简称为“大模型”，是指具有数十亿甚至数千亿参数的深度神经网络模型。这些模型因为其庞大的参数量和复杂的计算结构，能够处理非常复杂的任务和数据，展现出强大的表达能力和预测性能。大模型在自然语言处理、计算机视觉、语音识别和推荐系统等多个领域都有广泛的应用。 大模型的设计目的是通过训练海量数据来学习复杂的模式和特征，从而具备强大的泛化能力，能够对未见过的数据做出准确的预测。与参数较少、层数较浅的小模型相比，大模型虽然需要更多的计算资源和时间来训练和推理，但它们展现出了一种被称为“涌现能力”的特性，即从原始训练数据中自动学习并发现新的、更高层次的特征和模式。 大模型的发展历程可以分为几个阶段，从早期的卷积神经网络模型发展到以Transformer为代表的全新神经网络模型，再到当前以GPT为代表的预训练大模型阶段。Transformer架构因其自注意力机制而成为大模型的核心技术，它使得模型能够更好地理解和处理文本数据中的长距离依赖关系。 大模型的分类主要包括语言大模型（NLP）、视觉大模型（CV）和多模态大模型。语言大模型专注于处理文本数据和理解自然语言，而视觉大模型用
唯美淡雅四页.zip
最新发布
07-30
【简历模板】工作总结、商业计划书、述职报告、读书分享、家长会、主题班会、端午节、期末、夏至、中国风、卡通、小清新、岗位竞聘、公司介绍、读书分享、安全教育、文明礼仪、儿童故事、绘本、防溺水、夏季安全、科技风、商务、炫酷、企业培训、自我介绍、产品介绍、师德师风、班主任培训、神话故事、巴黎奥运会、世界献血者日、防范非法集资、3D快闪、毛玻璃、人工智能等等各种样式的ppt素材风格。 设计模板、图片素材、PPT模板、视频素材、办公文档、小报模板、表格模板、音效配乐、字体库。 广告设计：海报，易拉宝，展板，宣传单，宣传栏，画册，邀请函，优惠券，贺卡，文化墙，标语，制度，名片，舞台背景，广告牌，证书，明信片，菜单，折页，封面，节目单，门头，美陈，拱门，展架等。 电商设计：主图，直通车，详情页，PC端首页，移动端首页，钻展，优惠券，促销标签，店招，店铺公告等。 图片素材：PNG素材，背景素材，矢量素材，插画，元素，艺术字，UI设计等。 视频素材：AE模板，会声会影，PR模板，视频背景，实拍短片，音效配乐。 办公文档：工作汇报，毕业答辩，企业介绍，总结计划，教学课件，求职简历等PPT/WORD模板。
仿真验证阵列可以提高接收信号信噪比
07-30
这是一份模拟了阵列输入信号及噪声，并验证了相对于阵列接收到的信号，阵列输出信号可以将信噪比提高M倍，其中M为阵列的阵元个数的代码。 代码中可以随意修改阵元个数、阵元间距、波束指向角度、信号频率等。 代码中关键部分均含有文字注释，完全不必担心看不懂。 无论是从仿真波形，还是计算的信噪比结果均能看出阵元数为M的阵列将信号的信噪比提高了M倍。
NOIP进阶：插头DP策略与连通性状态压缩解析
例如，课件中给出的Formula1(Ural1519)问题是一个m*n的棋盘，目标是找到经过所有非障碍格子的哈密顿回路的个数，其中m和n的值不超过12。由于数据规模较小，动态规划是解决此类问题的理想选择，避免了暴力搜索可能...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值