#260635 (Rust) No.5002 stick xor

• 問題: No.5002 stick xor
• ユーザー: hatoo
• 提出日時: 2018-05-26 18:10:03
• 言語: Rust (1.77.0 + proconio)
• 結果: AC
• 実行時間: 983 ms / 1,000 ms
• コード長: 11,565 bytes
• コンパイル時間: 35,912 ms
• 実行使用メモリ: 1,320 KB
• スコア: 41,355
• 最終ジャッジ日時: 2018-05-26 18:10:41
• ジャッジサーバーID （参考情報）: judge8 /

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	982 ms; 1,316 KB
testcase_01	AC	982 ms; 1,316 KB
testcase_02	AC	982 ms; 1,316 KB
testcase_03	AC	982 ms; 1,316 KB
testcase_04	AC	982 ms; 1,312 KB
testcase_05	AC	982 ms; 1,316 KB
testcase_06	AC	982 ms; 1,312 KB
testcase_07	AC	982 ms; 1,312 KB
testcase_08	AC	982 ms; 1,312 KB
testcase_09	AC	981 ms; 1,316 KB
testcase_10	AC	982 ms; 1,320 KB
testcase_11	AC	982 ms; 1,312 KB
testcase_12	AC	982 ms; 1,320 KB
testcase_13	AC	982 ms; 1,316 KB
testcase_14	AC	982 ms; 1,316 KB
testcase_15	AC	982 ms; 1,316 KB
testcase_16	AC	982 ms; 1,312 KB
testcase_17	AC	982 ms; 1,316 KB
testcase_18	AC	982 ms; 1,316 KB
testcase_19	AC	982 ms; 1,316 KB
testcase_20	AC	983 ms; 1,316 KB
testcase_21	AC	982 ms; 1,320 KB
testcase_22	AC	982 ms; 1,316 KB
testcase_23	AC	981 ms; 1,316 KB
testcase_24	AC	982 ms; 1,312 KB
testcase_25	AC	982 ms; 1,312 KB
testcase_26	AC	982 ms; 1,316 KB
testcase_27	AC	983 ms; 1,316 KB
testcase_28	AC	982 ms; 1,312 KB
testcase_29	AC	982 ms; 1,316 KB
testcase_30	AC	982 ms; 1,312 KB
testcase_31	AC	982 ms; 1,316 KB

コンパイルメッセージ

warning: function is never used: `score`
   --> Main.rs:141:1
    |
