結果

問題 No.430 文字列検索
コンテスト
ユーザー henoc on demand
提出日時 2019-05-05 22:04:26
言語 Rust
(1.83.0 + proconio)
結果
AC  
実行時間 44 ms / 2,000 ms
コード長 3,481 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 11,161 ms
コンパイル使用メモリ 402,500 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-11-10 00:24:50
合計ジャッジ時間 12,402 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 14
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

fn main() {
    use std::io::Read;
    let mut buf = String::new();
    std::io::stdin().read_to_string(&mut buf).unwrap();
    let mut iter = buf.split_whitespace();
    
    let s: String = iter.next().unwrap().parse::<String>().unwrap();
    let m: usize = iter.next().unwrap().parse::<usize>().unwrap();
    
    let suffix_array = SuffixArray::new(s.as_bytes());
    let mut cnt = 0;
    for _ in 0..m {
        let c: String = iter.next().unwrap().parse::<String>().unwrap();
        cnt += suffix_array.count(c.as_bytes());
    }
    println!("{}", cnt);
}

#[allow(dead_code)]
#[derive(Debug)]
pub struct SuffixArray<'a> {
    sa: Vec<usize>,
    base: &'a [u8],
}
#[allow(dead_code)]
impl<'a> SuffixArray<'a> {
    fn cmp(rank: &[Option<u64>], n: usize, k: usize, i: usize, j: usize) -> std::cmp::Ordering {
        if rank[i] != rank[j] {
            rank[i].cmp(&rank[j])
        } else {
            let ri = if i + k <= n { Some(rank[i + k]) } else { None };
            let rj = if j + k <= n { Some(rank[j + k]) } else { None };
            ri.cmp(&rj)
        }
    }
    pub fn new(base: &'a [u8]) -> SuffixArray {
        let n = base.len();
        let mut sa = vec![0usize; n + 1];
        let mut rank = vec![None; n + 1];
        for i in 0..=n {
            sa[i] = i;
            rank[i] = if i < n { Some(base[i] as u64) } else { None };
        }
        {
            let mut k = 1;
            while k <= n {
                sa.sort_by(|&i, &j| SuffixArray::cmp(&rank, n, k, i, j));
                let mut tmp = vec![0u64; n + 1];
                for i in 1..=n {
                    tmp[sa[i]] = tmp[sa[i - 1]]
                        + (if SuffixArray::cmp(&rank, n, k, sa[i - 1], sa[i])
                            == std::cmp::Ordering::Less
                        {
                            1
                        } else {
                            0
                        });
                }
                for i in 0..=n {
                    rank[i] = Some(tmp[i]);
                }
                k *= 2;
            }
        }
        SuffixArray { sa: sa, base: base }
    }
    pub fn count(&self, key: &[u8]) -> usize {
        self.upper_bound(key) - self.lower_bound(key)
    }
    #[doc = " returns unordered indices"]
    pub fn find(&self, key: &[u8]) -> Vec<usize> {
        let mut ret = vec![];
        for i in self.lower_bound(key)..self.upper_bound(key) {
            ret.push(self.sa[i]);
        }
        ret
    }
    #[doc = " lower bound of rank"]
    pub fn lower_bound(&self, key: &[u8]) -> usize {
        let mut low = 0;
        let mut high = self.base.len() + 1;
        while low != high {
            let mid = (low + high) / 2;
            let z = std::cmp::min(self.sa[mid] + key.len(), self.base.len());
            match self.base[self.sa[mid]..z].cmp(key) {
                std::cmp::Ordering::Less => low = mid + 1,
                _ => high = mid,
            };
        }
        low
    }
    #[doc = " upper bound of rank"]
    pub fn upper_bound(&self, key: &[u8]) -> usize {
        let mut low = 0;
        let mut high = self.base.len() + 1;
        while low != high {
            let mid = (low + high) / 2;
            let z = std::cmp::min(self.sa[mid] + key.len(), self.base.len());
            match self.base[self.sa[mid]..z].cmp(key) {
                std::cmp::Ordering::Greater => high = mid,
                _ => low = mid + 1,
            }
        }
        low
    }
}
0