| # | Time | Username | Problem | Language | Result | Execution time | Memory | 
|---|---|---|---|---|---|---|---|
| 160054 | RedNextCentury | Mobitel (COCI19_mobitel) | C++14 | 1771 ms | 65540 KiB | 
This submission is migrated from previous version of oj.uz, which used different machine for grading. This submission may have different result if resubmitted.
#include<bits/stdc++.h>
using namespace std;
#define ll long long
ll mod = 1000000007;
ll Power(ll a,ll b) {
    if (b==0)return 1;
    ll x = Power(a,b/2);
    x = (x*x)%mod;
    if (b%2)x=(x*a)%mod;
    return x;
}
const int problemsettersIQ = 404; // not found
vector<pair<int,int> > dp[2][problemsettersIQ];
int a[problemsettersIQ][problemsettersIQ];
ll inv(ll x) {return Power(x,mod-2);}
int main()
{
    //ios_base::sync_with_stdio(0);cin.tie(0);cout.tie(0);
    int n,m,k;
    scanf("%d%d%d",&n,&m,&k);
    for (int i=1;i<=n;i++){
        for (int j=1;j<=m;j++) {
            scanf("%d",&a[i][j]);
        }
    }
    if (a[1][1]>k){
        printf("0\n");
        return 0;
    }
    dp[1][1].emplace_back(a[1][1],1);
    for (int i=1;i<=n;i++) {
        for (int j=1;j<=m;j++) {
            if (i==1 && j==1)continue;
            if (i==1){
                for (auto x:dp[i%2][j-1]){
                    if (k/x.first < a[i][j])
                        break;
                    dp[i%2][j].emplace_back(x.first*a[i][j],x.second);
                }
            }
            else if (j==1){
                for (auto x:dp[(i+1)%2][j]){
                    if (k/x.first < a[i][j])
                        break;
                    dp[i%2][j].emplace_back(x.first*a[i][j],x.second);
                }
            }
            else {
                int l=0,r=0;
                vector<pair<int,int> >& L = dp[(i+1)%2][j];
                vector<pair<int,int> >& R = dp[i%2][j-1];
                while(l<L.size() && r<R.size()) {
                    if (L[l].first==R[r].first) {
                        if (k/L[l].first < a[i][j])break;
                        dp[i%2][j].emplace_back(L[l].first*a[i][j],(L[l].second+R[r].second)%mod);
                        l++,r++;
                    } else if (L[l].first<R[r].first) {
                        if (k/L[l].first < a[i][j])break;
                        dp[i%2][j].emplace_back(L[l].first*a[i][j],L[l].second);
                        l++;
                    } else {
                        if (k/R[r].first < a[i][j])break;
                        dp[i%2][j].emplace_back(R[r].first*a[i][j],R[r].second);
                        r++;
                    }
                }
                while(l<L.size()) {
                    if (k/L[l].first < a[i][j])break;
                    dp[i%2][j].emplace_back(L[l].first*a[i][j],L[l].second);
                    l++;
                }
                while(r<R.size()) {
                    if (k/R[r].first < a[i][j])break;
                    dp[i%2][j].emplace_back(R[r].first*a[i][j],R[r].second);
                    r++;
                }
            }
        }
        for (int j=1;j<=m;j++) dp[(i+1)%2][j].clear();
    }
    ll ret=1;
    for (int i=n+m-2 - (n-1) +1;i<=n+m-2;i++)ret=(ret*i)%mod;
    for (int i=1;i<=(n-1);i++)ret=(ret*inv(i))%mod;
    for (auto x:dp[n%2][m]) {
        if (x.first==k)break;
        ret-=x.second;
        ret+=mod;
        ret%=mod;
    }
    printf("%lld\n",ret);
}
Compilation message (stderr)
| # | Verdict | Execution time | Memory | Grader output | 
|---|---|---|---|---|
| Fetching results... | ||||
