CodeVS1515跳

最新推荐文章于 2021-09-17 22:27:09 发布

*ACoder*

最新推荐文章于 2021-09-17 22:27:09 发布

阅读量503

点赞数

分类专栏： # 排列组合

本文链接：https://blog.csdn.net/FSAHFGSADHSAKNDAS/article/details/58607304

版权

排列组合专栏收录该内容

54 篇文章 0 订阅

订阅专栏

链接

　　http://codevs.cn/problem/1515/

题解

　　首先这是一张组合数表， $f(i,j)=C_{i+j}^i$
　　设 $m\le n$ ，你会发现最优策略总是从 $(0,0)$ 走到 $(n,0)$ ，然后再去 $(n,m)$ ，上述都是说的走直线。
　　将这张表用组合数的形式写出来，会发现
　　

a n s = n + C 0 n + C 1 n + 1 + C 2 n + 2 + . . . + C m n + m = n + C n n + C n n + 1 + C n n + 2 + . . . + C n n + m = n + C n + 1 n + 1 + C n n + 1 + C n n + 2 + . . . + C n n + m = n + C n + 1 n + 2 + C n n + 2 + . . . C n n + m = n + C n + 1 n + 3 + C n n + 3 + . . . + C n n + m = . . . = n + C n + 1 n + m + 1 = n + C m n + m + 1

$ans=n+C_n^0+C_{n+1}^1+C_{n+2}^2+...+C_{n+m}^m\\=n+C_n^n+C_{n+1}^n+C_{n+2}^n+...+C_{n+m}^n\\=n+C_{n+1}^{n+1}+C_{n+1}^n+C_{n+2}^n+...+C_{n+m}^n\\=n+C_{n+2}^{n+1}+C_{n+2}^n+...C_{n+m}^n\\=n+C_{n+3}^{n+1}+C_{n+3}^n+...+C_{n+m}^n\\=...\\=n+C_{n+m+1}^{n+1}\\=n+C_{n+m+1}^m$
　　因此对于这道题目
　　

a n s = m a x (n, m) + C m i n (n, m) n + m + 1

$ans=max(n,m)+C_{n+m+1}^{min(n,m)}$
　　网上老用什么lacus什么的，但完全没有必要…，题目中有个条件

nm≤1012 $nm\le 10^{12}$ ，那么

min(n,m)≤106 $min(n,m)\le 10^6$ ，使用高中数学中所说的“组合数的计算式”，是可以

O(min(n,m)) $O(min(n,m))$ 求出的。

代码

//组合数 
#include <cstdio>
#include <algorithm>
#define ll long long
#define maxn 3000000
#define mod 1000000007
using namespace std;
ll n, m;
inline void exgcd(ll a, ll b, ll &x, ll &y)
{
    if(!b){x=1;y=0;return;}
    ll xx, yy;
    exgcd(b,a%b,xx,yy);
    x=yy, y=xx-a/b*yy;
}
ll inv(ll a, ll p)
{
    ll x, y;
    exgcd(a,p,x,y);
    return (x+p)%p;
}
ll C(ll n, ll m)
{
    ll ans=1, i, t;
    for(i=n;i>=n-m+1;i--)ans=ans*(i%mod)%mod;
    for(t=1,i=1;i<=m;i++)t=t*(i%mod)%mod;
    ans=ans*inv(t,mod)%mod;
    return ans;
}
ll calc()
{
    ll ans;
    ans=max(n,m)+C(n+m+1,min(n,m));
    return ans%mod;
}
int main()
{
    scanf("%lld%lld",&n,&m);
    printf("%lld",calc());
    return 0;
}