牛客網noip賽前集訓 提高組第4場 T2 區間 單調求解法,卡常
阿新 • • 發佈：2018-11-10
文章目錄
非常簡單的一場了,然而我T1巧妙地寫萎把 給爆了,掛成了 分.
不過T2非常好玩,我介紹介紹.
題意
分做法.
一開始寫了個
分.可以發現的性質是,一個序列
的
.
所以一個序列的
如果要出現在這個序列中,必然是最小的數值.
我們
預處理每個區間的
和最小值,再
求出答案.
#include<bits/stdc++.h> //Ithea Myse Valgulious
namespace chtholly{
typedef long long ll;
#define re0 register int
#define rel register ll
#define rec register char
#define gc getchar
#define pc putchar
#define p32 pc(' ')
#define pl puts("")
/*By Citrus*/
inline int read(){
int x=0,f=1;char c=gc();
for (;!isdigit(c);c=gc()) f^=c=='-';
for (;isdigit(c);c=gc()) x=(x<<3)+(x<<1)+(c^'0');
return f?x:-x;
}
template <typename mitsuha>
inline bool read(mitsuha &x){
x=0;int f=1;char c=gc();
for (;!isdigit(c)&&~c;c=gc()) f^=c=='-';
if (!~c) return 0;
for (;isdigit(c);c=gc()) x=(x<<3)+(x<<1)+(c^'0');
return x=f?x:-x,1;
}
template <typename mitsuha>
inline int write(mitsuha x){
if (!x) return 0&pc(48);
if (x<0) pc('-'),x=-x;
int bit[20],i,p=0;
for (;x;x/=10) bit[++p]=x%10;
for (i=p;i;--i) pc(bit[i]+48);
return 0;
}
inline char fuhao(){
char c=gc();
for (;isspace(c);c=gc());
return c;
}
}using namespace chtholly;
using namespace std;
const int yuzu=4e6,aoi=2038;
typedef ll fuko[yuzu|10];
fuko a,tmp;
ll gcd[aoi][aoi],xiao[aoi][aoi];
int main(){
int i,j,k,n=read(),llx=0;
memset(xiao,0x3f,sizeof xiao);
for (i=1;i<=n;++i) read(a[i]);
for (i=1;i<=n;++i){
for (j=i;j<=n;++j){
gcd[i][j]=__gcd(gcd[i][j-1],a[j]);
xiao[i][j]=min(xiao[i][j-1],a[j]);
}
}
for (i=1;i<=n;++i){
for (j=i;j<=n;++j){
if (gcd[i][j]==xiao[i][j]) llx=max(llx,j-i+1);
}
} write(llx);
}
分做法
我們定義
為從
向左能到達的離它最近的不能被它整除的數的位置
.
同理定義
為從
向右能到達的離它最近的不能被它整除的數的位置
.
這樣答案就是
.
如何
求出
這兩個陣列呢?
可以注意到,如果
,
可以直接跳到
.
這樣我們可以像
一樣利用已經求出來的
求出接下來的
.
均攤複雜度是
.
for (l[i]=i;l[i]>1&&a[l[i]-1]%a[i]==0;l[i]=l[l[i]-1]);
這樣便可以獲得
分.
#include<bits/stdc++.h> //Ithea Myse Valgulious
namespace chtholly{
typedef long long ll;
#define re0 register int
#define rel register ll
#define rec register char
#define gc getchar
#define pc putchar
#define p32 pc(' ')
#define pl puts("")
/*By Citrus*/
inline int read(){
int x=0,f=1;char c=gc();
for (;!isdigit(c);c=gc()) f^=c=='-';
for (;isdigit(c);c=gc()) x=(x<<3)+(x<<1)+(c^'0');
return f?x:-x;
}
template <typename mitsuha>
inline bool read(mitsuha &x){
x=0;int f=1;char c=gc();
for (;!isdigit(c)&&~c;c=gc()) f^=c=='-';
if (!~c) return 0;
for (;isdigit(c);c=gc()) x=(x<<3)+(x<<1)+(c^'0');
return x=f?x:-x,1;
}
template <typename mitsuha>
inline int write(mitsuha x){
if (!x) return 0&pc(48);
if (x<0) pc('-'),x=-x;
int bit[20],i,p=0;
for (;x;x/=10) bit[++p]=x%