Pagini recente » Cod sursa (job #665734) | Cod sursa (job #1078086) | Cod sursa (job #1426758) | Cod sursa (job #1408713) | Cod sursa (job #2932424)
#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;
}