CCF-NOIP-2018 提高组(复赛) 模拟试题(七)
T1 Adjoin
【问题描述】
定义一种合法的\(0-1\)串:串中任何一个数字都与\(1\)相邻。例如长度为$ 3 的 0-1 $串中,\(101\)是非法的,因为两边的\(1\)没有相邻的\(1,011\)是合法的,因为三个数都有\(1\)相邻。现在问,长度为\(N\)的\(0-1\)中有多少是合法的。
【输入格式】
一行,一个整数\(N\)
【输出格式】
一行,合法\(0-1\)串的个数,答案对\(1000000007\)取模。
【样例1】
样例输入 |
---|
3 |
样例输出 |
---|
3 |
样例说明
\(110、111、011\)是所有长度为\(3\)的合法\(0-1\)串
数据规模与约定
所有测试点的数据规模与约定如下:
\(30\%\)输入数据,保证\(n<=50。\)
\(60\%\)输入数据,保证$ n<=5000$
\(80\%\)输入数据,保证$ n<=1000000$
\(100\%\)输入数据,保证$ 1<=n<=10^{18}$
题解
一道很经典的需要反向优化的题目。我们首先考虑暴搜得到较小范围内每一个\(n\)所对应的答案,如下所示
i | f[i] |
---|---|
1 | 0 |
2 | 1 |
3 | 3 |
4 | 4 |
5 | 5 |
6 | 9 |
7 | 16 |
8 | 25 |
9 | 29 |
10 | 54 |
然而直接观察数据似乎没有什么明显的规律。于是我们选择将奇偶数分开判断。经过一段时间的观察,似乎所有的数满足这样一个规律:$$f[i]=\begin{cases}
0 & i=1\
1 & i=2\
f[i-1]+f[i-2] & i%2=0\
f[i-1]+f[i-2]+(flag=1)?-2:2 & i%2=1\
\end{cases}
```cpp
#include<bits/stdc++.h>
#define maxn 1000005
const long long mod = 1000000007;
using namespace std;
long long n;
long long dp[maxn];
bool flag=1;
int main(){
//freopen("adjoin.in","r",stdin);
//freopen("adjoin.out","w",stdout);
ios::sync_with_stdio(false);
cin.tie(0);
cin>>n;
dp[1]=0;
dp[2]=1;
for(register long long i=3;i<=n;i++){
if(i%2==0)dp[i]=dp[i-1]+dp[i-2];
else{
if(flag){
dp[i]=dp[i-1]+dp[i-2]+2;
flag=0;
}
else{
dp[i]=dp[i-1]+dp[i-2]-2;
flag=1;
}
}
//cout<<dp[i]<<endl;
dp[i]%=mod;
}
//for(register long long i=3;i<=n;i++)cout<<dp[i]<<endl;
cout<<dp[n]%mod<<endl;
return 0;
}
```
这时我们就拥有了85分的总分。最大数据范围为$n\le 10^{18}$,早已超出了$O(n)$的复杂度能到达的极限。此时,我们思考所有和动态规划有关的优化。经过一番思索后,我们会发现只有矩阵优化稍微与目前的$dp$有点联系。然而矩阵优化要求使用通项公式,而我们当前只有一个递推式。那么我们现在考虑将方程式**反向优化,从一维方程变为三维方程,使整个式子具有通项公式,再使用矩阵优化来降低整体的复杂度**。想到了这一点后,实现起来应该并不难。
```cpp
#include <bits/stdc++.h>
#define RG register
using namespace std;
typedef long long ll;
const int maxn = 1000010;
const int mod = 1e9+7;
const int P = 1e9+7;
ll n;
struct Matrix {
ll val[4][4];
} A, I, ans;
Matrix operator*(const Matrix &A,const Matrix &B){
Matrix C;
for(int i = 0;i < 4;i++)
for(int j = 0;j < 4;j++){
C.val[i][j] = 0;
for(int k = 0;k < 4;k++)
C.val[i][j] = (1ll * A.val[i][k] * B.val[k][j] + C.val[i][j]) % P;
}
return C;
}
Matrix fpm(Matrix x, long long y){
Matrix ret = I;
while(y){
if(y & 1) ret = ret * x;
x = x*x;
y >>= 1;
}
return ret;
}
int main(){
freopen("adjoin.in","r",stdin);
freopen("adjoin.out","w",stdout);
scanf("%lld",&n);
if(n == 1){
puts("0");
return 0;
}
for(int i = 0;i < 4;i++){
for(int j = 0;j < 4;j++){
I.val[i][j] = (i == j ? 1 : 0);
A.val[i][j] = 0;
}
}
ans.val[0][0] = 0;
ans.val[0][1] = 1;
ans.val[0][2] = 0;
ans.val[0][3] = 1;
A.val[2][0] = A.val[0][1] = A.val[2][1] = A.val[3][2] = A.val[1][3] = A.val[3][3] = 1;
A = fpm(A,n-2);
ans = ans * A;
ll ansss = (ans.val[0][2] + ans.val[0][3]) % mod;
printf("%lld",ansss);
return 0;
}
```
### T2 Sorting
#### 【问题描述】
小F不喜欢递减,他会想尽一切办法将看到的一切东西排序!
现在小F得到了一个数列,他当然要将这个数列排序了,但他太累了,以至于最多只能交换其中两个元素,如果这样不能使得这个数列不递减,他就要先睡觉了。你能告诉他是否可行吗?
#### 【输入格式】
第一行:一个整数N表示小F的数列中数的个数。
第二行,N个正整数,描述小F的数列。
#### 【输出格式】
一行,YES或NO,表示通过一次“最多交换其中两个元素(可以不交换)”的操作,是否可使得小F的数列不递减。
#### 【样例1】
|样例输入
|-
|3
|1 3 1
|样例输出
|-
|YES
#### 数据规模与约定
$30\%的数据,N ≤ 10^2 。$
$60\%的数据,N ≤ 10^4 。$
$100\%的数据,N ≤ 10^5 ,所有数为正整数且在longint范围内。$
#### 题解
是的,本蒟蒻又一次和AC插肩而过。拿到这道题时,大多数人都会联想到逆序对和最长不降子序列问题,然而几组充分设置的样例就可以卡掉这两种思路,使得其得分甚至不如60分的暴力。
附六十分的暴力写法
```cpp
#include<bits/stdc++.h>
using namespace std;
int main(){
puts("YES");
return 0;
}
```
对于**LIS**来说,改变一个数有时可以导致多对大小关系的改变;对于逆序对来说,逆序对个数不一定可以决定最终大小关系。对两种思路最好的反驳样例如下:
|5 3 2
|-
在排除了这样的思路后,蒟蒻的我开始思考自己的做法。我们直接从头往后寻找第一对不满足条件的组合$a_i,a_{i-1}$。此时我们取出$a_i$,从头往后将其与第一个大于他的值交换。此时我们再重新在原串中查找是否存在不合法情况,若存在则输出NO,否则输出YES。
```cpp
#include<bits/stdc++.h>
#define maxn 100005
using namespace std;
inline char get(){
static char buf[30],*p1=buf,*p2=buf;
return p1==p2 && (p2=(p1=buf)+fread(buf,1,30,stdin),p1==p2)?EOF:*p1++;
}
inline long long read(){
register char c=get();register long long f=1,_=0;
while(c>'9' || c<'0')f=(c=='-')?-1:1,c=get();
while(c<='9' && c>='0')_=(_<<3)+(_<<1)+(c^48),c=get();
return _*f;
}
long long n;
long long a[maxn];
long long ans=0;
bool cd=1;
int main(){
//freopen("sorting.in","r",stdin);
//freopen("sorting.out","w",stdout);
n=read();
for(register long long i=1;i<=n;i++)a[i]=read();
for(register long long i=2;i<=n;i++){
if(a[i]<a[i-1]){
for(register long long j=1;j<=n;j++){
if(a[j]>a[i]){
swap(a[i],a[j]);
goto next;
}
}
}
}
next:;
for(register long long i=2;i<=n;i++){
if(a[i]<a[i-1]){
puts("NO");
return 0;
}
}
puts("YES");
return 0;
}
```
为什么我们选择了这样一个奇怪的算法呢?事实上,在比赛中选择算法的第一原则是能否证明其错误性。在这道题中,蒟蒻无法证明该算法是错误的,于是就这么得到了85分的安慰分。
其实题目已经提示得很明显了。Sorting,就是在暗示我们进行一次排序操作。我们只需要比较排序前后的两个两个序列,若其中同一位置不一样的元素的个数在两个以内(一次交换最多导致4对大小关系发生改变),则输出YES,否则就记其为非法情况,输出NO。
```cpp
#include <bits/stdc++.h>
using namespace std;
int b[2111111],a[2111111];
int n;
int main(){
freopen("sorting.in","r",stdin);
freopen("sorting.out","w",stdout);
cin>>n;
for(int i=1;i<=n;i++)cin>>a[i],b[i]=a[i];
sort(b+1,b+1+n);
int len=0;
for(int i=1;i<=n;i++){
if(a[i]!=b[i])len++;
if(len>2){
cout<<"NO"<<endl;
return 0;
}
}
cout<<"YES"<<endl;
return 0;
}
```
### T3 Editor
#### 【问题描述】
小$F$有一个梦想:为数列写一个最强大的编辑器!
一开始,数列为空,光标在开头位置,小F的编辑器要对这个数列作如下五种操作:
**I x**:在光标的后面插入一个数字x,并将光标移到这个新加入的$x$后。
**D**:删除光标前的最后一个数字(保证存在),光标位置不变。
**L**:光标左移一位,如果已经在开头则不做任何事。
**R**:光标右移一位,如果已经在结尾则不做任何事。
**Q k**:编辑器需要给出$A_1,A_2 ··· A_k$的最大前缀和(前缀长度不能为0),保证$1 ≤ k ≤ N$,其中$N$为当
前光标前的数字个数。
#### 【输入格式】
第一行,一个整数$Q$,表示操作的总次数。
后$Q$行,每行是上列五种操作中的一种。
#### 【输出格式】
对每个$Q$操作,输出一行一个整数,表示答案。
#### 【样例输入1】
|
|-
|8
|I 2
|I -1
|I 1
|Q 3
|L
|D
|R
|Q 2
#### 【样例输出1】
|样例输出
|:-
|2
|3
#### 【数据规模与约定】
$ 1\le q \le 1000000,-1000 \le m \le 1000 $
#### 【题解】
模拟题,注意一下优化就好。本蒟蒻的代码风格太丑因此在此不予贴出。\]
CCF-NOIP-2018 提高组(复赛) 模拟试题(七)的更多相关文章
- CCF-NOIP-2018 提高组(复赛) 模拟试题(九)(2018 CSYZ长沙一中)
T1 Circle [问题描述] 小 w 的男朋友送给小 w 一个 n 个点 m 条边的图,并且刁难小 w 要她找出点数最少的正环. 小 w 不会做,于是向你求助. [输入格式] 第一行两个整数\(n ...
- 冲刺NOIP2015提高组复赛模拟试题(五)2.道路修建
2.道路修建 描述 Description liouzhou_101最悲痛的回忆就是NOI2011的道路修建,当时开了系统堆栈,结果无限RE… 出于某种报复心理,就把那题神奇了一下: 在 Z星球上有N ...
- CCF-NOIP-2018 提高组(复赛) 模拟试题(四)
T1 贪吃蛇 [问题描述] 贪吃蛇是一个好玩的游戏.在本题中,你需要对这个游戏进行模拟. 这个游戏在一个 \(n\) 行 \(m\) 列的二维棋盘上进行. 我们用 \((x, y)\) 来表示第 \( ...
- CCF-NOIP-2018 提高组(复赛) 模拟试题(一)
T1 帽子戏法 问题描述 小 Y 有一个\(n*n*n\)的"帽子立方体" ,即一个\(n\)层的立方体,每层的帽子都 可以排成\(n*n\)的矩阵. "帽子立方体&qu ...
- 破译情报-NOIP2016提高组复赛模拟试题
[题目描述] 最近国安人员截获了一份 RB 国的秘密情报, 全文都是经过加密的,每个单 词都很长.破译人员想到先把单词化简一下,方法是把每个单词尽量取短些的前 缀,但所取的前缀不能是其他单词的前缀. ...
- 冲刺NOIP2015提高组复赛模拟试题(五) 3.破坏基地
3.破坏基地 描述 Description 在Z国和W国之间一直战火不断. 好不容易,W国的间谍把完整的Z国的军事基地的地图到手了. 于是W国决定再次出击,一举击破Z国的防线. W国认真研究了Z国的地 ...
- 冲刺NOIP2015提高组复赛模拟试题(五)1.数学作业
1. 数学作业 [问题描述] 路人丙的数学老师非常乏力,他喜欢出一些非常乏力的数学题来为难乏力的学生们.这次数学老师布置了一堆的数学题作为作业,而且这些数学题有个共同的特点是都求C(N,M)中不同质因 ...
- CCF-NOIP-2018 提高组(复赛) 模拟试题(五)
T1 相遇 [问题描述] 在一场奇怪的梦里,小 Y 来到了一个神奇的国度.这个国度可以用一根数轴表示,小 Y 在 N 处,而小 Y 想吃的美食在 K 处.小 Y 有两种方式移动, 一种叫做步行, 一种 ...
- CCF-NOIP-2018 提高组(复赛) 模拟试题(三)
T1 取球游戏 问题描述 现有\(N\)个小球,依次编号为\(1\)到\(N\),这些小球除了编号以外没有任何区别.从这\(N\)个小球中取出\(M\)个,请问有多少种取球方案使得在取出的\(M\)个 ...
随机推荐
- 【luogu P1186 玛丽卡】 题解
题目链接:https://www.luogu.org/problemnew/show/P1186 邻接表开大开小真的奇妙,毒瘤玩意,再您妈的见. #include <queue> #inc ...
- lucene&solr学习——创建和查询索引(理论)
1.Lucene基础 (1) 简介 Lucene是apache下的一个开放源代码的全文检索引擎工具包.提供完整的查询引擎和索引引擎:部分文本分析引擎. Lucene的目的是为软件开发人员提供一个简单易 ...
- How to Create a Basic Plugin
Sometimes you want to make a piece of functionality available throughout your code. For example, per ...
- data-ng-hide指令
<!DOCTYPE html><html><head><meta http-equiv="Content-Type" content=&q ...
- [oracle]分区表学习
(一)什么是分区 所谓分区,就是将一张巨型表或巨型索引分成若干个独立的组成部分进行存储和管理,每一个相对小的,可独立管理的部分,称为分区. (二)分区的优势 提高数据可管理性.对表进行分区,数据的加载 ...
- iOS之查看代码运行的时间
有时候我们想要准确的知道某段代码.某个循环执行的时间,然后分析效率等问题,这个时候就需要执行时间是多少.正好看到网上已经有人做了这个工作,我就直接摘下来了.正好也用了宏的方式计算时间,我们只要在需要计 ...
- zepto 基础知识(2)
20.append append(content) 类型:self 在每个匹配的元素末尾插入内容(内部插入).内容可以为html 字符串.dom节点,或者节点组成的数组. $('ul').append ...
- awk分隔符
最近需要检测日志,shell中用到了awk,因为分割条件不止一个,并且包括了中括号.在此记录一下关于多分隔符并且包含中括号的情况 awk -F'[=,]|[][]+' '{print $6}'
- python 用selenuim判断页面是否全部加载完成,并且加上最大时长,超过时长报错
STR_READY_STATE = '' time_start = time.time() while STR_READY_STATE != 'complete': time.sleep(0.001) ...
- MHA实现mysql高可用复制集群
MHA简述 MHA(Master High Availability)目前在MySQL高可用方面是一个相对成熟的解决方案,是一套优秀的作为MySQL高可用性环境下故障切换和主从提升的高可用软件.在My ...