题目



解析

很神奇的一道题

显然,对于一种排列,相当于给出了数字 \(1..n\) 的对应关系,且不重复不遗漏,刚好把 \(1\) 到 \(n\) 又包含了一遍。

对,连边!

每个数向它对应的数连边,这样我们就得到了一幅图,且这幅图很有特点——全是简单环。

因为每个数的对应有且仅有一个,且不重复不遗漏,所以不存在大环套小环的情况,不存在复杂图形,而且,由于连了边,每个点是什么数字已经没有意义了。

显然,对于一个排列,最大的秩就是所有环大小的最小公倍数。

把 \(n\) 拆成几个正整数的和(正整数可以为 \(1\)),这些正整数的最小公倍数就是我们要的最大的秩,而字典序最小的拆法就是我们要输出的东西。

既然我们已经知道答案就是某个最小公倍数,那我们为什么不直接构造这个最小公倍数?

最小公倍数归根到底是很多质数的乘积,因此我们直接用质数来构造。

假设现在有三个质数 \(p_1、p_2、p_3\),它们的和 \(\leq n\),那么 \(p_1*p_2*p_3\) 一定是一个合法的答案。我们可以先弄一个大小为 \(p_1\) 的环,然后弄一个大小为 \(p_2\) 的环,再弄大小为 \(p_3\) 的环,如果n还有剩余,那么剩下的通通自环,这样一定是合法的。

推广:现在把三个质数变成 \(p1^{c1}、p2^{c2}、p3^{c3}\),他们乘起来依然是合法的。跟上面其实是一样的道理。

再推广:设环的大小为 \(w\),那么 \(w\) 可不可以包含两种或以上素数?可以,但是完全可以转化成上面的情况处理。设 \(w=p1^{c1} * p2^{c2} * p3^{c3}\),它对最小公倍数的贡献就是 \(w\),但是如果我把它拆成三个环来处理,使得每个环大小仅包含一种质数,效果是一样的。而且拆了之后,环的大小也变小了,显然答案的排列会更优。

因为我们的构造方式是:对于一个大小为 \(k\) 的环,我们先输出 \(k - 1\) 个顺序为 \(l+1..l+k-1\) 的数,然后再输 \(l\)。其中 \(l\) 为当前可填的最小的数

问题再次转化

现在问题变成:

现在有一堆质数(\(n\) 以内),每个质数的使用次数上限是已知的(就是上面的 \(c\) 的大小上限)。我要选择一些质数(或它的幂),使他们的和 \(\leq n\),然后这些质数(或它的幂)的乘积要最大。

而这其实就是有限背包问题。

设 \(f_{i,j}\) 表示我们处理到第i个质数、当前和为j所能获得的最大秩。则 \(f_{i,j}=\max(f_{i-1,j-w}*w)\)。而我们要求排列的话,只需记录一下每个状态是由哪个状态转移过来的,最后还原即可。

于是我们似乎做完了

嗯,确实是似乎

因为你会发现一堆质数的最小公倍数会灰常灰常大

所以题解给了一个很牛逼的方法

那就是给 \(f\) 值转为自然对数来做。

即 \(f_{i,j} \Longleftrightarrow \ln{f_{i,j}}\)

于是转移时 \(\times w\) 变为 \(+ \ln{w}\)

\(Code\)

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
using namespace std; const int N = 1e4 + 5;
const double INF = 1e40;
int t , n , vis[N] , tot , pr[1300] , d[1300] , cnt , g[1300][N];
double f[1300][N] , ans; inline void getprime(int Mx)
{
vis[1] = 1;
for(register int i = 2; i <= Mx; i++)
{
if (!vis[i]) pr[++tot] = i;
for(register int j = 1; j <= tot && pr[j] * i <= Mx; j++)
{
vis[pr[j] * i] = 1;
if (i % pr[j] == 0) break;
}
}
} int main()
{
getprime(1e4);
scanf("%d" , &t);
for(; t; t--)
{
scanf("%d" , &n);
memset(f , 0 , sizeof f);
f[0][0] = log(1.0);
for(register int i = 1; i <= tot; i++)
for(register int j = 0; j <= n; j++)
{
f[i][j] = f[i - 1][j];
int w = pr[i];
while (j - w >= 0)
{
if (f[i - 1][j - w] + log(1.0 * w) > f[i][j]) f[i][j] = f[i - 1][j - w] + log(1.0 * w) , g[i][j] = w;
w *= pr[i];
}
}
ans = -INF;
int m , sum = 0;
for(register int i = 0; i <= n; i++)
if (f[tot][i] > ans) ans = f[tot][i] , m = i;
cnt = 0;
for(register int i = tot; i; i--)
{
if (g[i][m]) d[++cnt] = g[i][m] , sum += d[cnt];
m -= g[i][m];
}
while (sum < n) d[++cnt] = 1 , sum++;
sort(d + 1 , d + cnt + 1);
int k = 1;
for(register int i = 1; i <= cnt; i++)
{
for(register int j = 1; j < d[i]; j++)
printf("%d " , k + j);
printf("%d " , k) , k += d[i];
}
printf("\n");
}
}

