题目链接:D. The Child and Sequence
题解
本题的难点是取模。
例如:
x
%
i
{x\%i}
x%i
- 如果 i ≤ ⌈ x 2 ⌉ {i≤\lceil{\frac{x}{2}} \rceil} i≤⌈2x⌉,那么余数必定小于除数 ⌈ x 2 ⌉ {\lceil{\frac{x}{2}} \rceil} ⌈2x⌉
- 如果 x ≥ i > ⌈ x 2 ⌉ {x≥i> \lceil{\frac{x}{2}} \rceil} x≥i>⌈2x⌉,那么余数为x-i仍然小于 ⌈ x 2 ⌉ {\lceil{\frac{x}{2}} \rceil} ⌈2x⌉
综上所述,只要除数≤x,余数必定减半。
所以就暴力更新即可,每个数最多被更新
l
o
g
(
a
i
)
{log(a_i)}
log(ai)次。
维护一个区间最大值,如果当前区间最大值小于取模的数就直接返回,因为取模比自身大的数还是自己本身。
时间复杂度 O ( ( n + m ) l o g a i l o g n ) {O((n+m)loga_i logn)} O((n+m)logailogn)
代码
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<bitset>
#include<cassert>
#include<cctype>
#include<cmath>
#include<cstdlib>
#include<ctime>
#include<deque>
#include<iomanip>
#include<list>
#include<map>
#include<queue>
#include<set>
#include<stack>
#include<vector>
#include<unordered_set>
#include<unordered_map>
using namespace std;
//extern "C"{void *__dso_handle=0;}
typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
#define fi first
#define se second
#define pb push_back
#define mp make_pair
#define pii pair<int,int>
#define lowbit(x) x&-x
const double PI=acos(-1.0);
const double eps=1e-6;
const ll mod=1e9+7;
const int inf=0x3f3f3f3f;
const int maxn=1e5+10;
const int maxm=1e6+10;
#define ios ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
ll a[maxn];
struct Tree
{
ll maxx,summ;
}tree[maxn<<2];
void pushup(int p,int l,int r)
{
tree[p].summ=tree[p*2].summ+tree[p*2+1].summ;
tree[p].maxx=max(tree[p*2].maxx,tree[p*2+1].maxx);
}
void build(int p,int l,int r)
{
if(l==r) { tree[p].summ=tree[p].maxx=a[l]; return ; }
int mid=(l+r)>>1;
build(p*2, l, mid);
build(p*2+1, mid+1, r);
pushup(p,l,r);
}
void update_m(int p,int l,int r,int ql,int qr,int x)
{
if(x>tree[p].maxx) return ;
if(l==r) { tree[p].summ%=x; tree[p].maxx%=x; return ;}
int mid=(l+r)>>1;
if(ql<=mid) update_m(p*2, l, mid, ql, qr, x);
if(qr>mid) update_m(p*2+1, mid+1, r, ql, qr, x);
pushup(p, l, r);
}
void update_c(int p,int l,int r,int k,int x)
{
if(l==r)
{
tree[p].summ=tree[p].maxx=(ll)x;
return ;
}
int mid=(l+r)>>1;
if(k<=mid) update_c(p*2, l, mid, k, x);
if(k>mid) update_c(p*2+1, mid+1, r, k, x);
pushup(p, l, r);
}
ll query(int p,int l,int r,int ql,int qr)
{
if(ql<=l && qr>=r) return tree[p].summ;
int mid=(l+r)>>1;
ll ans=0;
if(ql<=mid) ans+=query(p*2, l, mid, ql, qr);
if(qr>mid) ans+=query(p*2+1, mid+1, r, ql, qr);
pushup(p, l, r);
return ans;
}
int main()
{
int n,m;
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++) scanf("%lld",&a[i]);
build(1, 1, n);
while (m--) {
int op; scanf("%d",&op);
if(op==1)
{
int l,r; scanf("%d%d",&l,&r);
printf("%lld\n",query(1, 1, n, l, r));
}
else if(op==2)
{
int l,r,x; scanf("%d%d%d",&l,&r,&x);
update_m(1, 1, n, l, r, x);
}
else {
int k,x; scanf("%d%d",&k,&x);
update_c(1, 1, n, k, x);
}
}
}