Description

AK:All kill
“你为什么没背书?”
“没有为什么,我就是没背书。”
“……我去年买了个表,G—U—N!”
头铁王InFleaKing把背书的时间都拿去列排列了......
n=3的排列一共有六个(顺序按字典序从小到大):
1 2 3
1 3 2
2 1 3
2 3 1
3 1 2
3 2 1
气不打一处来的InFleaKing把n的排列打乱了。
他想知道字典序第k小的n的排列是什么?
由于InFleaKing被捉去背书了,所以这个问题只能交给被万人顶礼膜拜的dalao您来解决了。

Input

一行,两个数,分别是n,k。
n,k的含义见题目描述。

Output

一行,n个数。
代表字典序第k小的n的排列。
注意两两数之间要用空格隔开。

Sample Input

Sample Input1:

1 1

Sample Input2:

2 2

Sample Output

Sample Output1:

1

Sample Output2:

2 1

Data Constraint

【数据约定】
对于10%的数据:1<=n<=3
对于20%的数据,1<=n<=9
对于30%的数据:1<=n<=18,1<=k<=10^6
对于60%的数据:1<=n<=18
对于80%的数据:1<=n<=100,1<=k<=10^100
对于90%的数据:1<=n<=1000,1<=k<=10^1000
对于100%的数据:1<=n<=100000,1<=k<=min(10^20000,n!)

思路解析

无力吐槽。
如果想恶心我的话成功了。

根据惯例,我们面向数据编程。

第一档

对于10%的数据:1<=n<=3
手算 + 打表
预期分数:10

第二档

对于30%的数据:1<=n<=18,1<=k<=10^6
手算 + 打表还是可以
直接暴力递归求解。
预期分数:20 ~ 30

第三档

对于60%的数据:1<=n<=18
要不要听一个看上去很高级的玩意叫康托展开?
这个题听说是逆康托展开。没兴趣讲,请自行baidu

说一下我的想法
因为是输出字典序,我姑且列出所有情况看一下吧
以4 6这组数据为例

1 2 3 4
1 2 4 3
1 3 2 4
1 3 4 2
1 4 2 3
1 4 3 2
2 1 3 4
2 1 4 3
2 3 1 4
2 3 4 1
2 4 1 3
2 4 3 1
........

可以找出一个规律:
第一位数是当前1~4中没用过的第k/6大的数
第二位是当前1~4中没用过的第k/2大的数

6 = 3!
2 = 2!
好了,写代码吧

#include<iostream>
#include<cstdio>
#include<cstdlib>
using namespace std;
const long long sum[] = {,,,,,,,,,,,,,,,,,,};
long long n,k;
bool vis[];
inline void find_kth(long long x) {
long long res = ;
long long i;
for(i = ;i <= n;i++) {
if(!vis[i]) res++;
if(res == x) break;
}
printf("%lld ",i);
vis[i] = ;
}
int main() {
freopen("array.in","r",stdin);
freopen("array.out","w",stdout);
scanf("%lld%lld",&n,&k);
long long now,tmp,cnt,x;
now = k;
cnt = n;
while(cnt--) {
tmp = now / sum[cnt];
x = now;
now = now % sum[cnt];
if(!now) now = sum[cnt];
if(tmp * sum[cnt] < x) tmp++;
find_kth(tmp);
}
return ;
}
第四档

对于90%的数据:1<=n<=1000,1<=k<=10^1000
高精不解释
每次的商是小于n的,所以可以二分商涉及到高精度减法,高精度乘法,还算好打。
一次的复杂度是O(n log n)
总复杂度O(n^2 log n)
预期分数:70~90

第五档

对于100%的数据:1<=n<=100000,1<=k<=min(10^20000,n!)
以下引自作者的题解》》

k<=10^20000
我可以很负责任的告诉你:6100!>10^20000,但是我们的n却有10^5那么大。
这说明什么?
这说明前100000-6100=93900个数一定是1、2、3、…、93900(你发现了
吗?)
因此只需要O(6100^2)的复杂度来完成这个部分。
即使你不会数据结构或者分块也可以解决这道题!

假数据范围
预计分数:90~100(需要常数优化)

总结

考场不要写高精,除非你时间充裕

 

