T1 Adjoin 【問題描述】 定義一種合法的$0 1$串:串中任何一個數字都與$1$相鄰。例如長度為$ 3 的 0 1 $串中,$101$是非法的,因為兩邊的$1$沒有相鄰的$1,011$是合法的,因為三個數都有$1$相鄰。現在問,長度為$N$的$0 1$中有多少是合法的。 【輸入格式】 一行, ...
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}
\]
其中我們將flag的初值賦為1,在每碰到一個奇數時為其取反即可。
#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\)有點聯繫。然而矩陣優化要求使用通項公式,而我們當前只有一個遞推式。那麼我們現在考慮將方程式反向優化,從一維方程變為三維方程,使整個式子具有通項公式,再使用矩陣優化來降低整體的複雜度。想到了這一點後,實現起來應該並不難。
#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分的暴力。
附六十分的暴力寫法
#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。
#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。
#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 $
【題解】
模擬題,註意一下優化就好。本蒟蒻的代碼風格太醜因此在此不予貼出。