JZOJ 3232. 【佛山市选2013】排列的更多相关文章

  1. 【佛山市选2013】JZOJ2020年8月7日T4 排列

    [佛山市选2013]JZOJ2020年8月7日T4 排列 题目 描述 一个关于n个元素的排列是指一个从{1, 2, -, n}到{1, 2, -, n}的一一映射的函数.这个排列p的秩是指最小的k,使 ...

  2. 【佛山市选2013】JZOJ2020年8月7日提高组T3 海明距离

    [佛山市选2013]JZOJ2020年8月7日提高组T3 海明距离 题目 描述 对于二进制串a,b,他们之间的海明距离是指两个串异或之后串中1的个数.异或的规则为: 0 XOR 0 = 0 1 XOR ...

  3. 【佛山市选2013】JZOJ2020年8月7日提高组T2 树环转换

    [佛山市选2013]JZOJ2020年8月7日提高组T2 树环转换 题目 描述 给定一棵N个节点的树,去掉这棵树的一条边需要消耗值1,为这个图的两个点加上一条边也需要消耗值1.树的节点编号从1开始.在 ...

  4. 【佛山市选2013】JZOJ2020年8月7日提高组T1 回文子序列

    [佛山市选2013]JZOJ2020年8月7日提高组T1 回文子序列 题目 描述 回文序列是指左右对称的序列.例如1 2 3 2 1是回文序列,但是1 2 3 2 2就不是.我们会给定一个N×M的矩阵 ...

  5. 纪中集训2020.02.05【NOIP提高组】模拟B 组总结反思——【佛山市选2010】组合数计算,生成字符串 PPMM

    目录 JZOJ2290. [佛山市选2010]组合数计算 比赛时 之后 JZOJ2291. [佛山市选2010]生成字符串 比赛时 之后 JZOJ2292. PPMM 比赛时 之后 JZOJ2290. ...

  6. JZOJ 1003 [ 东莞市选 2007 ] 拦截导弹 —— 递推

    题目:https://jzoj.net/senior/#main/show/1003 n^2 的话递推就可以啦. 代码如下: #include<iostream> #include< ...

  7. [JZOJ4024] [佛山市选2015] 石子游戏 解题报告

    Description     Alice 和 Bob 总喜欢聚在一起玩游戏(T_T),今天他(她)们玩的是一款新型的取石子游戏.游戏一开始有N堆石子,Alice 和 Bob 轮流取出石子.在每次操作 ...

  8. 【原创】开源.NET排列组合组件KwCombinatorics使用(二)——排列生成

           本博客所有文章分类的总目录:本博客博文总目录-实时更新 本博客其他.NET开源项目文章目录:[目录]本博客其他.NET开源项目文章目录 KwCombinatorics组件文章目录: 1. ...

  9. 生成n个数的全排列【递归、回溯】

    下面讨论的是n个互不相同的数形成的不同排列的个数.毕竟,假如n个数当中有相同的数,那n!种排列当中肯定会有一些排列是重复的,这样就是一个不一样的问题了. /*===================== ...

  10. xamarin for vs2013

    安装需求(下载的包及版本) 先安装VS2013 然后到官网下Xamarin,运行后会自动下载以下文件 这是下载的详细列表 jdk-6u39-windows-i586.exe(69.73M) Andro ...

随机推荐

  1. redis集群之主从复制集群的原理和部署

    最近在复盘redis的知识,所以本文开始希望介绍下redis的集群架构.原理以及部署:本文主要介绍redis的主从复制集群,包括其架构模型,原理,高可用等: 一.主从集群的介绍   redis的主从复 ...

  2. MySQL进阶实战2,那些年学过的事务

    @ 目录 一.MySQL服务器逻辑架构 二.并发控制 1.读写锁 2.锁粒度 3.表锁 4.行级锁 三.事务 1.原子性(atomicity) 2.一致性(consistency) 3.隔离性(iso ...

  3. USB限流,短路保护芯片IC

    USB口的输出电压一般是5V,在一些电源中,由于总电源5V是一个很大的总电源,再分别出很多路输出负载出来,例如5V10A,分成4个USB输出口,如果没加其他限流和保护的话,任意一个USB口的输出电流都 ...

  4. 基于MATLAB的人民币识别系统

    1.选图背景 随着科技的发展,图像识别技术在计算机视觉领域起着越来越重要的作用,现已被应用到金融.军事.公安刑侦.生物医学.机器视觉等各个领域: 随着国民经济的快速发展,人民币大量流通市场,商业交易都 ...

  5. python2和python3的区别(1)

    1.python2和python3的解释器的默认编码不同 python2解释器的编码默认用的是 ascii python3解释器的编码默认用的 utf-8 2.python2和python3输入的表示 ...

  6. Python全栈工程师之从网页搭建入门到Flask全栈项目实战(6) - Flask表单的实现

    1.表单介绍 1.1.表单知识回顾 常见的表单元素: 表单标签<form> action:表单提交的URL地址 method:表单请求的方式(GET/POSt) enctype:请求内容的 ...

  7. 重学c#系列—— 反射深入一点点[三十三]

    前言 在上一章中介绍了什么是反射: https://www.cnblogs.com/aoximin/p/16440966.html 正文 上一节讲述反射的基本原理和为什么要用反射,还用反射的优缺点这些 ...

  8. [图像处理] YUV图像处理入门3

    5 yuv420格式的灰阶测试图 本程序中的函数主要是为YUV420P视频数据流的第一帧图像添加边框.函数的代码如下所示: /** * @file 5 yuv_graybar.cpp * @autho ...

  9. Python实现单项链表

    单向链表 单向链表也叫单链表,是链表中最简单的一种形式,它的每个节点包含两个域,一个信息域(元素域)和一个链接域.这个链接指向链表中的下一个节点,而最后一个节点的链接域则指向一个空值. 表元素域ele ...

  10. Java环境搭载的低级错误

    关于Java环境搭载的问题 在按照狂神的教程在Dos窗口执行了javac hello.java后在执行java hello报错"" 发现路径文件夹内并未生成class文件 检查了环 ...