const MOD: usize = 1e9 as usize + 7; fn power(base: usize, times: usize) -> usize { if times == 0 { return 1usize; } if times == 1 { return base; } let temp = power(base, times/2); temp * temp % MOD * power(base, times%2) % MOD } fn main() { let mut hw = String::new(); std::io::stdin().read_line(&mut hw).ok(); let hw: Vec = hw.trim().split_whitespace().map(|s| s.parse().unwrap()).collect(); let h = hw[0]; let w = hw[1]; let grid = (0..h).map(|_| { let mut temp = String::new(); std::io::stdin().read_line(&mut temp).ok(); let temp: Vec = temp.trim().split_whitespace().map(|s| s.parse().unwrap()).collect(); temp }) .collect::>(); let mut q = String::new(); std::io::stdin().read_line(&mut q).ok(); let q: usize = q.trim().parse().unwrap(); let queries = (0..q).map(|_| { let mut temp = String::new(); std::io::stdin().read_line(&mut temp).ok(); let temp: Vec = temp.trim().split_whitespace().map(|s| s.parse().unwrap()).collect(); (temp[0]-1, temp[1]-1) }) .collect::>(); let prod = grid.iter().map(|v| v.iter().filter(|&&vv| vv > 0).fold(1usize, |x, &y| x * y % MOD)).fold(1usize, |x, y| x * y % MOD); let hmap = (0..h).map(|i| (0..w).filter(|&j| grid[i][j] > 0).fold(1usize, |x, y| x * grid[i][y] % MOD)).collect::>(); let wmap = (0..w).map(|j| (0..h).filter(|&i| grid[i][j] > 0).fold(1usize, |x, y| x * grid[y][j] % MOD)).collect::>(); let hzeros = (0..h).map(|i| (0..w).filter(|&j| grid[i][j] == 0).count()).collect::>(); let wzeros = (0..w).map(|j| (0..h).filter(|&i| grid[i][j] == 0).count()).collect::>(); let zero_cnt = hzeros.iter().sum::(); for &(x, y) in queries.iter() { let zero_cnt = zero_cnt + if grid[x][y] == 0 { 1 } else { 0 } - hzeros[x] - wzeros[y]; if zero_cnt > 0 { println!("0"); continue; } let point = grid[x][y].max(1); let result = prod * point % MOD * power(hmap[x], MOD-2) % MOD * power(wmap[y], MOD-2) % MOD; println!("{}", result); } }