背景

贪玩的sunnypig请Charles为他打造一个奇幻世界,Charles欣然答应了。然而一向善于出难题的Charles是决不会轻易让sunnypig轻松拥有一个奇幻世界的,于是Charles在建造过程中设置了重重机关,只有在sunnypig破解了这些障碍之后,才能尝试到奇幻世界中最有玩头的终极宝贝——时空穿梭机。虽然奇幻世界中其他的宝贝也很有趣,但贪玩的sunnypig怎能放过打boss的机会呢?于是他开始了破解障碍的旅程。

描述

第二道障碍来源于一种古老的数学发现——杨辉三角,不过应该是倒过来的杨辉三角。若给出1~n的一个排列A,则将A1、A2相加,A2、A3相加……An-1、An相加,则得到一组n-1个元素的数列B;再将B1、B2相加,B2、B3相加,Bn-2、Bn-1相加,则得到一组n-2个元素的数列……如此往复,最终会得出一个数T。而Charles给sunnypig出的难题便是,给出n和T,再尽可能短的时间内,找到能通过上述操作得到T且字典序最小的1~n的排列。经过汉诺塔问题的训练,sunnypig开始沉着的思考。。。

格式

输入格式

本题有多组数据,对于每组数据:
一行两个整数n(0<n<=20),t即最后求出来的数。

用文件结尾符判断输入结束。

输出格式

对于每组测试数据输出一行n个整数,用空格分开,行尾无多余空格,表示求出来的满足要求的1~n的一个排列。

样例1

样例输入1

4 16
3 9

样例输出1

3 1 2 4
1 3 2

限制

各个测试点2s
不同测试点分数可能不同

建立一个数组b[i][j][k],表示杨辉三角第i行从第j个数第k小的数,比如第4行,1 3 3 1,则b[4][1][1]=1,b[4][1][2]=1;b[4][1][3]=3;b[4][2][1]=1;b[4][2][2]=3;b[4][2][3]=3

搜索剪枝:每次搜索时判断剩余的数和剩余的杨辉数相乘的最小值和最大值,如果最小值加上当前和>T则剪枝,如果最大值加上当前和<T则剪枝

数据太弱了,可以找到大量的算不出来的数据比如 n=20 T=5904462  6143311

 #include<iostream>
#include<string>
using namespace std; int n,t,a[][]={},b[][][];
struct{int x,y;}order[][],tr;
bool over=; void search(int x,int sum,int ans[],bool vis[]){
if(over) return ;
if(x>n)
{
for(int i=;i<n;++i) cout<<ans[i]<<" ";
cout<<ans[n]<<endl;
over=;
return ;
} int min=,max=,d=;
for(int i=;i<=n;++i)
if(vis[i]==)
{max+=i*b[n][x][d];min+=i*b[n][x][n-x-d+];d++;}
if(sum+min>t||sum+max<t) return ;
for(int i=;i<=n;++i)
if(vis[i]==)
{
vis[i]=;
ans[x]=i;
search(x+,sum+a[n][x]*i,ans,vis);
vis[i]=;
} } int main()
{
for(int i=;i<=;++i)
for(int j=;j<=i;++j)
{
if(j==||j==i) a[i][j]=;
else a[i][j]=a[i-][j-]+a[i-][j];
order[i][j].x=a[i][j];
order[i][j].y=j;
} for(int i=;i<=;++i)
for(int j=;j<=i;++j)
for(int k=j+;k<=i;++k)
if(order[i][j].x>order[i][k].x)
{tr=order[i][j];order[i][j]=order[i][k];order[i][k]=tr;}
else if(order[i][j].x==order[i][k].x&&(order[i][j].y>order[i][k].y))
{tr=order[i][j];order[i][j]=order[i][k];order[i][k]=tr;} for(int i=;i<=;++i)
for(int j=;j<=i;++j)
for(int k=;k<=i-j+;++k)
{
int z;int tot=;
for(z=;z<=i&&tot<k;++z)
if(order[i][z].y>=j) tot++;
b[i][j][k]=order[i][z-].x;
} while(cin>>n>>t)
{
over=;
int ans[]={};
bool vis[]={};
search(,,ans,vis); }
// system("pause");
}

最新文章

  1. C# 5.0 异步编程
  2. Python第一模块
  3. webdriver中处理alert
  4. SocketServer model_use
  5. &lt;?php&gt;慢慢写一些php的cookie问题&lt;?&gt;
  6. 【LeetCode】204 - Count Primes
  7. &lt;一道题&gt;求1 + 2! + 3! + .... + N!
  8. heritrix
  9. A+B问题(java)
  10. 转:MFC网络编程学习
  11. Mysql笔记之 -- 开启Mysql慢查询
  12. 计算机网络——DNS协议的学习与实现
  13. ffdshow 源代码分析 3: 位图覆盖滤镜(设置部分Settings)
  14. BZOJ_2662_[BeiJing wc2012]冻结_分层图最短路
  15. shell编程练习(二): 笔试11-20
  16. kubernetes云平台管理实战: 最小的资源pod(二)
  17. IntelliJ Idea 跳出括号并且光标移到末尾的快捷键
  18. [转]Spring MVC之@RequestMapping 详解
  19. Linux账号和密码文件 /etc/passwd和/etc/shadow
  20. python后端工程师 数据爬虫

热门文章

  1. IP V4 和 IP V6 初识
  2. Weblogic 启动慢解决方法
  3. 【WIP】Bootstrap&#160;modal
  4. Road Construction(无向图的双连通分量)
  5. BZOJ 4811 树链剖分+线段树
  6. Centos7基本命令
  7. EasyUI系列学习(三)-Draggable(拖动)
  8. [转]python模块全面
  9. SQL练习题_用户购买收藏记录合并(拼多多)
  10. win32动态库