[ZJOJ] 5772【NOIP2008模拟】今天你AK了吗的更多相关文章

  1. JZOJ 5777. 【NOIP2008模拟】小x玩游戏

    5777. [NOIP2008模拟]小x玩游戏 (File IO): input:game.in output:game.out Time Limits: 1000 ms  Memory Limits ...

  2. JZOJ 5809. 【NOIP2008模拟】数羊

    5809. [NOIP2008模拟]数羊 (File IO): input:sheep.in output:sheep.out Time Limits: 1000 ms  Memory Limits: ...

  3. JZOJ 5793. 【NOIP2008模拟】小S练跑步

    5793. [NOIP2008模拟]小S练跑步 (File IO): input:run.in output:run.out Time Limits: 2000 ms  Memory Limits:  ...

  4. JZOJ 5791. 【NOIP2008模拟】阶乘

    5791. [NOIP2008模拟]阶乘 (File IO): input:factorial.in output:factorial.out Time Limits: 1000 ms  Memory ...

  5. JZOJ 5776. 【NOIP2008模拟】小x游世界树

    5776. [NOIP2008模拟]小x游世界树 (File IO): input:yggdrasil.in output:yggdrasil.out Time Limits: 1500 ms  Me ...

  6. JZOJ 5775. 【NOIP2008模拟】农夫约的假期

    5775. [NOIP2008模拟]农夫约的假期 (File IO): input:shuru.in output:shuru.out Time Limits: 1000 ms  Memory Lim ...

  7. JZOJ 5773. 【NOIP2008模拟】简单数学题

    5773. [NOIP2008模拟]简单数学题 (File IO): input:math.in output:math.out Time Limits: 1000 ms  Memory Limits ...

  8. JZOJ 5771. 【NOIP2008模拟】遨游

    5771. [NOIP2008模拟]遨游 (File IO): input:trip.in output:trip.out Time Limits: 2000 ms  Memory Limits: 2 ...

  9. JZOJ5776. 【NOIP2008模拟】小x游世界树

    题目:[NOIP2008模拟]小x游世界树: 题目的附加题解给的很清楚,这里只给一个代码: #include<iostream> #include<cstdio> #inclu ...

随机推荐

  1. Hibernate基于注解的双向one-to-many映射关系的实现

    在项目中用到了一对多的实体类关系映射,之前接触的都是基于配置文件的映射实现.可是公司的大部分都是基于注解的.因此自己參考之前的代码捣鼓了基于注解的一对多的映射关系实现. 背景: 一的一端:QingAo ...

  2. CSDN博客怎样加入量子恒道统计?

    CSDN博客申请专家通过后,能够加入量子恒道统计,这样就能够查看更加具体的訪问统计信息,而不不过一个訪问次数.达到专家级别了可能都是电脑高手了.这里班门弄斧了. 登录博客进入个人中心首页,点击管理博客 ...

  3. charts柱状图,定时刷新

    option: var dataoptionone = { title : { text: '数据存储情况', subtext: '数据表', x:'center' }, tooltip: { sho ...

  4. MBLAST - BLAST

    There are given two strings, A and B. An expansion of some string X is a string created by adding or ...

  5. 2-1 创建第一个Vue实例

    https://cn.vuejs.org/ 生产版本是经过压缩后的版本,它里面把一些没用的警告的代码删除掉了所以它会更精悍一些. 直接引入CDN的网址我们就可以使用Vuejs了.使用生产版本或者CDN ...

  6. hibernate字段名和属性

    字段名和属性名相同 Annotation:默认为@Basic 注意:如果在成员属性没有加入任何注解,则默认在前面加入了@Basic Xml中不用写column 字段名和属性名不同 Annotation ...

  7. PCB InCAM 获取 JOB STEP 实现外挂脚本调试功能实现

    PCB CAM自动化基于Incam 打造,在测试时经常遇到调试障碍,每次自行对功能测试时,生成了exe脚本后,再到Incam里面运行,发现问题,再回来修改代码,非常不爽, 参考Genesis调试运行模 ...

  8. 2.android

    ImageButton action_btn = (ImageButton) findViewById(R.id.action_btn);action_btn.setOnClickListener(n ...

  9. codevs3728联合权值(LCA)

    3728 联合权值  时间限制: 1 s  空间限制: 128000 KB  题目等级 : 黄金 Gold 题解  查看运行结果     题目描述 Description 输入描述 Input Des ...

  10. P3194 [HNOI2008]水平可见直线

    传送门 我们把所有的直线按斜率从小到大排序,然后用单调栈维护 发现,如果当前直线与\(st[top-1]\)直线的交点的横坐标大于等于与\(st[top]\)的交点的横坐标,当前直线可以覆盖掉\(st ...