HDU 6274 二分+预处理（CCPC K题

 

#include"bits/stdc++.h"
#define db double
#define ll long long
#define vec vector<ll>
#define Mt  vector<vec>
#define ci(x) scanf("%d",&x)
#define cd(x) scanf("%lf",&x)
#define cl(x) scanf("%lld",&x)
#define pi(x) printf("%d\n",x)
#define pd(x) printf("%f\n",x)
#define pl(x) printf("%lld\n",x)
using namespace std;
const int N   = 1e6 + 5;
int a[N],b[N],s[1010][1010],n; //s[x][y]表示对于所有a[i]==x时的b[i]%a[i]余数大于等于y的数目
ll ret;
bool check(ll x,ll k)
{
    ll ans=ret*-1;
    for(int i=1; i<=1000; i++)
    {
        ans+=x/i*s[i][0]; //s[i][0]为a数组中等于i的数的个数
        ans-=s[i][x%i+1];  //减掉所有 c1<c2 的情况
    }
    return ans>=k;
}
ll cal(ll x)
{
    ll l=1,r=1e13,mid;
    ll res=l;
    while(l<=r)
    {
        mid=(l+r)>>1;
        if(check(mid,x)) res=mid,r=mid-1;
        else l=mid+1;
    }
    return res;
}
int main()
{
    int TA,T,x,y,z;
    scanf("%d",&TA);
    while(TA--)
    {
        scanf("%d%d",&n,&T);
        ret=0;
        memset(s,0, sizeof(s));
        for(int i=1; i<=n; i++) scanf("%d",&a[i]);
        for(int i=1; i<=n; i++) scanf("%d",&b[i]),ret+=(b[i]/a[i]),s[a[i]][b[i]%a[i]]++;
        for(int i=1; i<=1000; i++)
            for(int j=i-1; j>=0; j--)
                s[i][j]+=s[i][j+1];
        while(T--)
        {
            scanf("%d%d",&z,&x);
            if(z==1)
            {
                scanf("%d",&y);
                for(int i=b[x]%a[x]; i>=0; i--) s[a[x]][i]--;
                for(int i=b[x]%y; i>=0; i--)    s[y][i]++;
                ret-=(b[x]/a[x]);
                ret+=(b[x]/y);
                a[x]=y;
            }
            else if(z==2)
            {
                scanf("%d",&y);
                for(int i=b[x]%a[x]; i>=0; i--) s[a[x]][i]--;
                for(int i=y%a[x]; i>=0; i--)    s[a[x]][i]++;
                ret-=(b[x]/a[x]);
                ret+=(y/a[x]);
                b[x]=y;
            }
            else
            {
                printf("%lld\n",cal(1ll*x));
            }
        }
    }
}