码迷,mamicode.com
首页 > 其他好文 > 详细

【codevs1515】 跳

时间:2016-12-25 01:12:12      阅读:251      评论:0      收藏:0      [点我收藏+]

标签:inf   组合   ret   棋盘   har   ring   highlight   print   $$   

http://codevs.cn/problem/1515/ (题目链接)

题意

  给出一个棋盘,规定走到(x,y)的花费C(x,y)=C(x-1,y)+C(x,y-1),x=0或y=0时C(x,y)=1。求从(0,0)走到(n,m)的路径上的花费和的最小值。

Solution

  画画图很容易就会发现是个杨辉三角,我们假设${n<=m}$于是问题就转化为了求解:

$${\sum^{n}_{i=0}{C^i_{m+i}}}$$

  然后就不会了,翻了题解,没想到还有这种公式。。

$${\sum^{n}_{i=0}{C^i_{m+i}}=C^n_{n+m+1}}$$

  来证明一下。对于一个组合${C^n_m}$我们可以这样来表示:格点图上从(0,0)到(n-m,m)的路径条数。(不好画图,尴尬→_→)

  那么${\sum^{n}_{i=0}{C^i_{m+i}}}$就可以分别表示为:从(0,0)到(m,0)的路径条数;从(0,0)到(m,1)的路径条数;从(0,0)到(m,2)的路径条数······从(0,0)到(m,n)的路径条数。

  而${C^n_{n+m+1}}$就可以表示为从(0,0)到(m+1,n)的路径条数。这与前者有什么关系呢?很简单,自己画个图就明白了。

  于是我们用Lucas定理求解组合数取模。(然而这数据范围是什么鬼)

代码

// codevs1515
#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<queue>
#define LL long long
#define inf 2147483640
#define MOD 1000000007
#define Pi acos(-1.0)
#define free(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout);
using namespace std;



LL power(LL a,LL b) {
	LL res=1;
	while (b) {
		if (b&1) res=res*a%MOD;
		b>>=1;a=a*a%MOD;
	}
	return res;
}
LL C(LL n,LL m) {
	LL s1=1,s2=1;
	for (int i=n-m+1;i<=n;i++) s1=s1*i%MOD;
	for (int i=1;i<=m;i++) s2=s2*i%MOD;
	return s1*power(s2,MOD-2)%MOD;
}
LL Lucas(LL n,LL m) {
	if (m==0) return 1;
	return C(n%MOD,m%MOD)*Lucas(n/MOD,m/MOD)%MOD;
}
int main() {
	LL n,m;
	scanf("%lld%lld",&n,&m);
	printf("%lld",(max(n,m)+Lucas(n+m+1,min(n,m)))%MOD);
	return 0;
}

  

【codevs1515】 跳

标签:inf   组合   ret   棋盘   har   ring   highlight   print   $$   

原文地址:http://www.cnblogs.com/MashiroSky/p/6218682.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!