装备——题解

题目大意：
给出两个序列，期间会对其中几个数进行修改
“不太多次”求两序列中的子序列（$[L1,R1]$,$[L2,R2]$）两两相乘得到的结果的第K大元素
$1≤n1≤250,1≤n2≤100000,1≤Q≤100000,询问次数≤200，4s$

将所有答案都算出来再求其中的第K大的$O(250*100000*200*log)$显然太大了
实际上这题就是典型的求两序列乘积的第K大值
我们可以像猜数一样直接二分“猜”最后的答案，修改直接暴力$O(1)$修改
然后进行快速check，检验它是否就是小于等于K个数
——那么，排序+二分就OK了：枚举一个序列，在另一个序列中二分“临界点”

int count(int x){
    int S=0;
    for(int i=1;i<=aa[0];i++){
        int L=1,R=bb[0],mid;
        while(L<=R){
            mid=L+R>>1;
            if(aa[i]*bb[mid]<x) L=mid+1;else R=mid-1;
        }
        S+=bb[0]-R;
    }
    return S;
}

再仔细一分析，aa,bb都单调不降，则我们每次二分出来的位置显然也是单调不升的
也就是说，随着aa[i]的增大，“临界点”所在位置只会在上一次基础上适当下降，于是成功又消掉一个$log$！！！
复杂度近似$O($$log_{ans}$$*N*Q)$

#pragma GCC optimize(6)
#include<cstdio>
using namespace std;
const int maxa=255,maxb=(1e5)+5;
int n,m,q,k,a[maxa],b[maxb],aa[maxa],bb[maxb],hsh[maxb*10];
char gt(){
    static char buf[100000],*p1=buf,*p2=buf;
    return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
int read(){
    int ret=0;bool f=0;char ch=gt();
    while(ch<'0'||ch>'9') f|=(ch=='-'),ch=gt();
    while(ch>='0'&&ch<='9') ret=ret*10+ch-'0',ch=gt();
    return f?-ret:ret;
}
//int count(int x){
//  int S=0;
//  for(int i=1;i<=aa[0];i++){
//      int L=1,R=bb[0],mid;
//      while(L<=R){
//          mid=L+R>>1;
//          if(aa[i]*bb[mid]<x) L=mid+1;else R=mid-1;
//      }
//      S+=bb[0]-R;
//  }
//  return S;
//}
int count(int x){
    int S=0;
    for(int i=1,j=bb[0];i<=aa[0];i++){
        while(j&&aa[i]*bb[j]>=x) j--;
        S+=bb[0]-j;
    }
    return S;
}
int work(){
    int L=aa[1]*bb[1],R=aa[aa[0]]*bb[bb[0]],mid,ans;
    while(L<=R){
        mid=(R-L>>1)+L;//小心爆出
        if(count(mid)>=k) ans=mid,L=mid+1;else R=mid-1;
    }
    return ans;
}
int main(){
    n=read(),m=read();
    for(int i=1;i<=n;i++) a[i]=read();
    for(int i=1;i<=m;i++) b[i]=read();
    q=read();
    while(q--){
        bool flg=read();
        if(!flg){
            flg=read();int x=read(),y=read();
            if(!flg) a[x]=y;else b[x]=y;
        }else{
            int L=read(),R=read(),l=read(),r=read();k=read();
            aa[0]=bb[0]=0;
            for(int i=L;i<=R;i++) hsh[a[i]]++;for(int i=1;i<=2000;i++) while(hsh[i]) aa[++aa[0]]=i,hsh[i]--;
            for(int i=l;i<=r;i++) hsh[b[i]]++;for(int i=1;i<=1000000;i++) while(hsh[i]) bb[++bb[0]]=i,hsh[i]--;
            printf("%d\n",work());
        }
    }
    return 0;
}

PS：二分的等号千万注意，绝对整死人！！！