jzoj5439【NOIP2017提高A组集训10.31】Calculate

题目

me Limits: 4000 ms Memory Limits: 131072 KB

Description

Input

Output

Sample Input

见下载

Sample Output

见下载

Data Constraint

题解

观察式子，感觉那个下取整很难搞啊，那么转化一下形式
原来的式子也可以表示为t/a[i]-b[i]/a[i]-[b[i]%a[i]>t%a[i]]（这条式子里的“/”都表示下取整）
证明的话我们可以划一条数轴，发现两个模完之后的值就相当于a[i],t相对于相同值顶点的距离，如果他们的距离小于a[i]那么就要多减去一个1
n有一点大，不好搞，但是观察到a十分的小，所以我们可以按照a的值分组
前面的两项直接维护就好了，第三项的话我们可以设sum[i][j]表示a[x]的值为i，且b[x]%a[x]>=j的个数，然后我们就可以通过O(1000)的时间算出一个S(T)
前面两种操作的话我们暴力的维护就可以了

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<cmath>
#define fo(i,a,b) for(i=a;i<=b;i++)
#define fo1(i,b,a) for(i=b;i>=a;i--)
#define ll long long

using namespace std;

const int maxn=1e5+5;
const int cs=1e9;
int p[1005][1005],sum[1005][1005];
ll sc[1005],s[1005],a[maxn],b[maxn];
ll i,j,k,l,r,m,n,t,o,x,y,z,an,mid;
char ch;

int read(){
    int x=0;
    ch=getchar();
    while (ch<'0' || ch>'9') ch=getchar();
    while (ch>='0' && ch<='9'){
        x=x*10+ch-48;
        ch=getchar();
    }
    return x;
}
void remake1(int x){
    int i;
    fo1(i,x-1,0) sum[x][i]=sum[x][i+1]+p[x][i];
}
void ge1(){
    sc[x]++;
    sc[a[i]]--;
    p[a[i]][b[i]%a[i]]--;
    p[x][y%x]++;
    s[a[i]]=s[a[i]]-b[i]/a[i];
    s[x]=s[x]+y/x;
    remake1(a[i]);
    remake1(x);
}
int main(){
    freopen("calculate.in","r",stdin);
    freopen("calculate.out","w",stdout);
    t=read();
    fo(o,1,t){
        memset(p,0,sizeof(p));
        memset(s,0,sizeof(s));
        memset(sc,0,sizeof(sc));
        n=read(); m=read();
        fo(i,1,n) a[i]=read();
        fo(i,1,n) b[i]=read();
        fo(i,1,n){
            s[a[i]]=s[a[i]]+b[i]/a[i];
            p[a[i]][b[i]%a[i]]++;
            sc[a[i]]++;
        }
        fo(i,1,1000)
            fo1(j,i-1,0) sum[i][j]=sum[i][j+1]+p[i][j];
        fo(j,1,m){
            x=read();
            if (x==1){
                x=read(); z=read(); i=x; y=b[x]; x=z; 
                ge1();
                a[i]=x;
            } else if (x==2){
                x=read(); y=read(); i=x; x=a[x];
                ge1();
                b[i]=y;
            } else{
                x=read();
                l=1;
                r=cs;
                while (l<r){
                    mid=(l+r)/2;
                    an=0;
                    fo(i,1,500)
                        an=an+sc[i]*(mid/i)-s[i]-sum[i][mid%i+1];
                    if (an>=x) r=mid; else l=mid+1;
                }
                printf("%d\n",l);
            }
        }
    }
    return 0;
}