ICode9

精准搜索请尝试: 精确搜索
首页 > 编程语言> 文章详细

第12届蓝桥杯省赛A组C++回路计数(统计哈密尔顿回路个数,状压dp,记忆化搜索,超详解)

2022-01-05 21:37:44  阅读:136  来源: 互联网

标签:状态 12 访问 11111 蓝桥 int state 回路 dp


在这里插入图片描述
答案:881012367360
题意:给一个无向图,求哈密顿回路数
分析:暴力思路就是数搜索路径,时间复杂度是O(n^n)级别的,超过15就跑不动了。所以需要优化。因为是计数问题,不是优化问题,带剪枝的回溯法也不适用。
于是想到了记忆化搜索(类dp),看到21能想到状压dp,正好2的21次int能存下,先尝试使用 dp [2 ^ 21]一维数组,转换成2进制后1代表访问过,0代表没访问过,可以表示所有点是否被访问的状态,dp[state]存有多少种方式可以形成状态state。然后发现状态转移方程无法表示,比如状态11101到状态11111,如果路径为1-3-5-4,则不能和2连通,而路径1-3-4-5就可以接着访问2,只有dp[11101]中的一部分可以为dp[11111]所用,因此dp[11111]不能用dp[11101]推出。
于是可以看出该问题不仅和是否访问有关,还和访问顺序有关,如果只把是否访问作为状态,则状态的粒度过大,无法利用是否联通的信息。我们需要重新设置含访问顺序的状态。但是如果把所有顺序存下,则dp会失效,因为dp的思想本来就是将暴力搜索树中的部分状态合并,如果所有的顺序都存下,和暴力搜索复杂度没有区别。
于是我们继续观察,根据题意互质的结点才相连,在状态10111下路径1-2-3-5和路径1-3-2-5虽然顺序不同,但是最后的5号点都是能和4相连,导出状态11111;还是状态10111,1-3-5-2和路径1-5-3-2都不能和4相连,不能导出状态11111,我们只需要再加一个最末位置,就可以分割开这两种情况,同时把1-2-3-5与1-3-2-5合并为一个状态 dp[10111][5] , 1-3-5-2与1-5-3-2合并为一个状态 dp[10111][2] ,而dp[11111]中的dp[11111][4]=dp[10111][5]+dp[10111][3],于是顺其自然的也推出了状态转移方程,在本例中这个式子的意思是,最后访问4的状态数量等于倒数第2个访问3的路径数加上和倒数第二个访问5的路径数。
总结:
1.看到N<=21左右,想到状态压缩
2.有限制的选择计数问题,考虑(类)dp和记忆化搜索
3.带着实例去分析问题往往更容易发现规律,更高效
4.状态的设置和状态转移方程紧密相连,往往设置状态要从一维开始,如果一维状态无法推出转移方程再考虑二维,高维。
5.记忆化搜索的核心是合并暴力搜索树上的结点,因此一定要朝合并状态方向去分析,写代码时也一定要用到存起来的状态。

#include <bits/stdc++.h>
const int maxn=21;
const int maxs=1<<maxn;
typedef long long ll;
using namespace std;
int N;
ll dp[maxs+100][maxn];//存状态
int g[maxn+5][maxn+5];//存图
ll dfs(int state,int pos){
    if(dp[state][pos]!=-1)return dp[state][pos];//一定要用到存储的状态!!!
    ll res=0;
    for(int i=1;i<N;i++){
        if(((state>>i)&1)==0)continue;
        if(g[i][pos]==1){
            dp[state-(1<<pos)][i]=dfs(state-(1<<pos),i);
            res+=dp[state-(1<<pos)][i];//不要直接用dp[state][pos]加因为初值是-1!!
        }
    }
    return dp[state][pos]=res;
}
int main()
{
    memset(dp,-1,sizeof dp);
    cin>>N;
    while(1);
    for(int i=1;i<=N;i++){
        for(int j=i;j<=N;j++){
            if(__gcd(i,j)==1)g[i-1][j-1]=1,g[j-1][i-1]=1;//1到N映射为0到N-1
        }
    }
    for(int i=1;i<N;i++){//初始化
        dp[(1<<i)+1][i]=1;
    }
    ll sum=0;
    for(int i=1;i<N;i++){
        dp[(1<<N)-1][i]=dfs((1<<N)-1,i);
        sum+=dp[(1<<N)-1][i];
        cout<<dp[(1<<N)-1][i]<<endl;
    }
    cout<<sum<<endl;
    return 0;
}

注意:dp数组要初始化为-1,这样可以判断state有没有更新过,如果更新过可以直接在函数开头return dp[state][pos]
同时记得不要直接dp[state][pos]+=x,因为初值是-1。

附两个有该题解的链接,可以参照着学习:
https://blog.csdn.net/weixin_46239370/article/details/105464885
中的https://blog.csdn.net/weixin_46239370/article/details/116805499
还要
https://blog.csdn.net/qq_36306833/article/details/121872050
这两个博客还有其他题的题解

标签:状态,12,访问,11111,蓝桥,int,state,回路,dp
来源: https://blog.csdn.net/tongjingqi_/article/details/122330386

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有