結果
問題 | No.1202 お菓子の食べ方 |
ユーザー |
![]() |
提出日時 | 2020-08-22 13:39:00 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
TLE
|
実行時間 | - |
コード長 | 3,548 bytes |
コンパイル時間 | 3,064 ms |
コンパイル使用メモリ | 204,776 KB |
最終ジャッジ日時 | 2025-01-13 07:52:33 |
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 2 RE * 18 TLE * 25 |
ソースコード
#include <iostream>#include <random>#include <bits/stdc++.h>using namespace std;using ull = __int128;using ll = long long;#define rep(i,a,b) for(int i=a;i<b;i++)#include <iostream>using namespace std;const int MAX = 2000010;const int mod = 1000000007;const int MOD = 1000000007;long long fac[MAX], finv[MAX], inv[MAX];// テーブルを作る前処理void COMinit() {fac[0] = fac[1] = 1;finv[0] = finv[1] = 1;inv[1] = 1;for (int i = 2; i < MAX; i++){fac[i] = fac[i - 1] * i % MOD;inv[i] = MOD - inv[MOD%i] * (MOD / i) % MOD;finv[i] = finv[i - 1] * inv[i] % MOD;}}// 二項係数計算long long COM(int n, int k){if (n < k) return 0;if (n < 0 || k < 0) return 0;return fac[n] * (finv[k] * finv[n - k] % MOD) % MOD;}class mint {long long x;public:mint(long long x=0) : x((x%mod+mod)%mod) {}mint operator-() const {return mint(-x);}mint& operator+=(const mint& a) {if ((x += a.x) >= mod) x -= mod;return *this;}mint& operator-=(const mint& a) {if ((x += mod-a.x) >= mod) x -= mod;return *this;}mint& operator*=(const mint& a) {(x *= a.x) %= mod;return *this;}mint operator+(const mint& a) const {mint res(*this);return res+=a;}mint operator-(const mint& a) const {mint res(*this);return res-=a;}mint operator*(const mint& a) const {mint res(*this);return res*=a;}mint pow(ll t) const {if (!t) return 1;mint a = pow(t>>1);a *= a;if (t&1) a *= *this;return a;}// for prime modmint inv() const {return pow(mod-2);}mint& operator/=(const mint& a) {return (*this) *= a.inv();}mint operator/(const mint& a) const {mint res(*this);return res/=a;}friend ostream& operator<<(ostream& os, const mint& m){os << m.x;return os;}friend istream& operator>>(istream& lhs,mint& rhs) noexcept {lhs >> rhs.x;return lhs;}};int main() {// 前処理COMinit();int n,m;cin>>n>>m;vector<vector<ll>> o(n, vector<ll>(m));assert(1<=n&&n<=1000&&1<=m&&m<=1000);for (int i = 0; i < n; i++) {for (int j = 0; j < m; j++) {cin >> o.at(i).at(j);assert(1<=o.at(i).at(j)&&o.at(i).at(j)<=1000000);}}mint ans=0;mint p=0;for(int i=1;i<n;i++){for(int j=1;j<m;j++){assert(o.at(i).at(j)>=o.at(i).at(j-1)&&o.at(i).at(j)>=o.at(i-1).at(j));}}for(int i=1;i<n;i++){assert(o.at(i).at(0)>=o.at(i-1).at(0));}for(int j=1;j<m;j++){assert(o.at(0).at(j)>=o.at(0).at(j-1));}for(int i=0;i<n;i++){for(int j=0;j<m;j++){p=0;p=p+COM(o.at(i).at(j)+i+j-1,o.at(i).at(j)-1);if(j!=m-1){p=p+(mint(((o.at(i).at(j)*COM(o.at(i).at(j)+i+j,o.at(i).at(j))-(o.at(i).at(j+1)*COM(o.at(i).at(j+1)+i+j,o.at(i).at(j+1))))))/mint((i+j+1)));}else{p=p+mint((COM(o.at(i).at(j)+i+j,o.at(i).at(j))*o.at(i).at(j)))/mint((i+j+1));}if(i!=n-1){p=p+(mint(((o.at(i).at(j)*COM(o.at(i).at(j)+i+j,o.at(i).at(j))-(o.at(i+1).at(j)*COM(o.at(i+1).at(j)+i+j,o.at(i+1).at(j))))))/mint(i+j+1));}else{p=p+mint((COM(o.at(i).at(j)+i+j,o.at(i).at(j))*o.at(i).at(j)))/mint(i+j+1);}ans=ans+p*COM(i+j,i);}}cout<<ans<<endl;}