题目链接:http://acm.swust.edu.cn/problem/0247/

Time limit(ms): 1000        Memory limit(kb): 65535
 
Description
在很久很久以前,有个臭美国王。一天,他得到了一件新衣,于是决定巡城看看老百姓的反应(囧)。于是他命令可怜的宰相数一下他有多少种巡城方案。 
城市是一个N*M的矩形,(N+1)*(M+1)条街把城市分成了N*M块。国王从左下角出发,每次只能向右或向上走,右上角是终点。 
请你帮帮可怜的宰相。

 
Input
每组测试数据有3个值 n m p(0< n <=1000000000,0 < m <=1000000000,p为质数且p<400). 
 
Output
多组测试数据 0 0 0结束。
输出方案数 因为数据太大 请 mod p

 
Sample Input
2 3 97
20 40 37
0 0 0
Sample Output
10
32
 
Hint
2009暑期ACM集训组队赛《一》 -- 剑问苍生
 
解题思路:从左下角到右上角的方案数,那么(假设方阵大小n*m)
     (1)向上必须移动n格,向右必然移动m格
     (2)只需要以向上或向右移动为基准求C(n+m) n 或者C(n+m) m的组合数即可,好好想想为什么吧~~~~
        可以有一个优化组合数上面的去n,m中的较小值(组合数性质)
     (3)发现n,m<=1000000000那么在求组合数的时候肯定会直接爆掉,那怕是在求的过程中不断MOD 数p,
       那么我们以p的上限值为界求组合数,运用Lucas定理求解
 
