有这么一个游戏:
  写出一个1~N的排列a[i],然后每次将相邻两个数相加,构成新的序列,再对新序列进行这样的操作,显然每次构成的序列都比上一次的序列长度少1,直到只剩下一个数字位置。下面是一个例子:
  3  1  2  4
   4  3  6
    7  9
     16
最后得到16这样一个数字。
  现在想要倒着玩这样一个游戏,如果知道N,知道最后得到的数字的大小sum,请你求出最初序列a[i],为1~N的一个排列。若答案有多种可能,则输出字典序最小的那一个。 (n<12)

首先我们通过手算及暴力程序应该可以发现

设 sum 为最后一行的值

总共有 1 行 则 sum=a1*1

总共有 2 行 则 sum=a1*1+a2*1

总共有 3 行 则 sum=a1*1+a2*2+a3*1

总共有 4 行 则 sum=a1*1+a2*3+a3*3+a4*1

不难发现sum等于二项式定理中的系数乘上第n个数

所以对于这题我们只需要先预处理出二项式定理的系数 再对1~n进行dfs排列 适当加上一个剪枝就能AC

#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;
long long a[13],n,m,sum,b[13],last_num[13];
bool k[13],p;
void dfs(long long now,long long num,long long now_num)//当前位置 到当前位置的总和 当前位置的数值
{
if(p)//取到解就不必再继续dfs
return ;
if(now==n)
{
if(num==m)
{
for(long long i=1;i<=n;i++)
printf("%lld ",a[i]);
p=1;
return ;
}
else
{
k[now_num]=0;
return ;
}
} for(long long i=1;i<=n;i++)
{
if(!k[i])
{
if(num+i*b[now+1]>m)//一个小剪枝 如果总和大于所要的解则不再扩展
{
k[now_num]=0;
return;
}
else
{
k[i]=1;
a[now+1]=i;
dfs(now+1,num+i*b[now+1],i);
k[i]=0;//根据函数定义的回溯法求排列
}
}
}
}
int main()
{
scanf("%lld%lld",&n,&m);
for(long long i=1;i<=n;i++)
{
b[i]=1;
for(long long j=i-1;j>=1;j--)
b[j]+=b[j-1];
}
for(long long i=n;i>=1;i--)
{
last_num[i]=b[i]+last_num[i+1];
}
dfs(0,0,0);
}

n<12

对于原题中的n<12 我们的总和超过m就剪枝的做法是完全可以的

但教练提了一个n<20 的问题 这时一个剪枝就显得不足了(只能拿20分)

我们可以用IDA*的算法

即在dfs扩展时计算未来可能的代价 如果代价超出能忍受的范围则剪枝

我们考虑下面三个剪枝

剪枝一:
计算当前的“累加和”,若超出Sum,则没必要继续搜索

剪枝二:
计算当前的“累加和”+“未来的最小累加和”,若超出Sum,则没必要继续搜索

剪枝三:
计算当前的“累加和”+“未来的最大累加和”,若小于Sum,则没必要继续搜索

其中剪枝 二 三 就是对未来的一个估价函数

#pragma GCC optimize(2)
#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;
long long a[21],n,m,sum,b[21],sort_num[21],pai[21],minm,maxm;
bool k[21],p,c;
void yugu(long long now)
{
pai[0]=0;
sort_num[0]=0;
maxm=0;
minm=0;
for(long long i=n/2+1;i>=now;i--)
pai[++pai[0]]=b[i];
for(long long i=n/2+2;i<=n;i++)
pai[++pai[0]]=b[i];
for(long long i=1;i<=n;i++)
{
if(!k[i])
sort_num[++sort_num[0]]=i;
}
for(long long i=1;i<=sort_num[0];i++)
{
minm+=sort_num[i]*pai[i];
maxm+=sort_num[i]*pai[sort_num[0]-i+1];
}
}
void dfs(long long now,long long num)
{ if(c)
return ;
if(num>m)
return ;
yugu(now+1);//进行预估代价
if(num+maxm<m||num+minm>m)//判断是否可行剪枝
return ;
if(now==n)
{
if(num==m)
{
for(long long i=1;i<=n;i++)
printf("%lld ",a[i]);
c=1;
return ;
}
else
return ;
}
for(long long i=1;i<=n;i++)
{ if(!c)
if(!k[i])
{
if(num+i*b[now+1]>m)
continue;
else
{
k[i]=1;
a[now+1]=i;
dfs(now+1,num+i*b[now+1]);
k[i]=0;
}
}
}
}
int main()
{
// freopen("szyx.in","r",stdin);
// freopen("szyx.out","w",stdout);
scanf("%lld%lld",&n,&m);
for(long long i=1;i<=n;i++)
{
b[i]=1;
for(long long j=i-1;j>=1;j--)
b[j]+=b[j-1];
}
dfs(0,0);
return 0;
}

启发式搜索 IDA*

