// begin おためしbeats pub trait Beats { type Value: Clone; fn merge(l: &Self::Value, r: &Self::Value) -> Self::Value; fn propagate(p: &mut Self::Value, c: &mut [Self::Value]); fn e() -> Self::Value; } pub trait Effect { type Value; type Query; fn break_cond(x: &Self::Value, q: &Self::Query) -> bool; fn tag_cond(x: &Self::Value, q: &Self::Query) -> bool; fn update(x: &mut Self::Value, q: &Self::Query); } struct SegmentTreeBeats { size: usize, val: Box<[R::Value]>, dfs: Vec<(usize, usize, usize)>, } impl SegmentTreeBeats { fn new(ini: &[R::Value]) -> Self { let size = ini.len().next_power_of_two(); let val = vec![R::e(); 2 * size]; let mut val = val.into_boxed_slice(); val[size..(size + ini.len())].clone_from_slice(ini); for i in (1..size).rev() { val[i] = R::merge(&val[2 * i], &val[2 * i + 1]); } SegmentTreeBeats { size, val, dfs: vec![], } } fn update>(&mut self, x: usize, y: usize, f: E::Query) { assert!(x <= y && y <= self.size); if x >= y { return; } let val = &mut self.val; self.dfs.push((1, 0, self.size)); while let Some((k, l, r)) = self.dfs.pop() { if k > 2 * self.size { let k = !k; val[k] = R::merge(&val[2 * k], &val[2 * k + 1]); continue; } if E::break_cond(&val[k], &f) { continue; } if x <= l && r <= y && E::tag_cond(&val[k], &f) { E::update(&mut val[k], &f); continue; } let (a, b) = val.split_at_mut(2 * k); R::propagate(&mut a[k], &mut b[..2]); self.dfs.push((!k, l, r)); let m = (l + r) / 2; if m < y { self.dfs.push((2 * k + 1, m, r)); } if x < m { self.dfs.push((2 * k, l, m)); } } } fn find(&mut self, x: usize, y: usize) -> R::Value { assert!(x <= y && y <= self.size); let mut res = R::e(); if x >= y { return res; } let val = &mut self.val; self.dfs.push((1, 0, self.size)); while let Some((k, l, r)) = self.dfs.pop() { if x <= l && r <= y { res = R::merge(&res, &val[k]); continue; } let (a, b) = val.split_at_mut(2 * k); R::propagate(&mut a[k], &mut b[..2]); let m = (l + r) / 2; if m < y { self.dfs.push((2 * k + 1, m, r)); } if x < m { self.dfs.push((2 * k, l, m)); } } res } } // end おためしbeats // ---------- begin binary_gcd ---------- pub fn gcd(a: u64, b: u64) -> u64 { if a == 0 || b == 0 { return a + b; } let x = a.trailing_zeros(); let y = b.trailing_zeros(); let mut a = a >> x; let mut b = b >> y; while a != b { let x = (a ^ b).trailing_zeros(); if a < b { std::mem::swap(&mut a, &mut b); } a = (a - b) >> x; } a << x.min(y) } // ---------- end binary_gcd ---------- // ---------- begin Scanner(require delimiter) ---------- mod scanner { pub struct Scanner { reader: R, buf: Vec, } #[allow(dead_code)] impl Scanner { pub fn new(reader: R) -> Self { Scanner { reader: reader, buf: Vec::with_capacity(1024), } } fn read(&mut self, del: u8) { self.buf.clear(); self.reader.read_until(del, &mut self.buf).ok(); assert!(self.buf.pop().unwrap() == del); } pub fn next(&mut self, del: u8) -> T { self.read(del); std::str::from_utf8(&self.buf) .unwrap() .trim() .parse::() .ok() .unwrap() } pub fn next_bytes(&mut self, del: u8) -> Vec { self.read(del); std::str::from_utf8(&self.buf) .unwrap() .trim() .bytes() .collect() } } } // ---------- end scanner(require delimiter) ---------- use std::io::Write; fn main() { let stdin = std::io::stdin(); let mut sc = scanner::Scanner::new(stdin.lock()); run(&mut sc); } const INF: u64 = 1_000_000_000 + 1; fn lcm(a: u64, b: u64) -> u64 { if a >= INF || b >= INF { return INF; } let g = gcd(a, b); std::cmp::min(INF, a * b / g) } #[derive(Clone, Debug)] struct Value { lcm: u64, max: u64, sum: u64, len: u64, } struct Chgcd; impl Effect for Chgcd { type Value = Value; type Query = u64; fn break_cond(x: &Self::Value, q: &Self::Query) -> bool { *q % x.lcm == 0 } fn tag_cond(x: &Self::Value, _: &Self::Query) -> bool { x.max * x.len == x.sum } fn update(x: &mut Self::Value, q: &Self::Query) { let v = gcd(*q, x.max); x.max = v; x.lcm = v; x.sum = v * x.len; } } struct Set; impl Effect for Set { type Value = Value; type Query = u64; fn break_cond(_: &Self::Value, _: &Self::Query) -> bool { false } fn tag_cond(_: &Self::Value, _: &Self::Query) -> bool { true } fn update(x: &mut Self::Value, q: &Self::Query) { x.max = *q; x.lcm = *q; x.sum = *q * x.len; } } struct Yuki880; impl Beats for Yuki880 { type Value = Value; fn merge(l: &Self::Value, r: &Self::Value) -> Self::Value { Value { lcm: lcm(l.lcm, r.lcm), max: std::cmp::max(l.max, r.max), sum: l.sum + r.sum, len: l.len + r.len, } } fn propagate(p: &mut Self::Value, c: &mut [Self::Value]) { let v = p.max; if v * p.len == p.sum { for c in c.iter_mut() { c.lcm = v; c.max = v; c.sum = v * c.len; } } } fn e() -> Self::Value { Value { lcm: 1, max: 0, sum: 0, len: 0, } } } fn run(sc: &mut scanner::Scanner) { let out = std::io::stdout(); let mut out = std::io::BufWriter::new(out.lock()); let n: usize = sc.next(b' '); let q: usize = sc.next(b'\n'); let mut ini = Vec::with_capacity(n); for i in 0..n { let d = if i == n - 1 {b'\n'} else {b' '}; let v: u64 = sc.next(d); ini.push(Value { lcm: v, max: v, sum: v, len: 1, }); } let mut seg = SegmentTreeBeats::::new(&ini); for _ in 0..q { let op: u8 = sc.next(b' '); let (l, r, x) = if op <= 2 { let l: usize = sc.next(b' '); let r: usize = sc.next(b' '); let x: u64 = sc.next(b'\n'); (l - 1, r, x) } else { let l: usize = sc.next(b' '); let r: usize = sc.next(b'\n'); (l - 1, r, 0) }; if op == 1 { seg.update::(l, r, x); } else if op == 2 { seg.update::(l, r, x); } else if op == 3 { let ans = seg.find(l, r).max; writeln!(out, "{}", ans).ok(); } else { let ans = seg.find(l, r).sum; writeln!(out, "{}", ans).ok(); } } }