結果

問題 No.3129 Multiple of Twin Subarray
コンテスト
ユーザー 👑 ArcAki
提出日時 2025-04-25 22:10:23
言語 Rust
(1.83.0 + proconio)
結果
AC  
実行時間 177 ms / 2,000 ms
コード長 5,606 bytes
コンパイル時間 13,804 ms
コンパイル使用メモリ 394,900 KB
実行使用メモリ 37,708 KB
最終ジャッジ日時 2025-04-25 22:11:23
合計ジャッジ時間 18,963 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 46
権限があれば一括ダウンロードができます

ソースコード

diff # 

#[allow(unused_imports)]
use std::{
    convert::{Infallible, TryFrom, TryInto as _}, fmt::{self, Debug, Display, Formatter,},
    fs::{File}, hash::{Hash, Hasher}, iter::{Product, Sum}, marker::PhantomData,
    ops::{Add, AddAssign, Sub, SubAssign, Div, DivAssign, Mul, MulAssign, Neg, RangeBounds},
    str::FromStr, sync::{atomic::{self, AtomicU32, AtomicU64}, Once},
    collections::{*, btree_set::Range}, mem::{swap},
    cmp::{self, Reverse, Ordering, Eq, PartialEq, PartialOrd},
    thread::LocalKey, f64::consts::PI, time::Instant, cell::RefCell,
    io::{self, stdin, Read, read_to_string, BufWriter, BufReader, stdout, Write},
};
#[allow(unused_imports)]
use proconio::{input, input_interactive, marker::{*}};
#[allow(unused_imports)]
//use rand::{thread_rng, Rng, seq::SliceRandom};
#[allow(unused_imports)]
//use ac_library::{*};

#[allow(dead_code)]
const INF: i64 = 1<<62;
#[allow(dead_code)]
const MOD: i64 = 998244353;
#[allow(dead_code)]
const D: [(usize, usize); 4] = [(1, 0), (0, 1), (!0, 0), (0, !0)];
#[allow(dead_code)]
const D2: [(usize, usize); 8] = [(1, 0), (1, 1), (0, 1), (!0, 1), (!0, 0), (!0, !0), (0, !0), (1, !0)];

pub trait SegTreeMonoid{
    type S: Clone;
    fn identity()->Self::S;
    fn op(a: &Self::S, b: &Self::S)->Self::S;
}

pub struct SegTree<M: SegTreeMonoid>{
    n: usize,
    data: Vec<M::S>,
}

impl<M: SegTreeMonoid> SegTree<M> {
    pub fn new(n: usize)->Self{
        let n = n.next_power_of_two();
        let data = vec![M::identity(); 2*n];
        SegTree{
            n, data
        }
    }

    pub fn set(&mut self, i: usize, x: M::S){
        let mut p = i+self.n;
        self.data[p] = x;
        while p > 0{
            p /= 2;
            self.data[p] = M::op(&self.data[p<<1], &self.data[(p<<1)|1]);
        }
    }

    pub fn push(&mut self, i: usize, x: M::S){
        let mut p = i+self.n;
        self.data[p] = M::op(&self.data[p], &x);
        while p > 0{
            p /= 2;
            self.data[p] = M::op(&self.data[p<<1], &self.data[(p<<1)|1]);
        }
    }

    pub fn prod(&mut self, l: usize, r: usize)->M::S{
        let mut p_l = l+self.n;
        let mut p_r = r+self.n;
        let mut res_l = M::identity();
        let mut res_r = M::identity();
        while p_l < p_r{
            if p_l&1==1{
                res_l = M::op(&res_l, &self.data[p_l]);
                p_l += 1;
            }
            if p_r&1==1{
                p_r -= 1;
                res_r = M::op(&self.data[p_r], &res_r);
            }
            p_l >>= 1;
            p_r >>= 1;
        }
        M::op(&res_l, &res_r)
    }

    pub fn all_prod(&mut self)->M::S{
        self.data[1].clone()
    }

    pub fn max_right<F>(&self, mut l: usize, f: F)->usize where F: Fn(&M::S)->bool{
        assert!(f(&M::identity())); // これはバグってくれないと多分デバックが悲惨
        if l==self.n{return self.n}
        l += self.n; let mut ac = M::identity();
        while {
            while l%2==0{
                l >>= 1;
            }
            if !f(&M::op(&ac, &self.data[l])){
                while l < self.n{
                    l <<= 1;
                    let res = M::op(&ac, &self.data[l]);
                    if f(&res){
                        ac = res;
                        l += 1;
                    }
                }
                return l-self.n;
            }
            ac = M::op(&ac, &self.data[l]);
            l += 1;
            let z = l as isize;
            (z & -z) != z
        } {}
        self.n
    }

    pub fn min_left<F>(&self, mut r: usize, f: F)->usize where F: Fn(&M::S)->bool{
        assert!(f(&M::identity()));r += self.n;let mut ac = M::identity();
        if r==0{return 0}
        while {
            while r > 1 && r%2 == 1{
                r >>= 1;
            }
            if !f(&M::op(&ac, &self.data[r])){
                while r < self.n{
                    r = 2*r+1;
                    let res = M::op(&ac, &self.data[r]);
                    if f(&res){
                        ac = res;
                        r -= 1;
                    }
                }
                return r+1-self.n;
            }
            ac = M::op(&ac, &self.data[r]);
            let z = r as isize;
            z & -z == z
        } {}
        0
    }
}

struct M1;
impl SegTreeMonoid for M1{
    type S = (i64, i64, i64, i64);

    fn identity() -> Self::S {
        (-INF, -INF, -INF, -INF)
    }

    fn op(&a: &Self::S, &b: &Self::S) -> Self::S {
        if b.0==-INF{
            a
        } else {
            let l = a.1.max(a.0+b.1);
            let r = b.2.max(a.2+b.0);
            (a.0+b.0, l, r, a.3.max(b.3).max(a.2+b.1).max(l).max(r))
        }
    }
}

struct M2;
impl SegTreeMonoid for M2 {
    type S = (i64, i64, i64, i64);

    fn identity() -> Self::S {
        (INF, INF, INF, INF)
    }

    fn op(&a: &Self::S, &b: &Self::S) -> Self::S {
        if b.0==INF{
            a
        } else {
            let l = a.1.min(a.0+b.1);
            let r = b.2.min(a.2+b.0);
            (a.0+b.0, l, r, a.3.min(b.3).min(a.2+b.1).min(l).min(r))
        }
    }
}

//use proconio::fastout;
//#[fastout]
fn main() {
    input! {
        n: usize,
        a: [i64; n],
    }
    let mut seg1 = SegTree::<M1>::new(n);
    let mut seg2 = SegTree::<M2>::new(n);
    for (i, &v) in a.iter().enumerate(){
        seg1.set(i, (v, v, v, v));
        seg2.set(i, (v, v, v, v));
    }
    let mut ans = -INF;
    for i in 1..n{
        ans = ans.max(seg1.prod(0,i).3*seg1.prod(i, n).3).max(seg2.prod(0, i).3*seg2.prod(i, n).3);
    }
    println!("{}", ans);
}
0