#260528 (Rust) No.5002 stick xor

• 問題: No.5002 stick xor
• ユーザー: hatoo
• 提出日時: 2018-05-26 12:05:45
• 言語: Rust (1.77.0 + proconio)
• 結果: AC
• 実行時間: 974 ms / 1,000 ms
• コード長: 10,793 bytes
• コンパイル時間: 39,499 ms
• 実行使用メモリ: 1,364 KB
• スコア: 42,319
• 最終ジャッジ日時: 2018-05-26 12:06:27
• ジャッジサーバーID （参考情報）: judge9 /

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	972 ms; 1,364 KB
testcase_01	AC	972 ms; 1,360 KB
testcase_02	AC	973 ms; 1,360 KB
testcase_03	AC	972 ms; 1,364 KB
testcase_04	AC	973 ms; 1,364 KB
testcase_05	AC	972 ms; 1,364 KB
testcase_06	AC	972 ms; 1,356 KB
testcase_07	AC	973 ms; 1,356 KB
testcase_08	AC	973 ms; 1,364 KB
testcase_09	AC	972 ms; 1,356 KB
testcase_10	AC	972 ms; 1,360 KB
testcase_11	AC	972 ms; 1,364 KB
testcase_12	AC	972 ms; 1,360 KB
testcase_13	AC	973 ms; 1,364 KB
testcase_14	AC	973 ms; 1,360 KB
testcase_15	AC	973 ms; 1,360 KB
testcase_16	AC	973 ms; 1,364 KB
testcase_17	AC	973 ms; 1,364 KB
testcase_18	AC	973 ms; 1,364 KB
testcase_19	AC	973 ms; 1,364 KB
testcase_20	AC	972 ms; 1,364 KB
testcase_21	AC	972 ms; 1,360 KB
testcase_22	AC	973 ms; 1,356 KB
testcase_23	AC	973 ms; 1,360 KB
testcase_24	AC	971 ms; 1,360 KB
testcase_25	AC	973 ms; 1,360 KB
testcase_26	AC	973 ms; 1,360 KB
testcase_27	AC	972 ms; 1,360 KB
testcase_28	AC	974 ms; 1,360 KB
testcase_29	AC	972 ms; 1,356 KB
testcase_30	AC	972 ms; 1,364 KB
testcase_31	AC	973 ms; 1,360 KB

コンパイルメッセージ

warning: unused variable: `l`
   --> Main.rs:157:43
    |
157 | fn score((t, y, x): (bool, usize, usize), l: usize, n: usize) -> f64 {
    |                                           ^
    |
    = note: #[warn(unused_variables)] on by default
    = note: to avoid this warning, consider using `_l` instead

warning: unused variable: `k`
   --> Main.rs:217:13
    |
217 |     let (n, k) = get!(usize, usize);
    |             ^
    |
    = note: to avoid this warning, consider using `_k` instead

warning: function is never used: `apply`
   --> Main.rs:127:1
    |
127 | fn apply((t, y, x): (bool, usize, usize), l: usize, field: &mut Vec<Vec<bool>>) {
    | ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
    |
    = note: #[warn(dead_code)] on by default

warning: function is never used: `count_b`
   --> Main.rs:139:1
    |
139 | fn count_b((t, y, x): (bool, usize, usize), l: usize, field: &Vec<Vec<bool>>) -> 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(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 = 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 apply((t, y, x): (bool, usize, usize), l: usize, field: &mut Vec<Vec<bool>>) {
    if t {
        for x in x..x + l {
            field[y][x] = !field[y][x];
        }
    } else {
        for y in y..y + l {
            field[y][x] = !field[y][x];
        }
    }
}

fn count_b((t, y, x): (bool, usize, usize), l: usize, field: &Vec<Vec<bool>>) -> f64 {
    let mut b = 0.0;
    if t {
        for x in x..x + l {
            if field[y][x] {
                b += 1.0;
            }
        }
    } else {
        for y in y..y + l {
            if field[y][x] {
                b += 1.0;
            }
        }
    }
    b
}

fn score((t, y, x): (bool, usize, usize), l: 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>>) -> Field {
        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;
                }
            }
        }

        Field { f1, f2 }
    }

    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;
            }
        }
    }

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

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

fn solve() {
    let start = Instant::now();
    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);

    // Find baseline with simple greedy

    let mut ans = Vec::new();
    for &l in &ls {
        let mut val = -1000.0;
        let mut tyx = (false, 0, 0);
        for y in 0..n {
            for x in 0..n - l {
                let b = field.count_b((true, y, x), l) + score((true, y, x), l, n); //count_b((true, y, x), l, &field);
                if b > val {
                    tyx = (true, y, x);
                    val = b;
                }
            }
        }
        for y in 0..n - l {
            for x in 0..n {
                let b = field.count_b((false, y, x), l) + score((false, y, x), l, n); // count_b((false, y, x), l, &field);
                if b > val {
                    tyx = (false, y, x);
                    val = b;
                }
            }
        }

        // apply(tyx, l, &mut field);
        field.apply(tyx, l);
        ans.push(tyx);
        /*
        let mut tyx = loop {
            let t = 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);
            }
        };
        //tyx = (true, 0, 0);
        field.apply(tyx, l);
        ans.push(tyx);
        */
    }

    let mut t = 1.0;
    for c in 0.. {
        if c % 10 == 0 && Instant::now() - start >= Duration::from_millis(970) {
            debug!(c);
            break;
        }

        let i = rng.rand(ans.len() 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.0 || rng.randf() <= (d / (t * 0.01)).exp() {
            ans[i] = action;
            field.apply(action, l);
        } else {
            field.apply(ans[i], l);
        }

        t *= 0.999995;

        /*
        let mut val = -1000.0;
        for y in 0..n {
            for x in 0..n {
                if x + l < n {
                    let b = field.count_b((true, y, x), l); // count_b((true, y, x), l, &field) + score((true, y, x), l, n);
                    if b > val {
                        ans[i] = (true, y, x);
                        val = b;
                    }
                }
                if y + l < n {
                    let b = field.count_b((false, y, x), l); // count_b((false, y, x), l, &field) + score((false, y, x), l, n);
                    if b > val {
                        ans[i] = (false, y, x);
                        val = b;
                    }
                }
            }
        }
        // apply(ans[i], l, &mut field);
        field.apply(ans[i], l);
        */
    }

    debug!(t);

    util::with_bufwriter(|mut out| {
        for (i, (t, y, x)) in 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();
            }
        }
    });
}