关于(A/B)%C 以及默慈金数&&HDU 5673 Robot
2016-04-23 10:06
471 查看
首先当b与c互素
=(a*kuaisum(b,c-2)%c)
这难道就是传说中的逆元
我不知道,
我只知道
当bc互素是
=a*b^(phi(c)-1)%c
其实上面2个式子是一样的,,
然后再说下默慈金数 :在一个圆上的n个点间,画出彼此不相交的弦的全部方法的总数
公式:
大概就是 1, 2, 4, 9, 21, 51, 127, 323, 835, 2188, 5798, 15511, 41835, 113634, 310572, 853467, 2356779, 6536382, 18199284, 50852019, 142547559,
400763223, 1129760415, 3192727797, 9043402501, 25669818476, 73007772802, 208023278209, 593742784829
HDU 5673 Robot
代码如下
仅供参考
=(a*kuaisum(b,c-2)%c)
这难道就是传说中的逆元
我不知道,
我只知道
当bc互素是
=a*b^(phi(c)-1)%c
其实上面2个式子是一样的,,
然后再说下默慈金数 :在一个圆上的n个点间,画出彼此不相交的弦的全部方法的总数
公式:
大概就是 1, 2, 4, 9, 21, 51, 127, 323, 835, 2188, 5798, 15511, 41835, 113634, 310572, 853467, 2356779, 6536382, 18199284, 50852019, 142547559,
400763223, 1129760415, 3192727797, 9043402501, 25669818476, 73007772802, 208023278209, 593742784829
HDU 5673 Robot
代码如下
仅供参考
#include<stdio.h> #include<string.h> #include<algorithm> using namespace std; #include<stdlib.h> #include<math.h> #include<queue> #define MOD 1000000007 long long a[1000005]; long long suan(long long a,long long b,long long m) { long long d,t; d=1; t=a; while (b>0) { if (b%2==1) d=(d*t)%m; b/=2; t=(t*t)%m; } return d; } int main() { //printf("%I64d\n",1129760415%MOD); a[1]=1; a[2]=2; for(long long i=2; i<=1000000; i++) { a[i+1]=((2*i+3)%MOD*a[i])%MOD; a[i+1]=(a[i+1]+(((3*i)%MOD)*a[i-1]%MOD)%MOD)%MOD; a[i+1]=(a[i+1]*(suan(i+3,MOD-2,MOD))%MOD)%MOD; a[i+1]=a[i+1]%MOD; } int test; scanf("%d",&test); while(test--) { int n; scanf("%d",&n); printf("%I64d\n",a %MOD); } }
相关文章推荐
- poj-3669-Meteor Shower
- 【杭电】[1872]稳定排序
- 用MATLAB设计一个超简单的变声GUI界面
- 360前端星计划
- 简单谈谈PHP中的include、include_once、require以及require_once语句
- 一个不会OOM的毛玻璃效果
- 最小二乘法与梯度下降法
- strcpy和memcpy的差别
- c++ 读书目录
- poj-3009-Curling2.0
- 如何在审查元素里边找文件位置,求大神帮忙
- Spring加载xml配置文件的方式
- iw工具安装和使用
- 24. Swap Nodes in Pairs
- HDU 4283 You Are the One 区间dp
- 采购批打ME9F
- git使用总结(一)
- POJ_2739_Sum_of_Consecutive_Prime_Numbers_(尺取法+素数表)
- [斜率优化 DP] BZOJ 4518 [Sdoi2016]征途
- UVA-445-Marvelous Mazes