NOIP模拟测试2「排列 (搜索)·APIO划艇」
排序
题目描述
输入格式
数据范围与提示
对于30%的数据,1<=N<=4; 对于全部的数据,1<=N<=12。
一群数论题中出了一个搜索。
甚至我考试时一点也不会。
自己算复杂度比较大然后交上去只跑了90ms。
这也太
不多说了。
引理1
当前可以将原数列变成排序的操作,其每一种其他排列都可以将原序列变成排序。
证明 不太好证,举几个例子
例如 3 4 2 1这个序列可以经过操作1 操作2 变成 1 2 3 4
那么显然 先进行操作2 再进行操作1 也可以变成 1 2 3 4
于是推而广之 (数学归纳)用cnt次操作可以达到目的,那么所有符合的方案数为cnt!
这样通过引理1 可以剪很多枝
我们从小到大枚举每一种方案,每次验证符合,若符合ans+cnt!就完了
引理2
假设枚举到第i个方案 设划分序列为区间长度为1<<(i+1)的序列,如果有三个及以上不符合递增且相连的区间,那么这个序列没救了。
这个我们可以验证
首先我们保证如果我们进行到第i个方案那么所有长度1<<i的序列都符合递增且相连
如果有一个长度1<<(i+1)不是递增且相连的区间,我们把它划分成两个1<<i,i<<i的序列,我们直接用第i种操作(可能)就行了
然后如果有两个长度那么把它划分成四个,然后也(可能)行
但如果有三个及以上我们最多保证经过i操作后,两个序列符合条件
所以得证
然后搜就完了
以下是本人丑陋的代码
#include<bits/stdc++.h>
#define ll long long
#define A 1000000
using namespace std;
ll xx[A],a[A],b[A],ans,n;
ll read()
{
ll x=0,f=1;char ch=getchar();
while(!isdigit(ch)){
if(ch=='-') f=-1;
ch=getchar();
}
while(isdigit(ch)){
x=x*10+ch-'0';
ch=getchar();
}
return f*x;
}
void SWAP(ll x,ll k,ll len)
{
for(ll i=x,j=k;len;len--,i++,j++)
swap(a[i],a[j]);
}
bool check(ll l,ll r)
{
for(ll i=l;i<r;i++)
if(a[i]+1!=a[i+1])
return 0;
return 1;
}
void shuchu()
{
for(ll i=1;i<=(1<<n);i++)
cout<<a[i]<<"";
cout<<endl;
} void dfs(ll k,ll cnt)
{
// shuchu();
// cout<<k<<endl;
if(k==n)
{
ans+=xx[cnt];
return ;
}
ll now=(1<<(k+1)),x=(1<<k),b[5],tot=0;
for(ll i=1;i<=(1<<n);i+=now)
{
// printf("检查i=%lld j=%lld\n",i,i+now-1);
if(!check(i,i+now-1))
{
b[++tot]=i;
if(tot==3) return ;
}
}
// cout<<"tot="<<tot<<"x="<<x<<"now="<<now<<endl;
// for(ll i=1;i<=n;i++)
// {
//
// }
if(tot==0) {dfs(k+1,cnt);return ;}
if(tot==1)
{
SWAP(b[1],b[1]+x,x);
if(check(b[1],b[1]+x)) dfs(k+1,cnt+1);
SWAP(b[1],b[1]+x,x);
}
if(tot==2)
{
SWAP(b[1],b[2],x);
if(check(b[1],b[1]+x)&&check(b[2],b[2]+x)) dfs(k+1,cnt+1);
SWAP(b[1],b[2],x);
SWAP(b[1]+x,b[2],x);
if(check(b[1],b[1]+x)&&check(b[2],b[2]+x)) dfs(k+1,cnt+1);
SWAP(b[1]+x,b[2],x);
SWAP(b[1],b[2]+x,x);
if(check(b[1],b[1]+x)&&check(b[2],b[2]+x)) dfs(k+1,cnt+1);
SWAP(b[1],b[2]+x,x);
SWAP(b[1]+x,b[2]+x,x);
if(check(b[1],b[1]+x)&&check(b[2],b[2]+x)) dfs(k+1,cnt+1);
SWAP(b[1]+x,b[2]+x,x);
}
// shuchu();
}
int main()
{
n=read();xx[0]=1;
for(ll i=1;i<=(1<<n);i++) a[i]=read();
for(ll i=1;i<=n;i++) xx[i]=xx[i-1]*i;
dfs(0,0);
cout<<ans<<endl;
}
划艇
有点难。
至少对于我来说有点难。弱的我在loj交了5页wa 才AC。
不多说了。
https://blog.csdn.net/qq_22541499/article/details/51674707
https://www.cnblogs.com/hzoi-DeepinC/p/11158757.html#4298030
https://www.cnblogs.com/heoitys/p/11166727.html
这三篇博客已经讲的很清楚了。
所以我鸽掉了
必须写出来,填掉这个这个坑,要不这个题就白做了。
首先这个题必须离散化,肯定开不了这么大的数组。
我们把区间离散化一下,把每段区间设置成左开右闭(或者左闭右开)的,把每个区间出现次序处理出来,左右端点处理出来,区间长度处理出来,
注意这里区间不可以设置成左闭右闭的,因为如果这样设置会出现重复枚举,前面这段区间会枚举最右面这个端点,后面这个区间会枚举最左面这个区间。
然后思考这个题肯定要dp
说明f数组含义。
设$f[i][j]$表示第i个学校选择第j个区间中某一个数
思考会被什么转移
首先如果前面所有学校都没有选择j之后的区间(要么没选要么选择了j之前的区间),那么这个区间可以选择任何一个数。
$\sum\limits_{l=1}^{i-1}\sum\limits_{r=1}^{j-1} f[l][r]*C_{len}^{1}$
然后如果前面有学校同样选择了j这个区间
先假设前面选择这个区间的学校为k
那么对当前这个区间的贡献即为
$\sum\limits_{l=1}^{k-1}\sum\limits_{r=1}^{j-1} f[l][r]$*(k~i之间一堆方案数)
$\sum\limits_{l=1}^{k-1}\sum\limits_{r=1}^{j-1} f[l][r]$$\sum\limits_{l=1}^{i-1}\sum\limits_{r=1}^{j-1} f[l][r]*C_{len}^{1}$这两个式子都可以用前缀和进行维护
那么我们现在要做的就是求出来(k~i之间一堆方案数)
剩下的明天写
NOIP模拟测试2「排列 (搜索)·APIO划艇」的更多相关文章
- NOIP模拟测试21「折纸·不等式」
折纸 题解 考试时无限接近正解,然而最终也只是接近而已了 考虑模拟会爆炸,拿手折纸条试一试,很简单 考你动手能力 代码 #include<bits/stdc++.h> using name ...
- NOIP模拟测试10「大佬·辣鸡·模板」
大佬 显然假期望 我奇思妙想出了一个式子$f[i]=f[i-1]+\sum\limits_{j=1}^{j<=m} C_{k \times j}^{k}\times w[j]$ 然后一想不对得容 ...
- NOIP模拟测试19「count·dinner·chess」
反思: 我考得最炸的一次 怎么说呢?简单的两个题0分,稍难(我还不敢说难,肯定又有人喷我)42分 前10分钟看T1,不会,觉得不可做,完全不可做,把它跳了 最后10分钟看T1,发现一个有点用的性质,仍 ...
- NOIP模拟测试30「return·one·magic」
magic 题解 首先原式指数肯定会爆$long$ $long$ 首先根据欧拉定理我们可以将原式换成$N^{\sum\limits_{i=1}^{i<=N} [gcd(i,N)==1] C_{G ...
- NOIP模拟测试24「star way to hevaen·lost my music」
star way to heaven 题解 大致尝试了一下并查集,记忆化搜索,最小生成树 最小生成树是正解,跑最小生成树然后找到最大的值 欧几里德距离最小生成树学习 prim楞跑 至于为什么跑最小生成 ...
- NOIP模拟测试18「引子·可爱宝贝精灵·相互再归的鹅妈妈」
待补 引子 题解 大模拟,注意细节 代码1 #include<bits/stdc++.h> using namespace std; int n,m;char a[1005][1005]; ...
- 20190803 NOIP模拟测试12「斐波那契(fibonacci)· 数颜色 · 分组 」
164分 rank11/64 这次考的不算太差,但是并没有多大的可能性反超(只比一小部分人高十几分而已),时间分配还是不均,T2两个半小时,T1半个小时,T3-额十几分钟吧 然额付出总是与回报成反比的 ...
- NOIP模拟测试8「寿司」
考试时打的类似$n^2$暴力,然后炸了只有10分 后来验证我的算法伪了. 题解 显然你有一种解法,假设你要在一个B点断开将R分别移向最左 最右,这样只用分别计算B点右面蓝色数量左面蓝色数量就得到了一个 ...
- NOIP模拟测试38「金·斯诺·赤」
金 辗转相减见祖宗 高精 #include<bits/stdc++.h> using namespace std; #define A 2000 #define P 1 #define N ...
随机推荐
- JSON简单了解
JSON简单了解 简介 JSON (JavaScript Object Notation):一种简单的数据格式,比xml更轻巧.JSON 是 JavaScript 原生格式,这意味着在 JavaScr ...
- .NET Core 服务诊断工具
前言: 前一篇文中介绍了.NET Core-全局性能诊断工具 的使用方法,那么接下来自己实现一个简单.NET Core的诊断工具. 该工具主要包括:.NET Core 程序进程信息查看.性能计数器结果 ...
- MySQL如何快速插入数据
前言: 日常学习和工作中,经常会遇到导数据的需求.比如数据迁移.数据恢复.新建从库等,这些操作可能都会涉及大量数据的导入.有时候导入进度慢,电脑风扇狂转真的很让人崩溃,其实有些小技巧是可以让导入更快速 ...
- (转)如何优雅的使用rabbit mq
RabbitMQ无疑是目前最流行的消息队列之一,对各种语言环境的支持也很丰富,作为一个.NET developer有必要学习和了解这一工具.消息队列的使用场景大概有3种: 1.系统集成,分布式系统的设 ...
- Linux的三剑客
首先,需要介绍一下管道和正则表达式,因为它经常和Linux三剑客一起使用. 一.管道Linux 提供管道符"|",将两个命令隔开,管道符左边命令的输出作为管道符右边命令的输入. c ...
- FFmpeg应用实践之命令查询
0. 前言 FFmpeg 中常用的工具有三个,分别是多媒体编解码工具ffmpeg.多媒体内容分析工具ffprobe和多媒体播放器ffplay.本文介绍的指令都是与编解码工具 ffmpeg 相关的. 学 ...
- .NetCore·集成Ocelot组件之完全解决方案
阅文时长 | 11.04分钟 字数统计 | 17672.8字符 主要内容 | 1.前言.环境说明.预备知识 2.Ocelot基本使用 3.Ocelot功能挖掘 4.Ocelot集成其他组件 5.避坑指 ...
- Tomcat的使用和配置
Tomcat的使用 安装 在tomcat官网找到你需要用的 Tomcat 版本对应的 zip 压缩包,解压到需要安装的目录即可 目录介绍 bin : 专门用来存放Tomcat服务器的可执行文件 con ...
- [bug] Maven [WARNING] 'parent.relativePath' of POM
参考 https://blog.csdn.net/simajinxiu/article/details/86667894
- 【转载】windows linux cent 7 制作U盘 启动盘
1 镜像iso文件存放在linux环境下用dd if=/dev/sdb of=/镜像存放路径/镜像iso文件 bs=1M u盘的盘符是/dev/sdb 2 镜像iso文件存放在windows环境下ul ...