“Earthson's Blog”的评论

Earthson 对《Lucas定理推广（组合数取模）》的评论

Earthson — Thu, 29 Dec 2011 04:04:18 +0000

@eggeek 最后那句说多了，这个其实和交换群并无关系。其实只要证明除法和乘法逆元等价，就无压力了。你的站改可能来自于对它的更直观理解，个人觉得，这种理解可能不存在。因为存在一个更匪夷所思的性质：只要结果是整数，则运算过程中即使交换了除法顺序，导致中间结果原本是不存在的情况（中间结果原本不是整数），结果也会是对的。

Earthson 对《Lucas定理推广（组合数取模）》的评论

Earthson — Thu, 29 Dec 2011 03:51:22 +0000

@eggeek 公式6其实很好理解，关键是对同余的性质深的理解。对于只包含加减乘除的同余式，首先保证原式确实是整数，并且，所有的数和模数互素。则在过程中对任意元素任意取模结果都等价（当然，模数不能是别的^^）。如果实在无法理解，请把除法看成是乘上逆元（这个过程可逆），并且，乘法情况下，这种等价性我相信你应该很清楚的（也很容易证明）。这个性质自己得好好理解才是，应该可以证明，交换群都有这个性质^^

eggeek 对《Lucas定理推广（组合数取模）》的评论

eggeek — Fri, 23 Dec 2011 01:32:39 +0000

公式6能不能再说详细点- -
按我的理解是从大往小数有k/p段被约掉了剩下k%p段应该是
(n-k +k%p)*...*(n-k+1)
___________________
(k%p)*...*1
这个与c(k%p,n%p) 同余是怎么推出来的？

Earthson 对《差分序列》的评论

Earthson — Fri, 21 Oct 2011 04:35:53 +0000

@bnuzhanyu 不好意思啊，好久没上blog了~ 你那个代码，在init的时候貌似写错了比较多的地方，包括循环参数什么的。我的代码，好多地方我用参数调的（懒得把循环参数想太清楚，所以就看运行结果，再调一下，看对了就不管了^^），所以有些循环参数看起来会比较别扭~

Earthson 对《差分序列》的评论

Earthson — Fri, 21 Oct 2011 04:31:01 +0000

@bnuzhanyu 你的代码里面，没有看到inv的初始化哎，只初始化了inv[0]~

Earthson 对《差分序列》的评论

Earthson — Fri, 21 Oct 2011 04:30:13 +0000

@bnuzhanyu 我是直接用的inv[i]是i对MOD的逆，浪费了inv[0]的^^

Earthson 对《差分序列》的评论

Earthson — Fri, 21 Oct 2011 04:25:58 +0000

@bnuzhanyu cof直接用前n+1项初始化~

bnuzhanyu 对《差分序列》的评论

bnuzhanyu — Sun, 09 Oct 2011 10:49:12 +0000

不知道怎么回事~ 那一段是inv[0]的值是1，inv[i]的值是i+1对于MOD的逆

bnuzhanyu 对《差分序列》的评论

bnuzhanyu — Sun, 09 Oct 2011 10:47:47 +0000

inv[0] = 1;
for(i = 1; i = i; j--) {
cof[j] -= cof[j - 1];
}
}

bnuzhanyu 对《差分序列》的评论

bnuzhanyu — Sun, 09 Oct 2011 10:46:40 +0000

好像懂了，不过您的calc函数貌似有点问题。
我改了一下：
int function(int k)
{
return k*k*k + k*k + 1;
}

void init()
{
int i, j;
for(int i=0; i<=order; ++i)
cof[i] = function(i);//e.g. f(i)=i^4

inv[0] = 1;
for(i = 1; i = i; j--) {
cof[j] -= cof[j - 1];
}
}
}

typec calc(typec n)
{
int i;
typec res = 0;
for(i = order; i >= 0; i--) {
res += cof[i], res %= MOD;
res *= (n - i + 1 - startfrom) % MOD * inv[i] % MOD;
res %= MOD;
}
while(res < 0) res += MOD;
return res;
}
欢迎指正~