19-10-29-Night-X
布谷。
欢迎大家来不耻下问。
这里是正解不会,暴力写跪,乱搞鸡肋的某虻 。
好想放假啊!!!!
话说猫厂现在产量低迷……
ZJ一下:
T1,三维偏序,只码了$\Theta(N^2)$暴力
T2,暴力愉快。
T3,输〇片分。
简单TJ:
T1
因为这个题的要求比较特殊,求解的答案可以和一个输入的顺序一起维护。
首先通过排序将$A$数组的偏序限制解决
然后用树状数组暴干$B$数组的大小,将位置插入树状数组比较维护。
最后的答案直接查即可。
#include <algorithm>
#include <iostream>
#include <cstring>
#include <climits>
#include <vector>
#include <cstdio>
#define N 555555
#define LL long long
using namespace std;
const int lim=500000;
LL pn;
struct PRE{
LL id,pa,pb;
}ps[N];
vector<LL>vls;
LL pre[N];
LL ans=0;
inline LL fvind(LL a){
return lower_bound(vls.begin(),vls.end(),a)-vls.begin()+1;
}
inline LL lowbit(int x){
return x&(-x);
}
void add(LL pos,LL val){
while(pos<=lim){
pre[pos]=min(pre[pos],val);
pos+=lowbit(pos);
}
}
LL query(LL pos){
LL res=LLONG_MAX;
while(pos){
res=min(res,pre[pos]);
pos-=lowbit(pos);
}
return res;
}
int main(){
LL a;
ios_base::sync_with_stdio(false);
cin>>pn;
for(int i=1;i<=pn;i++)ps[i].id=i;
for(int i=1;i<=pn;i++){
cin>>a;
ps[i].pa=ps[i-1].pa+a;
}
vls.push_back(0);
for(int i=1;i<=pn;i++){
cin>>a;
ps[i].pb=ps[i-1].pb+a;
vls.push_back(ps[i].pb);
}
// for(auto i:vls){\
cout<<i<<" ";\
}cout<<endl;
sort(vls.begin(),vls.end());
vls.erase(unique(vls.begin(),vls.end()),vls.end());
// for(auto i:vls){\
cout<<i<<" ";\
}cout<<endl;
for(int i=0;i<=pn;i++){
ps[i].pb=fvind(ps[i].pb);
}
sort(ps,ps+pn+1,[](const PRE &a,const PRE &b)
{return a.pa==b.pa?a.id<b.id:a.pa<b.pa;});
// for(int i=1;i<=pn;i++){\
cout<<ps[i].pa<<" "<<ps[i].pb<<" "<<ps[i].id<<endl;\
}
memset(pre,0x7f,sizeof pre);
// add(ps[0].pb,0);
for(int i=0;i<=pn;i++){
LL vl=query(ps[i].pb);//cout<<ps[i].id<<"val:"<<vl<<endl;
add(ps[i].pb,ps[i].id);
ans=max(ans,ps[i].id-vl);
}
cout<<ans<<endl;
}
T2
证明先咕下。
利用决策单调性,将$\Theta(N^3)$优化成$\Theta(N^2)$
#include <iostream>
#include <climits>
#include <cstring>
#include <cstdio>
#define N 5555
#define LL long long
using namespace std;
int pn;
LL arr[N],pre[N],
dp[N][N],
ap[N][N];
int main(){
cin.sync_with_stdio(false);
cin>>pn;
for(int i=1;i<=pn;i++){
cin>>arr[i];
pre[i]=pre[i-1]+arr[i];
}
for(int i=1;i<=pn;i++){
for(int j=1;j<=pn;j++){
if(i<=j)dp[i][j]=LLONG_MAX/3;
else dp[i][j]=0;
}
}
for(int i=1;i<=pn;i++){
dp[i][i]=arr[i];
ap[i][i]=i;
}
for(int i=1;i<pn;i++){
for(int j=1;j+i<=pn;j++){
int l=j,r=j+i;
// cout<<l<<" "<<r<<endl;
// cout<<ap[l][r-1]<<"~"<<ap[l+1][r]<<endl;
for(int k=ap[l][r-1];k<=ap[l+1][r];k++){
// cout<<dp[l][r]<<" <<<< "<<dp[l][k-1]<<"+"<<dp[k+1][r]<<endl;
if(dp[l][r]>dp[l][k-1]+dp[k+1][r]){
dp[l][r]=dp[l][k-1]+dp[k+1][r];
ap[l][r]=k;
}
}
dp[l][r]+=pre[r]-pre[l-1];
}
}
cout<<dp[1][pn]<<endl;
}
T3
大神题,好像和一些神奇的高斯消元有关。
Miemeng真的蒻

浙公网安备 33010602011771号