141 | fn score((t, y, x): (bool, usize, usize), n: usize) -> f64 {
    | ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
    |
    = note: #[warn(dead_code)] on by default

warning: method is never used: `zero`
   --> Main.rs:170:5
    |
170 |     fn zero() -> Self {
    |     ^^^^^^^^^^^^^^^^^

warning: method is never used: `or`
   --> Main.rs:196:5
    |
196 |     fn or(&mut self, (t, y, x): (bool, usize, usize), l: usize) {
    |     ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

warning: method is never used: `count_w`
   --> Main.rs:246:5
    |
246 |     fn count_w(&self, (t, y, x): (bool, usize, usize), l: usize) -> i64 {
    |     ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

warning: method is never used: `adj`
   --> Main.rs:250:5
    |
250 |     fn adj(&self, (t, y, x): (bool, usize, usize), l: usize, n: usize) -> f64 {
    |     ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

ソースコード

/**
*  _           _                 __                            _   _ _   _                                 _                    _                  _
* | |         | |               / /                           | | (_) | (_)                               | |                  (_)                | |
* | |__   __ _| |_ ___   ___   / /__ ___  _ __ ___  _ __   ___| |_ _| |_ ___   _____ ______ _ __ _   _ ___| |_ ______ ___ _ __  _ _ __  _ __   ___| |_ ___
* | '_ \ / _` | __/ _ \ / _ \ / / __/ _ \| '_ ` _ \| '_ \ / _ \ __| | __| \ \ / / _ \______| '__| | | / __| __|______/ __| '_ \| | '_ \| '_ \ / _ \ __/ __|
* | | | | (_| | || (_) | (_) / / (_| (_) | | | | | | |_) |  __/ |_| | |_| |\ V /  __/      | |  | |_| \__ \ |_       \__ \ | | | | |_) | |_) |  __/ |_\__ \
* |_| |_|\__,_|\__\___/ \___/_/ \___\___/|_| |_| |_| .__/ \___|\__|_|\__|_| \_/ \___|      |_|   \__,_|___/\__|      |___/_| |_|_| .__/| .__/ \___|\__|___/
*                                                  | |                                                                           | |   | |
*                                                  |_|                                                                           |_|   |_|
*
* https://github.com/hatoo/competitive-rust-snippets
*/
#[allow(unused_imports)]
use std::cmp::{max, min, Ordering};
#[allow(unused_imports)]
use std::collections::{BTreeMap, BTreeSet, BinaryHeap, HashMap, HashSet, VecDeque};
#[allow(unused_imports)]
use std::io::{stdin, stdout, BufWriter, Write};
#[allow(unused_imports)]
use std::iter::FromIterator;
mod util {
    use std::fmt::Debug;
    use std::io::{stdin, stdout, BufWriter, StdoutLock};
    use std::str::FromStr;
    #[allow(dead_code)]
    pub fn line() -> String {
        let mut line: String = String::new();
        stdin().read_line(&mut line).unwrap();
        line.trim().to_string()
    }
    #[allow(dead_code)]
    pub fn chars() -> Vec<char> {
        line().chars().collect()
    }
    #[allow(dead_code)]
    pub fn gets<T: FromStr>() -> Vec<T>
    where
        <T as FromStr>::Err: Debug,
    {
        let mut line: String = String::new();
        stdin().read_line(&mut line).unwrap();
        line.split_whitespace()
            .map(|t| t.parse().unwrap())
            .collect()
    }
    #[allow(dead_code)]
    pub fn with_bufwriter<F: FnOnce(BufWriter<StdoutLock>) -> ()>(f: F) {
        let out = stdout();
        let writer = BufWriter::new(out.lock());
        f(writer)
    }
}
#[allow(unused_macros)]
macro_rules ! get { ( $ t : ty ) => { { let mut line : String = String :: new ( ) ; stdin ( ) . read_line ( & mut line ) . unwrap ( ) ; line . trim ( ) . parse ::<$ t > ( ) . unwrap ( ) } } ; ( $ ( $ t : ty ) ,* ) => { { let mut line : String = String :: new ( ) ; stdin ( ) . read_line ( & mut line ) . unwrap ( ) ; let mut iter = line . split_whitespace ( ) ; ( $ ( iter . next ( ) . unwrap ( ) . parse ::<$ t > ( ) . unwrap ( ) , ) * ) } } ; ( $ t : ty ; $ n : expr ) => { ( 0 ..$ n ) . map ( | _ | get ! ( $ t ) ) . collect ::< Vec < _ >> ( ) } ; ( $ ( $ t : ty ) ,*; $ n : expr ) => { ( 0 ..$ n ) . map ( | _ | get ! ( $ ( $ t ) ,* ) ) . collect ::< Vec < _ >> ( ) } ; ( $ t : ty ;; ) => { { let mut line : String = String :: new ( ) ; stdin ( ) . read_line ( & mut line ) . unwrap ( ) ; line . split_whitespace ( ) . map ( | t | t . parse ::<$ t > ( ) . unwrap ( ) ) . collect ::< Vec < _ >> ( ) } } ; ( $ t : ty ;; $ n : expr ) => { ( 0 ..$ n ) . map ( | _ | get ! ( $ t ;; ) ) . collect ::< Vec < _ >> ( ) } ; }
#[allow(unused_macros)]
macro_rules ! debug { ( $ ( $ a : expr ) ,* ) => { eprintln ! ( concat ! ( $ ( stringify ! ( $ a ) , " = {:?}, " ) ,* ) , $ ( $ a ) ,* ) ; } }
const BIG_STACK_SIZE: bool = false;
#[allow(dead_code)]
fn main() {
    use std::thread;
    if BIG_STACK_SIZE {
        thread::Builder::new()
            .stack_size(32 * 1024 * 1024)
            .name("solve".into())
            .spawn(solve)
            .unwrap()
            .join()
            .unwrap();
    } else {
        solve();
    }
}

#[derive(Eq, PartialEq, Clone, Debug)]
/// Equivalent to std::cmp::Reverse
pub struct Rev<T>(pub T);
impl<T: PartialOrd> PartialOrd for Rev<T> {
    fn partial_cmp(&self, other: &Rev<T>) -> Option<Ordering> {
        other.0.partial_cmp(&self.0)
    }
}
impl<T: Ord> Ord for Rev<T> {
    fn cmp(&self, other: &Rev<T>) -> Ordering {
        other.0.cmp(&self.0)
    }
}

#[derive(Debug)]
#[allow(dead_code)]
pub struct Xorshift {
    seed: u64,
}
impl Xorshift {
    #[allow(dead_code)]
    pub fn new() -> Xorshift {
        Xorshift {
            seed: 0xf0fb588ca2196dac,
        }
    }
    #[allow(dead_code)]
    pub fn with_seed(seed: u64) -> Xorshift {
        Xorshift { seed: seed }
    }
    #[inline(always)]
    #[allow(dead_code)]
    pub fn next(&mut self) -> u64 {
        self.seed = self.seed ^ (self.seed << 13);
        self.seed = self.seed ^ (self.seed >> 7);
        self.seed = self.seed ^ (self.seed << 17);
        self.seed
    }
    #[inline(always)]
    #[allow(dead_code)]
    pub fn rand(&mut self, m: u64) -> u64 {
        self.next() % m
    }
    #[inline(always)]
    #[allow(dead_code)]
    pub fn randf(&mut self) -> f64 {
        use std::mem;
        const UPPER_MASK: u64 = 0x3FF0000000000000;
        const LOWER_MASK: u64 = 0xFFFFFFFFFFFFF;
        let tmp = UPPER_MASK | (self.next() & LOWER_MASK);
        let result: f64 = unsafe { mem::transmute(tmp) };
        result - 1.0
    }
}

fn rand_action(rng: &mut Xorshift, n: usize, l: usize) -> (bool, usize, usize) {
    loop {
        let t = true; //rng.rand(2) == 0;
        let y = rng.rand(n as u64) as usize;
        let x = rng.rand(n as u64) as usize;

        if (t && x + l <= n) || (!t && y + l <= n) {
            break (t, y, x);
        }
    }
}

fn score((t, y, x): (bool, usize, usize), n: usize) -> f64 {
    let s = if t { min(y, n - y) } else { min(x, n - x) };

    0.01 * s as f64 / n as f64
}

struct Field {
    f1: [u64; 60],
    f2: [u64; 60],
}

impl Field {
    fn new(init: &Vec<Vec<bool>>) -> Self {
        let n = init.len();
        let mut f1 = [0; 60];
        let mut f2 = [0; 60];

        for y in 0..n {
            for x in 0..n {
                if init[y][x] {
                    f1[y] |= 1 << x;
                    f2[x] |= 1 << y;
                }
            }
        }

        Self { f1, f2 }
    }

    fn zero() -> Self {
        Self {
            f1: [0; 60],
            f2: [0; 60],
        }
    }

    #[inline]
    fn apply(&mut self, (t, y, x): (bool, usize, usize), l: usize) {
        let w = (1 << l) - 1;
        if t {
            self.f1[y] ^= w << x;

            for x in x..x + l {
                self.f2[x] ^= 1 << y;
            }
        } else {
            self.f2[x] ^= w << y;

            for y in y..y + l {
                self.f1[y] ^= 1 << x;
            }
        }
    }

    #[inline]
    fn or(&mut self, (t, y, x): (bool, usize, usize), l: usize) {
        let w = (1 << l) - 1;
        if t {
            self.f1[y] |= w << x;

            for x in x..x + l {
                self.f2[x] |= 1 << y;
            }
        } else {
            self.f2[x] |= w << y;

            for y in y..y + l {
                self.f1[y] |= 1 << x;
            }
        }
    }

    fn bits(&self, (t, y, x): (bool, usize, usize), l: usize) -> u64 {
        let w = (1 << l) - 1;
        if t {
            (self.f1[y] >> x) & w
        } else {
            (self.f2[x] >> y) & w
        }
    }

    fn score(&self, (t, y, x): (bool, usize, usize), l: usize) -> f64 {
        let b = self.bits((t, y, x), l);

        let mut res = 0.0;

        for i in 0..l {
            if (b >> i) & 1 != (b >> (i + 1)) & 1 {
                res -= 1.0;
            }
        }

        res
    }

    #[inline]
    fn count_b(&self, (t, y, x): (bool, usize, usize), l: usize) -> i64 {
        let w = (1 << l) - 1;
        if t {
            (self.f1[y] & (w << x)).count_ones() as i64
        } else {
            (self.f2[x] & (w << y)).count_ones() as i64
        }
    }
    #[inline]
    fn count_w(&self, (t, y, x): (bool, usize, usize), l: usize) -> i64 {
        l as i64 - self.count_b((t, y, x), l)
    }

    fn adj(&self, (t, y, x): (bool, usize, usize), l: usize, n: usize) -> f64 {
        let mut res = 0.0;
        if t {
            if y > 0 {
                res += self.count_w((t, y - 1, x), l) as f64 / l as f64;
            }

            if y + 1 < n {
                res += self.count_w((t, y + 1, x), l) as f64 / l as f64;
            }
        } else {
            if x > 0 {
                res += self.count_w((t, y, x - 1), l) as f64 / l as f64;
            }

            if x + 1 < n {
                res += self.count_w((t, y, x + 1), l) as f64 / l as f64;
            }
        }

        res
    }
}

use std::time::{Duration, Instant};

fn solve() {
    let start = Instant::now();
    let lim = Duration::from_millis(980);
    let (n, k) = get!(usize, usize);
    let ls = get!(usize;;);
    let field: Vec<Vec<bool>> = (0..n)
        .map(|_| util::line().chars().map(|c| c == '1').collect())
        .collect();

    let mut rng = Xorshift::new();
    let mut field = Field::new(&field);

    let mut ans = vec![(true, 0, 0); k];
    let idx: Vec<usize> = (0..k).collect();

    for i in idx {
        let l = ls[i];
        let mut val = -1000.0;
        let mut tyx = (true, 0, 0);
        for y in 0..n {
            for x in 0..n - l {
                let a = (true, y, x);
                let b = (2 * field.count_b(a, l) - l as i64) as f64 + field.score(a, l) * 0.01;
                if b > val {
                    val = b;
                    tyx = a;
                }
            }
        }
        for y in 0..n - l {
            for x in 0..n {
                let a = (false, y, x);
                let b = (2 * field.count_b(a, l) - l as i64) as f64 + field.score(a, l) * 0.01;
                if b > val {
                    val = b;
                    tyx = a;
                }
            }
        }

        field.apply(tyx, l);
        ans[i] = tyx;
        /*
        let tyx = rand_action(&mut rng, n, l);
        field.apply(tyx, l);
        ans.push(tyx);
        */
    }

    let mut total = 0;
    let mut best = 0;
    let mut best_ans = ans.clone();
    let mut t = 1.0;
    for c in 0.. {
        if c % 1000 == 0 && Instant::now() - start >= lim {
            debug!(c);
            break;
        }

        let i = rng.rand(k as u64) as usize;

        let l = ls[i];

        field.apply(ans[i], l);
        let before = field.count_b(ans[i], l);
        let action = rand_action(&mut rng, n, l);
        let after = field.count_b(action, l);
        let d = after - before;

        if d > 0 || rng.randf() <= (d as f64 / (t * 0.001)).exp() {
            ans[i] = action;
            field.apply(action, l);
            total -= 2 * before - l as i64;
            total += 2 * after - l as i64;
        } else {
            field.apply(ans[i], l);
        }

        if total > best {
            best = total;
            best_ans = ans.clone();
        }

        t *= 0.99999999;
    }

    debug!(t);
    debug!(total);
    debug!(best);

    util::with_bufwriter(|mut out| {
        for (i, (t, y, x)) in best_ans.into_iter().enumerate() {
            if t {
                writeln!(out, "{} {} {} {}", y + 1, x + 1, y + 1, x + ls[i]).unwrap();
            } else {
                writeln!(out, "{} {} {} {}", y + 1, x + 1, y + ls[i], x + 1).unwrap();
            }
        }
    });
}