Gerald and Giant Chess

Problem's Link: http://codeforces.com/contest/559/problem/C

Mean:

一个n*m的网格，让你从左上角走到右下角，有一些点不能经过，问你有多少种方法。

analyse:

BZOJ上的原题。

首先把坏点和终点以

　　fi=Cxixi+yi−∑(xj<=xi,yj<=yi)C(xi−xj)(xi−xj)+(yi−yj)∗fj

当于枚举第一个遇到的坏点是啥，然后从总方案里减掉。

然后再利用组合数取模就行。

（http://blog.csdn.net/popoqqq/article/details/46121519）

Time complexity: O(N*N)

Source code:

/*
* this code is made by crazyacking
* Verdict: Accepted
* Submission Date: 2015-07-22-23.43
* Time: 0MAXMS
* MAXMemory: 137KB
*/
#include <queue>
#include <cstdio>
#include <set>
#include <string>
#include <stack>
#include <cmath>
#include <climits>
#include <map>
#include <cstdlib>
#include <iostream>
#include <vector>
#include <algorithm>
#include <cstring>
#define mod 1000000007
#define  Ulong long unsigned long long
using namespace std;

const int MAXN = 200005;
long long h, w, n, ans;
long long q1[MAXN], q2[MAXN], f[MAXN];
struct point
{
      int x, y;
      friend bool operator<( point x1, point y1 )
      {
            if( x1.x != y1.x )return x1.x < y1.x;
            return x1.y < y1.y;
      }
} a[MAXN];

long long quick_pow( long long x, long long y )
{
      long long ans = 1;
      while( y )
      {
            if( y & 1 )ans = ( ans * x ) % mod;
            y >>= 1;
            x = ( x * x ) % mod;
      }
      return ans;
}

int main()
{
      cin >> h >> w >> n;
      for( int i = 1; i <= n; i++ )
            cin >> a[i].x >> a[i].y;
      for( int i = 0; i < MAXN; ++i )
            q1[i] = q2[i] = f[i] = 0;
      q1[0] = q2[0] = 1;
      sort( a + 1, a + n + 1 );
      for( int i = 1; i <= h + w; ++i )
            q1[i] = q1[i - 1] * i % mod;
      q2[h + w] = quick_pow( q1[h + w], mod - 2 );
      for( int i = h + w - 1; i >= 1; i-- )
            q2[i] = q2[i + 1] * ( i + 1 ) % mod;
      ans = q1[h + w - 2] * q2[h - 1] % mod * q2[w - 1] % mod;
      for( int i = 1; i <= n; i++ )
      {
            f[i] = q1[a[i].x + a[i].y - 2] * q2[a[i].x - 1] % mod * q2[a[i].y - 1] % mod;
            for( int j = 1; j < i; j++ )
                  if( a[j].x <= a[i].x && a[j].y <= a[i].y )
                        f[i] -= f[j] * q1[a[i].x + a[i].y - a[j].x - a[j].y] % mod * q2[a[i].x - a[j].x] % mod * q2[a[i].y - a[j].y] % mod;
            f[i] = ( f[i] % mod + mod ) % mod;
            ans = ans - f[i] * q1[h + w - a[i].y - a[i].x] % mod * q2[h - a[i].x] % mod * q2[w - a[i].y] % mod;
      }
      ans = ( ans % mod + mod ) % mod;
      cout << ans << endl;
      return 0;
}

来源：https://www.cnblogs.com/crazyacking/p/4669486.html

标签

include

dp - Codeforces Round #313 (Div. 1) C. Gerald and Giant Chess

Gerald and Giant Chess

Problem's Link: http://codeforces.com/contest/559/problem/C