【2020-8-21】【数字游戏】【启发式搜索IDA*】
有这么一个游戏:
写出一个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*】的更多相关文章
- Python代码实现猜数字游戏
1 # -*- coding:utf-8 -*- 2 # 作者:IT小学生蔡坨坨 3 # 时间:2020/12/9 23:02 4 5 # 猜数字游戏: 6 # 代码中生成一个随机整数. 7 # 然后 ...
- NOIP2003pj数字游戏[环形DP]
题目描述 丁丁最近沉迷于一个数字游戏之中.这个游戏看似简单,但丁丁在研究了许多天之后却发觉原来在简单的规则下想要赢得这个游戏并不那么容易.游戏是这样的,在你面前有一圈整数(一共n个),你要按顺序将其分 ...
- [Noip2003 PJ] 数字游戏
Description & Range 丁丁最近沉迷于一个数字游戏之中.这个游戏看似简单,但丁丁在研究了许多天之后却发觉原来在简单的规则下想要赢得这个游戏并不那么容易.游戏是这样的,在你面前有 ...
- codevs 1085 数字游戏 dp或者暴搜
1085 数字游戏 2003年NOIP全国联赛普及组 时间限制: 1 s 空间限制: 128000 KB 题目描述 Description 丁丁最近沉迷于一个数字游戏之中.这个游戏看似简单 ...
- Poj 2328 Guessing Game(猜数字游戏)
一.题目大意 两个小盆友玩猜数字游戏,一个小盆友心里想着1~10中的一个数字,另一个小盆友猜.如果猜的数字比实际的大,则告诉他"too high",小则"too low& ...
- Codevs 1085 数字游戏
1085 数字游戏 2003年NOIP全国联赛普及组 时间限制: 1 s 空间限制: 128000 KB 题目等级 : 黄金 Gold 题目描述 Description 丁丁最近沉迷于一个数字游戏之中 ...
- 【dp】数字游戏&寒假祭
区间DP 题目描述 丁丁最近沉迷于一个数字游戏之中.这个游戏看似简单,但丁丁在研究了许多天之后却发觉原来在简单的规则下想要赢得这个游戏并不那么容易.游戏是这样的,在你面前有一圈整数(一共n个),你要按 ...
- luogu P1043 数字游戏
题目描述 丁丁最近沉迷于一个数字游戏之中.这个游戏看似简单,但丁丁在研究了许多天之后却发觉原来在简单的规则下想要赢得这个游戏并不那么容易.游戏是这样的,在你面前有一圈整数(一共n个),你要按顺序将其分 ...
- 数字游戏(codevs 1085)
题目描述 Description 丁丁最近沉迷于一个数字游戏之中.这个游戏看似简单,但丁丁在研究了许多天之后却发觉原来在简单的规则下想要赢得这个游戏并不那么容易.游戏是这样的,在你面前有一圈整数(一共 ...
随机推荐
- JVM详解(五)——运行时数据区-方法区
一.概述 1.介绍 <Java虚拟机规范>中明确说明:尽管所有的方法区在逻辑上属于堆的一部分,但一些简单的实现可能不会选择去进行垃圾收集或者进行压缩.但对于HotSpot JVM而言,方法 ...
- NXOpen.UF.UFView.CycleObjects 的使用
Public Sub CycleObjects(ByVal view As NXOpen.Tag, ByVal type As NXOpen.UF.UFView.CycleObjectsEnum, B ...
- spring提供的可拓展接口
接口:SmartLifecycle(https://www.jianshu.com/p/7b8f2a97c8f5)
- Java:容器类线程不安全
Java:容器类线程不安全 本笔记是根据bilibili上 尚硅谷 的课程 Java大厂面试题第二季 而做的笔记 1. Collection 线程不安全的举例 前言 1.当我们执行下面语句的时候,底层 ...
- Spring DeferredResult 异步请求
Spring DeferredResult 异步请求 一.背景 二.分析 三.实现要求 四.后端代码实现 五.运行结果 1.超时操作 2.正常操作 六.DeferredResult运行原理 六.注意事 ...
- PromQL的简单使用
PromQL的简单使用 一.背景 二.PromQL的数据类型 三.字面量 1.字符串字面量 2.浮点数字面量 四.时间序列选择器 1.即时向量选择器 1.组成部分 2.指标名称和匹配器的组合 3.匹配 ...
- 大厂面试题系列:重载(Overload)和重写(Override)的区别。重载的方法能否根据返回类型进行区分
面试题:重载(Overload)和重写(Override)的区别.重载的方法能否根据返回类型进行区分 面试官考察点猜想 这道题纯粹只是考查基础理论知识,对实际开发工作中没有太多的指导意义,毕竟编辑器都 ...
- c语言编程基础入门必备知识
数据类型 基本数据类型 类型名称说明char字符类型存放字符的ASCII码int整型存放有符号整数short短整型存放有符号整数long长整型存放有符号整数long long存放有符号整数float单 ...
- DDR3 IP和CIC IP仿真问题解决记录
1.更新vivado的仿真库(data/secureip和verilog和vhdl文件夹)至最新的vivado库和生成IP的版本匹配: 2.vcs编译脚本里面把仿真库地址指向匹配的仿真库版本: 3.v ...
- cf13B Letter A(分类+简单计算几何,,)
题意: 给三个线段(每个线段的两个端点的坐标),问这三个线段能否组成字母A. 组成字母A的条件: 1.两个线段有公共端点. 2.这两个线段夹角小于等于90度. 3.第三个线段的两个端点分别在这两个线段 ...