题意:
给你n,m。
接下来n个数,给出a[i]。再接下来n个数,给出b[i]。
m次查询,每次查询:
若 id=1 则给出 x y 把a[i]的值修改为y。
若 id=2 则给出 x y 把b[i]的值修改为y。
若 id=3 则给出 k 求,其中
思路:
求S(t)>=k最小的t 显然是二分。
注意S(t)这个式子。我们每次从1~n求的话,复杂度n*m会爆炸。
注意一个诡异的数据范围:1<=ai<=1000
这是解决此题的关键!!!
我们将上式拆开得到t/ai 和 bi/ai 则可以O(n)求出。之后便可以O(1)修改了。
先确定答案为。然后看向下取整的影响。
而影响答案的只是 ,t%ai 和t%bi。不难发现,对于每个i,当t%ai<t%bi答案需要减1。
因此按a[i]建立1000个树状数组,每个树状数组维护分母为a[i],%a[i]>=x的数字个数即可。
同时预处理每个a[i]出现的次数,这样总时间复杂度就降为O(T*m*1000*log(1000)*log(t))。
10s的时限,足以通过此题。如果还不懂,看代码就很容易理解了:
代码:
#include<bits/stdc++.h>
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
const int maxn=1010;
const int maxm=100010;
int n,m;
int lb(int x){return x&(-x);}
struct BIT
{
ll c[maxn];
void init(){memset(c,0,sizeof(c));}
void add(int i,int v)
{
while(i<maxn)
{
c[i]+=v;
i+=lb(i);
}
}
ll query(int i)
{
ll sum=0;
while(i>0)
{
sum+=c[i];
i-=lb(i);
}
return sum;
}
}bit[maxn];
ll sm,b[maxm],a[maxm];
ll cnt[maxn];
bool jud(ll x,ll y)
{
ll sum=0;
for(ll i=1;i<=1000;i++)
{
sum+=(x/i)*cnt[i]-(cnt[i]-bit[i].query(x%i+1));
if(sum>=y) return 1;
}
return 0;
}
int main()
{
int T,cas=1;
scanf("%d",&T);
while(T--)
{
scanf("%d%d",&n,&m);
memset(cnt,0,sizeof(cnt));
for(int i=1;i<=1000;i++)
bit[i].init();
for(int i=1;i<=n;i++)
{
scanf("%lld",&a[i]);
cnt[a[i]]++;
}
sm=0;
for(int i=1;i<=n;i++)
{
scanf("%lld",&b[i]);
sm+=(b[i]/a[i]);
bit[a[i]].add(b[i]%a[i]+1,1);
}
while(m--)
{
int id;
int x,y;
scanf("%d",&id);
if(id==1)
{
scanf("%d%d",&x,&y);
cnt[a[x]]--;
sm=sm-(b[x]/a[x])+(b[x]/(ll)y);
bit[a[x]].add(b[x]%a[x]+1,-1);
cnt[y]++;
bit[y].add(b[x]%y+1,1);
a[x]=y;
}
else if(id==2)
{
scanf("%d%d",&x,&y);
sm=sm-(b[x]/a[x])+((ll)y/a[x]);
bit[a[x]].add(b[x]%a[x]+1,-1);
bit[a[x]].add(y%a[x]+1,1);
b[x]=y;
}
else
{
ll y;
scanf("%lld",&y);
y+=sm;// !!!
ll l=0,r=1e15;
while(l<r)
{
ll mid=(l+r)>>1;
if(jud(mid,y)) r=mid;
else l=mid+1;
}
printf("%lld\n",r);
}
}
}
return 0;
}