Lucas定理:(详情可以戳戳这里:http://www.cnblogs.com/zyxStar/p/4588362.html)

      Lucas定理用来求C(a,b)mod p的值,其中p为素数。
      数学表达式为:
      Lucas(a,b,q)=C(a%q,b%q)*Lucas(a/p,b/p,p);
      Lucas(a,0,q)=0;

      模运算与基本四则运算有些相似,但是除法例外。其规则如下:
      (a + b) % p = (a % p + b % p) % p (1)
      (a - b) % p = (a % p - b % p) % p (2)
      (a * b) % p = (a % p * b % p) % p (3)
      a ^ b % p = ((a % p)^b) % p (4)

代码如下:

 #include <iostream>
#include <cstring>
#include <algorithm>
#define maxn 405
using namespace std;
long long dp[maxn][maxn], n, m, p;
long long Lucas(long long n, long long m, long long p){
//Lucas定理
if (n <= p && m <= p)
return dp[n][m];
else
return (Lucas(n / p, m / p, p)*dp[n%p][m%p]) % p;
} //杨辉三角求组合数
void init(){
int i, left, right;
for (i = ; i <= p; i++){
dp[i][] = dp[i][i] = ;
left = , right = i - ;
while (left <= right){
dp[i][left] = (dp[i - ][left - ] + dp[i - ][left]) % p;
dp[i][right--] = dp[i][left++];//组合数性质dp[i][j]=dp[i][i-j];
}
}
}
int main(){
while (cin >> n >> m >> p && n && m && p){
memset(dp, , sizeof(dp));
init();
cout << Lucas(n + m, min(n, m), p) << endl;
}
return ;
}

乱搞的Java代码如下:

 import java.math.*;
import java.io.*;
import java.util.*;
public class Main{
static int DP[][] = new int[401][401];
static int p;
public static void main(String[] args)
{
Scanner cin = new Scanner(System.in);
int a, b;
while (cin.hasNext()){
a = cin.nextInt();
b = cin.nextInt();
p = cin.nextInt();
if (a == 0)
break;
initDp();
System.out.println(sloveRe(a + b, b > a ? a : b));
}
}
private static void initDp(){
int i, j;
for (i = 0; i <= p; i++){
DP[i][0] = 1;
}
for (i = 1; i <= p; i++){
for (j = 1; j <= p; j++){
DP[i][j] = (DP[i - 1][j] + DP[i - 1][j - 1]) % p;
}
}
}
private static int sloveRe(int n, int m){
if (n <= p&&m <= p)
return DP[n][m];
else
return (sloveRe(n / p, m / p)*DP[n%p][m%p]) % p;
}
}

[Swust OJ 247]--皇帝的新衣(组合数+Lucas定理)的更多相关文章

  1. uoj86 mx的组合数 (lucas定理+数位dp+原根与指标+NTT)

    uoj86 mx的组合数 (lucas定理+数位dp+原根与指标+NTT) uoj 题目描述自己看去吧( 题解时间 首先看到 $ p $ 这么小还是质数,第一时间想到 $ lucas $ 定理. 注意 ...

  2. 【BZOJ-4591】超能粒子炮·改 数论 + 组合数 + Lucas定理

    4591: [Shoi2015]超能粒子炮·改 Time Limit: 10 Sec  Memory Limit: 256 MBSubmit: 95  Solved: 33[Submit][Statu ...

  3. luogu4345 [SHOI2015]超能粒子炮·改(组合数/Lucas定理)

    link 输入\(n,k\),求\(\sum_{i=0}^k{n\choose i}\)对2333取模,10万组询问,n,k<=1e18 注意到一个2333这个数字很小并且还是质数这一良好性质, ...

  4. 【(好题)组合数+Lucas定理+公式递推(lowbit+滚动数组)+打表找规律】2017多校训练七 HDU 6129 Just do it

    http://acm.hdu.edu.cn/showproblem.php?pid=6129 [题意] 对于一个长度为n的序列a,我们可以计算b[i]=a1^a2^......^ai,这样得到序列b ...

  5. 组合数(Lucas定理) + 快速幂 --- HDU 5226 Tom and matrix

    Tom and matrix Problem's Link:   http://acm.hdu.edu.cn/showproblem.php?pid=5226 Mean: 题意很简单,略. analy ...

  6. 【组合数+Lucas定理模板】HDU 3037 Saving

    acm.hdu.edu.cn/showproblem.php?pid=3037 [题意] m个松果,n棵树 求把最多m个松果分配到最多n棵树的方案数 方案数有可能很大,模素数p 1 <= n, ...

  7. CodeForces-451E:Devu and Flowers (母函数+组合数+Lucas定理)

    Devu wants to decorate his garden with flowers. He has purchased n boxes, where the i-th box contain ...

  8. HDU3037Saving Beans(组合数+lucas定理)

    Problem Description Although winter is far away, squirrels have to work day and night to save beans. ...

  9. 组合数取模Lucas定理及快速幂取模

    组合数取模就是求的值,根据,和的取值范围不同,采取的方法也不一样. 下面,我们来看常见的两种取值情况(m.n在64位整数型范围内) (1)  , 此时较简单,在O(n2)可承受的情况下组合数的计算可以 ...

随机推荐

  1. BZOJ 1189: [HNOI2007]紧急疏散evacuate( BFS + 二分答案 + 匈牙利 )

    我们可以BFS出每个出口到每个人的最短距离, 然后二分答案, 假设当前答案为m, 把一个出口拆成m个表示m个时间, 点u到出口v的距离为d, 那么u->v的[d, m]所有点连边, 然后跑匈牙利 ...

  2. 用Dart开发Android应用(sky应用)

    此文运用的是优雅的Markdown而书 流弊的Google公司,为了解决Android的卡顿问题,为此下了是不少的决心,这不最新研发的Sky移动应用UI框架,在用户界面和高质量的UI交互设计,很是让人 ...

  3. PHP的环境搭建

    下载开发环境 wampserver 下载sublime text 2 sublime使用技巧 1:安装漂亮的编程字体http://pan.baidu.com/s/1xMex9 下载"程序编写 ...

  4. 正则语法笔记-regular expression note

    参考文档:python正则表达式 正则表达式定义:正则是一门高度专业编程语言,内嵌在其他语言(python re模块)中使用.正则表达式包含元字符(metacharacter)列表,列表如下: . ^ ...

  5. web应用中Spring ApplicationContext的动态更新

    在web应用中时常需要修改配置,并动态的重新加载ApplicationContext.比如,设置和切换数据库.以下给出一个方法,并通过代码验证可行性. 方法的基本思路是,为WebApplication ...

  6. Azure 网站的新增功能:可配置的环境变量

     编辑人员注释:本文章由 WindowsAzure 网站团队的项目经理Erez Benari撰写. Azure最常用的强大功能之一是 XML文档转换 (XDT),通过此功能,您可以在Windows ...

  7. HDU 1722 Cake

    #include<cstdio> int gcd(int m, int n) { ?n:gcd(n % m, m); } int main() { int m, n; while(scan ...

  8. Network of Schools(强连通分量+缩点) (问添加几个点最少点是所有点连接+添加最少边使图强连通)

    Network of Schools Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 13801   Accepted: 55 ...

  9. android经常使用的电话操作

    给大家分享一下我的一个Android工具类,能够获取手机里面的各种信息,包含拨打电话. 获取全部联系人姓名及电话,插入联系人姓名及电话,插入联系人姓名及电话.插入通话记录.获取用户全部短信.批量插入短 ...

  10. HA for openstack

    mysql ha instance ha openstack博客:http://blog.csdn.net/tantexian/article/list/2 使用eclipse远程调试openstac ...