hdu5698百度之星瞬间移动(数论/卢卡斯定理)
2016-05-23 21:45
323 查看
这道题我拿道题是先写了个4维dp打表找规律
memset(dp,0,sizeof(dp));
int n,i,j;
int l,m;
dp[1][1]=1;
for(i=2;i<=20;i++)
{
for(j=2;j<=20;j++)
{
for(l=1;l<i;l++)
{
for(m=1;m<j;m++)
{
dp[i][j]+=dp[l][m];
}
}
}
}
然后发现规律就是一个斜着的杨辉三角,杨辉三角的公式就是 第 n 行第i个数就是 (a+b)的n次方的第i个多项式的系数。
然后这里就要用到卢卡斯定理了,套了个模板就好。
#include <iostream>
#include <string.h>
#include <stdio.h>
using namespace std;
typedef __int64 LL;
LL n,m,p;
LL quick_mod(LL a, LL b)
{
LL ans = 1;
a %= p;
while(b)
{
if(b & 1)
{
ans = ans * a % p;
b--;
}
b >>= 1;
a = a * a % p;
}
return ans;
}
LL C(LL n, LL m)
{
if(m > n) return 0;
LL ans = 1;
for(int i=1; i<=m; i++)
{
LL a = (n + i - m) % p;
LL b = i % p;
ans = ans * (a * quick_mod(b, p-2) % p) % p;
}
return ans;
}
LL Lucas(LL n, LL m)
{
if(m == 0) return 1;
return C(n % p, m % p) * Lucas(n / p, m / p) % p;
}
int main()
{
int T;
p=1000000007;
while(~scanf("%I64d%I64d", &n, &m))
{
LL maxi,mini;
maxi=n+m-3;
// if(n==m)maxi+=1;
mini=m-1;
// printf("%I64d %I64d\n",maxi,mini);
printf("%I64d\n", Lucas(maxi-1,mini-1));
}
return 0;
}
memset(dp,0,sizeof(dp));
int n,i,j;
int l,m;
dp[1][1]=1;
for(i=2;i<=20;i++)
{
for(j=2;j<=20;j++)
{
for(l=1;l<i;l++)
{
for(m=1;m<j;m++)
{
dp[i][j]+=dp[l][m];
}
}
}
}
然后发现规律就是一个斜着的杨辉三角,杨辉三角的公式就是 第 n 行第i个数就是 (a+b)的n次方的第i个多项式的系数。
然后这里就要用到卢卡斯定理了,套了个模板就好。
#include <iostream>
#include <string.h>
#include <stdio.h>
using namespace std;
typedef __int64 LL;
LL n,m,p;
LL quick_mod(LL a, LL b)
{
LL ans = 1;
a %= p;
while(b)
{
if(b & 1)
{
ans = ans * a % p;
b--;
}
b >>= 1;
a = a * a % p;
}
return ans;
}
LL C(LL n, LL m)
{
if(m > n) return 0;
LL ans = 1;
for(int i=1; i<=m; i++)
{
LL a = (n + i - m) % p;
LL b = i % p;
ans = ans * (a * quick_mod(b, p-2) % p) % p;
}
return ans;
}
LL Lucas(LL n, LL m)
{
if(m == 0) return 1;
return C(n % p, m % p) * Lucas(n / p, m / p) % p;
}
int main()
{
int T;
p=1000000007;
while(~scanf("%I64d%I64d", &n, &m))
{
LL maxi,mini;
maxi=n+m-3;
// if(n==m)maxi+=1;
mini=m-1;
// printf("%I64d %I64d\n",maxi,mini);
printf("%I64d\n", Lucas(maxi-1,mini-1));
}
return 0;
}
相关文章推荐
- C++多任务编程简明教程 (1) - C++的多任务其实很简单
- microstation level2 cell-2
- JavaScript中的对象
- nodejs模仿优酷网站,有后端
- Javascript中的shift() 和 pop()比较
- 数值的整数次方
- hdu_3182_Hamburger Magi(状压DP)
- 实时系统动态内存算法分析dsa(一)
- Android多点触摸交互处理
- hdu_3182_Hamburger Magi(状压DP)
- lua元表
- 使用AFNetworking
- SSO单点登录
- C++:STL标准入门汇总
- 面试题24
- 由两个栈组成的队列
- SpringMVC深入理解
- 2016 百度之星 ProblemA -- (Hash+乘法逆元,利用两种方法求逆元)
- zzuoj 10403: D.山区修路
- spring 注解注入构造函数