【2020-8-21】【数字游戏】【启发式搜索IDA*】的更多相关文章

  1. Python代码实现猜数字游戏

    1 # -*- coding:utf-8 -*- 2 # 作者:IT小学生蔡坨坨 3 # 时间:2020/12/9 23:02 4 5 # 猜数字游戏: 6 # 代码中生成一个随机整数. 7 # 然后 ...

  2. NOIP2003pj数字游戏[环形DP]

    题目描述 丁丁最近沉迷于一个数字游戏之中.这个游戏看似简单,但丁丁在研究了许多天之后却发觉原来在简单的规则下想要赢得这个游戏并不那么容易.游戏是这样的,在你面前有一圈整数(一共n个),你要按顺序将其分 ...

  3. [Noip2003 PJ] 数字游戏

    Description & Range 丁丁最近沉迷于一个数字游戏之中.这个游戏看似简单,但丁丁在研究了许多天之后却发觉原来在简单的规则下想要赢得这个游戏并不那么容易.游戏是这样的,在你面前有 ...

  4. codevs 1085 数字游戏 dp或者暴搜

    1085 数字游戏 2003年NOIP全国联赛普及组  时间限制: 1 s  空间限制: 128000 KB     题目描述 Description 丁丁最近沉迷于一个数字游戏之中.这个游戏看似简单 ...

  5. Poj 2328 Guessing Game(猜数字游戏)

    一.题目大意 两个小盆友玩猜数字游戏,一个小盆友心里想着1~10中的一个数字,另一个小盆友猜.如果猜的数字比实际的大,则告诉他"too high",小则"too low& ...

  6. Codevs 1085 数字游戏

    1085 数字游戏 2003年NOIP全国联赛普及组 时间限制: 1 s 空间限制: 128000 KB 题目等级 : 黄金 Gold 题目描述 Description 丁丁最近沉迷于一个数字游戏之中 ...

  7. 【dp】数字游戏&寒假祭

    区间DP 题目描述 丁丁最近沉迷于一个数字游戏之中.这个游戏看似简单,但丁丁在研究了许多天之后却发觉原来在简单的规则下想要赢得这个游戏并不那么容易.游戏是这样的,在你面前有一圈整数(一共n个),你要按 ...

  8. luogu P1043 数字游戏

    题目描述 丁丁最近沉迷于一个数字游戏之中.这个游戏看似简单,但丁丁在研究了许多天之后却发觉原来在简单的规则下想要赢得这个游戏并不那么容易.游戏是这样的,在你面前有一圈整数(一共n个),你要按顺序将其分 ...

  9. 数字游戏(codevs 1085)

    题目描述 Description 丁丁最近沉迷于一个数字游戏之中.这个游戏看似简单,但丁丁在研究了许多天之后却发觉原来在简单的规则下想要赢得这个游戏并不那么容易.游戏是这样的,在你面前有一圈整数(一共 ...

随机推荐

  1. 你对微信小程序的理解?优缺点?

    一.是什么 2017年,微信正式推出了小程序,允许外部开发者在微信内部运行自己的代码,开展业务 截至目前,小程序已经成为国内前端的一个重要业务,跟 Web 和手机 App 有着同等的重要性 小程序是一 ...

  2. VS2019中安装2017,2015

    VS2019中安装2017,2015

  3. .jar文件没有Java(TM) Platform SE binary打开方式解决办法

    下面是我个人在打开.jar文件时候的一些小问题: 明明已经配置好了环境变量.jar文件却没有 Java(TM) Platform SE binary 的打开方式, 网上查了资料点明是环境变量的问题,后 ...

  4. selenium3 利用cookie实现免登陆

    1.首先访问要操作的页面 2.登陆一次,使用Fiddle等工具抓取出cookie 3.按照如下代码,即可成功登陆 from selenium import webdriver url = " ...

  5. logstash的安装和简单使用

    logstash的安装和简单使用 一.安装 1.下载并解压 2.logstash 一些命令行参数 1.查看帮助信息 2.加载指定pipeline文件路径 3.检测配置文件语法是否有错误 4.热加载pi ...

  6. CSP2021 翻车记

    DAY - INF 日常模拟赛被吊打,不知道为啥总是出一些小问题导致正解gg,成绩的话也就是中游吧,不过方差不小 DAY - 2 感冒了,头疼得很,签到题甚至也签到失败了,烦得很 DAY -1 全真体 ...

  7. HDI PCB一阶和二阶和三阶如何区分??

      一阶板,一次压合即成,可以想像成最普通的板二阶板,两次压合,以盲埋孔的八层板为例,先做2-7层的板,压好,这时候2-7的通孔埋孔已经做好了,再加1层和8层压上去,打1-8的通孔,做成整板.三阶板就 ...

  8. gawk使用方法简介

    转载:gawk 使用方法简介 - 简书 (jianshu.com) gawk 是最初 Unix 系统上 awk 程序的 GNU 版本.相对于作为流式编辑器的 sed 而言,它提供了更为强大的编程语言特 ...

  9. Linux Ubuntu stty 使用

    stty(set tty)命令用于显示和修改当前注册的终端的属性. 该命令是一个用来改变并打印终端行设置的常用命令. stty -a #将所有选项设置的当前状态写到标准输出中 old_stty_set ...

  10. 理解前端blob和ArrayBuffer,前端接受文件损坏的问题

    1 downloadTemplate().then(res =>{ 2 3 const data = res.data 4 const url = window.URL.createObject ...