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 ...
随机推荐
- 电脑提示无法装入/加载SolidWorks DLL文件:sldshellutils如何解决
电脑提示无法装入/加载SolidWorks DLL文件:sldshellutils如何解决 参考资料:http://www.xitongcheng.com/jiaocheng/dnrj_article ...
- CCNA 第四章 轻松划分子网
1:划分子网的的好处: (1):减少网络流量 (2):优化网络性能 (3):简化管理 (4):有助于覆盖大型地理区域 2:CIDR和ISP的概念 (1):CIDR:Classless Inter-Do ...
- 关于HTTP的一些概念
各种概念 HTTP HTTP(HyperText Transfer Protocol) -- 超文本传输协议 它可以拆成三个部分:"超文本"."传输".&quo ...
- v-for详解
v-for的引入 当我们需要对一组数据进行渲染时,我们就可以使用v-for来完成 v-for遍历数组 格式:v-for="(item, index) in items".(也许是因 ...
- Redis数据结构—整数集合与压缩列表
目录 Redis数据结构-整数集合与压缩列表 整数集合的实现 整数集合的升级 整数集合不支持降级 压缩列表的构成 压缩列表节点的构成 小结 Redis数据结构-整数集合与压缩列表 大家好,我是白泽.今 ...
- Docker系列——Grafana+Prometheus+Node-exporter服务器监控平台(一)
在最近的博文中,都是介绍监控平台的搭建,其实并不难,主要是需要自己动手操作,实践一番就会了. 有天在想,云上的服务器,是不是也可以搭建一个监控平台,所以就捣鼓了一下,不过遗憾的是,使用阿里云开源的插件 ...
- Linux_rpm包管理
一.rpm包命令规范 1.包的组成 主包:bind-9.7.1-1.el5.i586.rpm 子包:bind-libs-9.7.1-1.el5.i586.rpm bind-utils-9.7.1-1. ...
- linux用户(组)及文件权限说明
写在前面: 在liunx系统中一切都是文件(文件夹和外设是特殊的文件),如果有可能尽量使用文本文件,文本文件是人.机都能够理解的文件.linux的用户(组)和权限完美的解决了linux系统的安全问题. ...
- 路由信息相关 route 网卡
目录 route命令 1.查看路由表 2.管理路由 基本网络配置 添加网卡地址 修改网卡UUID route命令 路由表管理命令,路由表主要构成: Destination: 目标网络ID,表示可以到达 ...
- 西门子S7200/300/400以太网通讯处理器选型分类
北京华科远创科技有限研发的远创智控转以太网模块适用于西门子S7-200/S7-300/S7-400.SMART S7-200.西门子数控840D.840DSL.合信.亿维PLC的PPI/MPI/PRO ...
