結果

問題 No.5002 stick xor
ユーザー hatoohatoo
提出日時 2018-05-26 14:26:06
言語 Rust
(1.77.0)
結果
AC  
実行時間 984 ms / 1,000 ms
コード長 9,447 bytes
コンパイル時間 37,951 ms
実行使用メモリ 3,368 KB
スコア 42,417
最終ジャッジ日時 2018-05-26 14:26:46
ジャッジサーバーID
(参考情報)
judge8 /
純コード判定しない問題か言語
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 984 ms
3,364 KB
testcase_01 AC 983 ms
3,364 KB
testcase_02 AC 983 ms
3,364 KB
testcase_03 AC 982 ms
3,368 KB
testcase_04 AC 983 ms
3,364 KB
testcase_05 AC 983 ms
3,360 KB
testcase_06 AC 983 ms
3,364 KB
testcase_07 AC 982 ms
3,360 KB
testcase_08 AC 983 ms
3,368 KB
testcase_09 AC 983 ms
3,364 KB
testcase_10 AC 983 ms
3,364 KB
testcase_11 AC 983 ms
3,364 KB
testcase_12 AC 983 ms
3,368 KB
testcase_13 AC 983 ms
3,364 KB
testcase_14 AC 983 ms
3,368 KB
testcase_15 AC 983 ms
3,364 KB
testcase_16 AC 983 ms
3,368 KB
testcase_17 AC 983 ms
3,364 KB
testcase_18 AC 983 ms
3,360 KB
testcase_19 AC 983 ms
3,368 KB
testcase_20 AC 983 ms
3,364 KB
testcase_21 AC 984 ms
3,360 KB
testcase_22 AC 983 ms
3,364 KB
testcase_23 AC 983 ms
3,368 KB
testcase_24 AC 983 ms
3,368 KB
testcase_25 AC 982 ms
3,360 KB
testcase_26 AC 983 ms
3,364 KB
testcase_27 AC 982 ms
3,360 KB
testcase_28 AC 983 ms
3,368 KB
testcase_29 AC 983 ms
3,368 KB
testcase_30 AC 982 ms
3,368 KB
testcase_31 AC 982 ms
3,364 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
warning: unused variable: `l`
   --> Main.rs:141:43
    |
141 | 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

ソースコード

diff #

/**
*  _           _                 __                            _   _ _   _                                 _                    _                  _
* | |         | |               / /                           | | (_) | (_)                               | |                  (_)                | |
* | |__   __ _| |_ ___   ___   / /__ ___  _ __ ___  _ __   ___| |_ _| |_ ___   _____ ______ _ __ _   _ ___| |_ ______ ___ _ __  _ _ __  _ __   ___| |_ ___
* | '_ \ / _` | __/ _ \ / _ \ / / __/ _ \| '_ ` _ \| '_ \ / _ \ __| | __| \ \ / / _ \______| '__| | | / __| __|______/ __| '_ \| | '_ \| '_ \ / _ \ __/ __|
* | | | | (_| | || (_) | (_) / / (_| (_) | | | | | | |_) |  __/ |_| | |_| |\ 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 = 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), 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 }
    }

    #[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 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);

    let mut ans = vec![(true, 0, 0); k];
    let idx: Vec<usize> = (0..k).collect();
    //idx.sort_by_key(|&i| Rev(ls[i]));

    for i in idx {
        let l = ls[i];
        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);
                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);
                if b > val {
                    tyx = (false, y, x);
                    val = b;
                }
            }
        }

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

    let mut t = 1.0;
    let lim = Duration::from_millis(980);
    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.0 || rng.randf() <= (d / (t * 0.01)).exp() {
            ans[i] = action;
            field.apply(action, l);
        } else {
            field.apply(ans[i], l);
        }

        t *= 0.999999;
    }

    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();
            }
        }
    });
}
0