May the force be with you!
littlekid

littlekid@R2
posts - 49,  comments - 17,  trackbacks - 0

问题来源:
HUST07校赛

原题见:http://acm.whu.edu.cn/oak/problem/problem.jsp?problem_id=1338

提交方式:WOJ1338

问题描述:
    对于N(1<=N<=80)个数A[1...N](1<=A[i]<=500),他们的和为sum,求sum!/(A[1]!*A[2]!*...*A[N]!)%40009。

解题过程:
    对于这个题目,我当时就推出了上面的公式,然后就卡了,不知道后面怎么办——这些数可是非常大。

    其实这个问题的重点就在于运用扩展欧几里德(感谢Felicia的指导):

对于 a/b%m = ans, 求 ans。

a = a%m, b = b%m
ans = (a % m)*(x % m) % m  (x为b的逆元)

求逆元则利用扩展欧几里德:
对于 b*x = 1(mod m)
可以求b*x + m*y = 1的解( 用extennd_Euclid(b, m, x, y) )!
然后把 x 映射到 [0,m)区间,带入上式, 即得解。


附代码:

 1 
 2 int extend_Euclid(int a, int b, int &x, int &y)
 3 {
 4     if (b == 0)
 5     {
 6         x = 1;
 7         y = 0;
 8          return a;
 9     }
10     else
11     {
12         int ans = extend_Euclid(b, a % b, x, y); /////
13         int t = x;
14         x = y;
15         y = t - (a / b) * y;
16         return ans;
17     }
18 }
19 
20 





posted on 2008-03-14 16:46 littlekid@R2 阅读(1235) 评论(4)  编辑 收藏 引用 所属分类: Problem SolvingPure TheoryMemo

FeedBack:
# re: 【数论】扩展欧几里德的一个妙用
2008-03-17 21:23 | 长江三峡
学习  回复  更多评论
  
# re: 【数论】扩展欧几里德的一个妙用
2008-03-17 21:42 | sdf
感觉ACM玩的东西都好复杂啊!!  回复  更多评论
  
# re: 【数论】扩展欧几里德的一个妙用
2008-04-07 21:48 | IceDragon
楼主你好,我也是用扩展欧几里得做的,但是我一直WA
可否借你的代码看下--不胜感激
这题目还有什么陷阱么??????
我邮箱--wangzhaoren@gmail.com
  回复  更多评论
  
# re: 【数论】扩展欧几里德的一个妙用
2008-04-09 12:53 | littlekid
@IceDragon
估计OJ出现不明错误,或者后来加强数据了(最近交的都没过,我贴了别人AC的代码都没过=,=!)

不过你可以考虑一下越界问题——用long long试试  回复  更多评论
  

标题  
姓名  
主页
验证码 *
内容(提交失败后,可以通过“恢复上次提交”恢复刚刚提交的内容)  
  登录  使用高级评论  新用户注册  返回页首  恢复上次提交      
[使用Ctrl+Enter键可以直接提交]




你是第 free hit counter 位访客


联系博主


<2008年5月>
27282930123
45678910
11121314151617
18192021222324
25262728293031
1234567

常用链接

留言簿

随笔分类(52)

随笔档案(49)

文章档案(1)

ACM/ICPC

技术综合

最新随笔

搜索

  •  

积分与排名

  • 积分 - 12180
  • 排名 - 126

最新评论

阅读排行榜

评论排行榜

60天内阅读排行