幂模算法

恶心的狐狸

呵呵,SB还知道这个成语?
行,模仿我,写一个生成很大质数的算法,这个对敢说这个话的人那肯定是小菜一碟了吧.

moranhuishou

下面引用由恶心的狐狸在 2009/06/17 11:45am 发表的内容：
呵呵,SB还知道这个成语?
行,模仿我,写一个生成很大质数的算法,这个对敢说这个话的人那肯定是小菜一碟了吧.

骂人可以，要拿出点真货才有资格！干嚎算个屁呀！！！
可以告诉你，“写一个生成很大质数的算法”是小菜一碟！！！！！！！！！！

恶心的狐狸

呵呵,编过程序吗?
知道什么叫算法吗?(我不需要你回答数学里的定义,因为这个对一般人来说太难了)

恶心的狐狸

我要告诉你的是,蠢材,
我一楼的那个东西叫算法,可以解决所有的
a^b%c
这叫解决一类问题,而不是解决单个问题,人蠢不能怨政府

moranhuishou

下面引用由恶心的狐狸在 2009/06/17 11:52am 发表的内容：
呵呵,编过程序吗?
知道什么叫算法吗?(我不需要你回答数学里的定义,因为这个对一般人来说太难了)

我都给你计算那么大的数字了，你还问这些？你真傻呀！

moranhuishou

不要再班门弄斧丢人了。我可以不客气地告诉你：
利用程序大数a^b%c运算是我的强项，你大概也就会纸上谈兵，实际排兵布阵还差点！！！

恶心的狐狸

呵呵,还真不服气?
判定一下,
9153837170793543235597585083490030479955331292055521168283445592938226952061614515816656193273439207474780594020068103215636057295338785038298689948720343
和
9153837170793543235597585083490030479955331292055521168283445592938226952061614515816656193273439207474780594020068103215636057295338785038298689948720341
和
9153837170793543235597585083490030479955331292055521168283445592938226952061614515816656193273439207474780594020068103215636057295338785038298689948720339
这三个数中谁满足2^(n-1)%2==1
的fermat小定理判定条件?

fleurly

告诉某个不懂装懂的家伙
1.
n=45646487979878974564446465465464646461313415464647989789652379时候
9^n mod n = 45646487979878974564446465465464646461313415464647989789652379
2.
n=29341213421321329341213421321329341213421321352525656547752525656547729341213421
的时候
911 ^n mod n = 27227059835155552772739572788935085975550644709642717731307390139894005227663848
3.
c=321315464567857987946584543541352132134165465458978979879879798797987987123546235413
a=113
b=c-1
a^b mod c = 1
附件是按照 狐狸 的算法写的程序， Python的。有懂计算机的朋友可以看一下究竟是谁在骗人。

fleurly

附件没法上传，就贴在这里了
#!/usr/bin/python
#a, b , c
# a^b mod c
import math
def m_mod(a,b,c) :
        l_v = math.log(c, a)
        if(l_v > b) :
                &#35;print "pow(%d, %d) < %d"%(a, b, c )
                return pow(a,b)
        i_l_v = int(l_v) + 1
        count=b / i_l_v
        mod_p = (pow(a,i_l_v) )% c
        &#35;print "pow(%d, %d) mod %d is %d" % (a, i_l_v, c, mod_p)
        &#35;print "Should pow %d times ."% (count)
        t_mod1 = 1
        _l_v = int (math.log(c, a)) +1
        if(count > _l_v * 2):
                t_mod1 = m_mod(mod_p, count, c)
        else:
                t_mod1 = pow(mod_p, count ) % c
        r= b % i_l_v
        t_mod2 = pow(a, (b % i_l_v))
        res = (t_mod1 * t_mod2 )% c
        &#35;print "%d ^ %d mod %d = %d " %(a,b,c,res)
        return res

a1=911
b1=45646487979878974564446465465464646461313415464647989789652379
c1=45646487979878974564446465465464646461313415464647989789652379
res1 = m_mod(a1, b1, c1)
a2=9
b2=29341213421321329341213421321329341213421321352525656547752525656547729341213421
c2=29341213421321329341213421321329341213421321352525656547752525656547729341213421
res2 = m_mod(a2, b2, c2)
c3=321315464567857987946584543541352132134165465458978979879879798797987987123546235413L
a3=113
b3=c3-1
res3 = m_mod(a3, b3, c3)
print "%d ^ %d mod %d = %d " %(a1,b1,c1,res1)
print "%d ^ %d mod %d = %d " %(a2,b2,c2,res2)
print "%d ^ %d mod %d = %d " %(a3,b3,c3,res3)

fleurly

下面引用由恶心的狐狸在 2009/06/17 00:11pm 发表的内容：
呵呵,还真不服气?
判定一下,
9153837170793543235597585083490030479955331292055521168283445592938226952061614515816656193273439207474780594020068103215636057295338785038298689948720343
和
...

狐狸别担心，我是软件工程师，专门写程序的。我给你写程序了。