結果
| 問題 |
No.1036 Make One With GCD 2
|
| コンテスト | |
| ユーザー |
👑 |
| 提出日時 | 2025-11-14 21:42:21 |
| 言語 | Rust (1.83.0 + proconio) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 9,521 bytes |
| コンパイル時間 | 12,848 ms |
| コンパイル使用メモリ | 395,792 KB |
| 実行使用メモリ | 19,664 KB |
| 最終ジャッジ日時 | 2025-11-14 21:42:50 |
| 合計ジャッジ時間 | 20,327 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 WA * 1 |
| other | AC * 31 WA * 10 |
ソースコード
#[allow(unused_imports)]
use std::{
convert::{Infallible, TryFrom, TryInto as _}, fmt::{self, Debug, Display, Formatter,},
fs::File, hash::{Hash, Hasher, BuildHasherDefault}, 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, btree_map::Range as BTreeRange}, 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},
};
pub trait SortD{ fn sort_d(&mut self); }
impl<T: Ord> SortD for Vec<T>{ fn sort_d(&mut self) {
self.sort_by(|u, v| v.cmp(&u));
} }
pub trait Mx{fn max(&self, rhs: Self)->Self;}
impl Mx for f64{ fn max(&self, rhs: Self)->Self{if *self < rhs{ rhs } else { *self } }}
pub trait Mi{ fn min(&self, rhs: Self)->Self; }
impl Mi for f64{ fn min(&self, rhs: Self)->Self{ if *self > rhs{ rhs } else { *self } } }
pub fn gcd(mut a: i64, mut b: i64)->i64{ if b==0{return a;}(a,b)=(a.abs(),b.abs());while b!=0{ let c = a;a = b;b = c%b; }a }
pub fn factorial_i64(n: usize)->(Vec<i64>, Vec<i64>){ let mut res = vec![1; n+1];let mut inv = vec![1; n+1];for i in 0..n{ res[i+1] = (res[i]*(i+1)as i64)%MOD; }inv[n] = mod_inverse(res[n], MOD);for i in (0..n).rev(){ inv[i] = inv[i+1]*(i+1) as i64%MOD; }(res, inv) }
pub fn floor(a:i64, b:i64)->i64{let res=(a%b+b)%b;(a-res)/b}
pub fn extended_gcd(a:i64,b:i64)->(i64,i64,i64)
{if b==0{(a,1,0)}else{let(g,x,y)=extended_gcd(b,a%b);(g,y,x-floor(a,b)*y)}}
pub fn mod_inverse(a:i64,m:i64)->i64{let(_,x,_) =extended_gcd(a,m);(x%m+m)%m}
pub fn comb(a: i64, b: i64, f: &Vec<(i64, i64)>)->i64{
if a<b{return 0;}else if b==0 || a==b{ return 1; }
else{let x=f[a as usize].0;
let y=f[(a-b) as usize].1;let z=f[b as usize].1;return((x*y)%MOD)*z%MOD;}}
pub fn factorial(x: i64)->Vec<(i64, i64)>{
let mut f=vec![(1i64,1i64),(1, 1)];let mut z = 1i64;
let mut inv = vec![0; x as usize+10];inv[1] = 1;
for i in 2..x+1{z=(z*i)%MOD;
let w=(MOD-inv[(MOD%i)as usize]*(MOD/i)%MOD)%MOD;
inv[i as usize] = w;
f.push((z, (f[i as usize-1].1*w)%MOD));}return f;}
pub fn fast_mod_pow(x: i64,p: usize, m: i64)->i64{
let mut res=1;let mut t=x;let mut z=p;while z > 0{
if z%2==1{res = (res*t)%m;}t = (t*t)%m;z /= 2; }res}
pub fn vec_out<T>(v: Vec<T>) where T: ToString{ println!("{}", v.iter().map(|x| x.to_string()).collect::<Vec<_>>().join(" ")); }
pub fn one_add_out(v: Vec<usize>){println!("{}", v.iter().map(|x| (x+1).to_string()).collect::<Vec<_>>().join(" "));}
pub trait Chmax{
fn chmax(&mut self, rhs: Self);
}
impl Chmax for i64{
fn chmax(&mut self, rhs: Self) {
*self = (*self).max(rhs)
}
}
impl Chmax for i32{
fn chmax(&mut self, rhs: Self) {
*self = (*self).max(rhs)
}
}
impl Chmax for f64{
fn chmax(&mut self, rhs: Self) {
*self = (*self).max(rhs)
}
}
impl Chmax for usize{
fn chmax(&mut self, rhs: Self) {
*self = (*self).max(rhs)
}
}
pub trait Chmin {
fn chmin(&mut self, rhs: Self);
}
impl Chmin for i64{
fn chmin(&mut self, rhs: Self) {
*self = (*self).min(rhs)
}
}
impl Chmin for i32{
fn chmin(&mut self, rhs: Self) {
*self = (*self).min(rhs)
}
}
impl Chmin for f64{
fn chmin(&mut self, rhs: Self) {
*self = (*self).min(rhs)
}
}
impl Chmin for usize{
fn chmin(&mut self, rhs: Self) {
*self = (*self).min(rhs)
}
}
#[derive(Debug, Clone)]
pub struct Counter<T: Ord>{
c: usize,
map: BTreeMap<T, usize>,
}
impl<T: Copy+Ord> Counter<T>{
pub fn new()->Self{
Counter{
c: 0,
map: BTreeMap::new(),
}
}
#[inline(always)]
pub fn range<R>(&self, range: R)->BTreeRange<'_, T, usize> where R: RangeBounds<T>{
self.map.range(range)
}
#[inline(always)]
pub fn mi(&self)->Option<T>{
if let Some((x, _)) = self.range(..).next(){
Some(*x)
} else {
None
}
}
#[inline(always)]
pub fn mx(&self)->Option<T>{
if let Some((x, _)) = self.range(..).next_back(){
Some(*x)
} else {
None
}
}
#[inline(always)]
pub fn one_add(&mut self, x: T){
*self.map.entry(x).or_insert(0) += 1;
self.c += 1;
}
#[inline(always)]
pub fn one_sub(&mut self, x: T){
if !self.map.contains_key(&x){return}
let e = self.map.entry(x).or_insert(0);
*e = e.saturating_sub(1);
if self.map[&x] <= 0{
self.map.remove(&x);
}
self.c = self.c.saturating_sub(1);
}
#[inline(always)]
pub fn one_update(&mut self, x: T, y: T){
self.one_sub(x);
self.one_add(y);
}
#[inline(always)]
pub fn del(&mut self, x: T){
self.c = self.c.saturating_sub(*self.map.get(&x).unwrap_or(&0));
self.map.remove(&x);
}
#[inline(always)]
pub fn add(&mut self, x: T, c: usize){
*self.map.entry(x).or_insert(0) += c;
self.c += c;
}
#[inline(always)]
pub fn sub(&mut self, x: T, c: usize){
let e = self.map.entry(x).or_insert(0);
*e = e.saturating_sub(c);
if self.map[&x] == 0{
self.map.remove(&x);
}
self.c = self.c.saturating_sub(c);
}
#[inline(always)]
pub fn include(&self, x: T)->bool{
self.map.contains_key(&x)
}
#[inline(always)]
pub fn cnt(&self, x: T)->usize{
*self.map.get(&x).unwrap_or(&0)
}
#[inline(always)]
pub fn is_empty(&self)->bool{
self.map.is_empty()
}
#[inline(always)]
pub fn len(&self)->usize{
self.map.len()
}
#[inline(always)]
pub fn clear(&mut self){
self.map.clear();
self.c = 0;
}
#[inline(always)]
pub fn merge(&mut self, rhs: &mut Counter<T>){
if self.len() < rhs.len(){
swap(self, rhs);
}
for (&k, &v) in rhs.map.iter(){
self.add(k, v);
}
rhs.clear();
}
}
#[allow(unused)]
mod fxhash{
use std::hash::BuildHasherDefault;
#[derive(Default)]
pub struct FxHasher{
pub hash: u64,
}
impl std::hash::Hasher for FxHasher{
#[inline(always)]
fn finish(&self) -> u64 {
self.hash
}
#[inline(always)]
fn write(&mut self, bytes: &[u8]) {
let mut h = self.hash;
for &b in bytes{
h = h.rotate_left(5)^(b as u64);
h = h.wrapping_mul(0x517cc1b727220a95);
}
self.hash = h;
}
}
pub type FxBuildHasher = BuildHasherDefault<FxHasher>;
pub type FxMap<K, V> = std::collections::HashMap<K, V, FxBuildHasher>;
pub type FxSet<K> = std::collections::HashSet<K, FxBuildHasher>;
}
#[allow(unused_imports)]
//use fxhash::{FxSet, FxMap, FxBuildHasher};
#[allow(unused_imports)]
use proconio::{input, input_interactive, marker::{*}, fastout};
/*
#[allow(unused_imports)]
use rustc_hash::FxHasher;
#[allow(dead_code)]
type FxMap<K, V> = HashMap<K, V, BuildHasherDefault<FxHasher>>;
#[allow(dead_code)]
type FxSet<K> = HashSet<K, BuildHasherDefault<FxHasher>>;
#[allow(unused_imports)]
use rand::{Rng, seq::SliceRandom, prelude::*};
#[allow(unused_imports)]
use itertools::Itertools;
#[allow(unused_imports)]
use ordered_float::OrderedFloat;
#[allow(unused_imports)]
use num_bigint::BigInt;
#[allow(unused_imports)]
use ac_library::{*, ModInt998244353 as mint};
*/
#[allow(dead_code)]
//type MI = StaticModInt<Mod998244353>;pub fn factorial_mint(n: usize)->(Vec<MI>, Vec<MI>){ let mut res = vec![mint::new(1); n+1];let mut inv = vec![mint::new(1); n+1];for i in 0..n{res[i+1] = res[i]*(i+1);}inv[n] = mint::new(1)/res[n];for i in (0..n).rev(){inv[i] = inv[i+1]*(i+1);}(res, inv)}
#[allow(dead_code)]
const INF: i64 = 1<<60;
#[allow(dead_code)]
const I: i32 = 1<<30;
#[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)];
//#[fastout]
fn main() {
input!{
n: usize,
a: [i64; n],
}
let mut r = 0;
let mut ans = 0;
let mut data = a.clone();
let mut ac = 0;
let mut pre = 0;
for i in 0..n{
if pre==i{
while r < n {
let ng = gcd(ac, a[r]);
if ng!=1{
ac = ng;
r += 1;
} else {
break;
}
}
ac = 0;
if pre+1 < r{
for j in (pre+1..r).rev(){
if j+1 < r{
data[j] = gcd(data[j], data[j+1]);
}
}
}
pre = r;
} else {
let mut g = gcd(ac, data[i]);
while r < n && g != 1{
let ng = gcd(g, a[r]);
if ng!=1{
g = ng;
ac = gcd(ac, a[r]);
r += 1;
} else {
break;
}
}
}
ans += n-r;
r.chmax(i+1);
}
println!("{}", ans);
}