Cod sursa(job #2932424)

Utilizator valentinchipuc123Valentin Chipuc valentinchipuc123 Data 2 noiembrie 2022 20:39:21
Problema Padure2 Scor 100
Compilator cpp-64 Status done
Runda Arhiva ICPC Marime 1.59 kb
#include <bits/stdc++.h>

using namespace std;

ifstream f("padure2.in");
ofstream g("padure2.out");

const int mod=2000003;
int n,m,c;
int ans=0,jos,sus,fact[2000005];

struct pct
{
    int a,b;
    int ways=0;
} v[1001];

bool cmp(pct a,pct b)
{
    if(a.b==b.b) return a.a<b.a;
    return a.b<b.b;
}

int lgput(int a,int b)
{
    int s=1;

    while(b)
    {
        if(b%2==1)
            s=1LL*s*a%mod;
        a=1LL*a*a%mod;
        b/=2;
    }
    return s;
}

int main()
{
    f>>n>>m;

    f>>c;
    for(int i=1; i<=c; i++)
        f>>v[i].a>>v[i].b;

    fact[0]=1;
    for(int i=1; i<=2*max(n,m); i++)
        fact[i]=1LL*fact[i-1]*i%mod;

    sort(v+1,v+c+1,cmp);

    jos=n+m-2;
    sus=n-1;
    ans=1LL*fact[jos]*lgput(fact[sus],mod-2)%mod*lgput(fact[jos-sus],mod-2)%mod;

    for(int i=1; i<=c; i++)
    {
        jos=v[i].a+v[i].b-2;
        sus=v[i].a-1;
        v[i].ways=1LL*fact[jos]*lgput(fact[sus],mod-2)%mod*lgput(fact[jos-sus],mod-2)%mod;

        for(int j=1; j<i; j++)
        {
            if(v[j].a<=v[i].a&&v[j].b<=v[i].b)
            {
                jos=v[i].a-v[j].a+v[i].b-v[j].b;
                sus=v[i].a-v[j].a;

                v[i].ways=v[i].ways+mod-1LL*fact[jos]*lgput(fact[sus],mod-2)%mod*lgput(fact[jos-sus],mod-2)%mod*v[j].ways%mod;
                v[i].ways%=mod;
            }
        }
        jos=n-v[i].a+m-v[i].b;
        sus=n-v[i].a;
        ans=ans+mod-1LL*fact[jos]*lgput(fact[sus],mod-2)%mod*lgput(fact[jos-sus],mod-2)%mod*v[i].ways%mod;
        ans%=mod;
    }

    g<<ans;

    return 0;
}