| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 1037 人关注过本帖
标题:快速幂,哪错了?
只看楼主 加入收藏
jianyuling00
Rank: 3Rank: 3
来 自:琼州学院
等 级:论坛游侠
帖 子:89
专家分:177
注 册:2012-3-17
收藏
得分:0 
LS的运行结果是一样的吧
2012-03-19 01:05
czz5242199
Rank: 11Rank: 11Rank: 11Rank: 11
等 级:小飞侠
威 望:4
帖 子:660
专家分:2400
注 册:2011-10-26
收藏
得分:0 
一开始没看懂你写的,你的思路一开始绕了我一下,这么改改

程序代码:
#include<cstdio>
__int64 quickpow(__int64 a,__int64 k,__int64 mod)
{
    __int64 temp=1;
    while(k>0){
        if(k&1)
            temp=(a%mod)*(temp%mod)%mod;  
        a=(a%mod)*(a%mod)%mod;      
        k/=2;
    }
    return temp;
}
int main()
{
    __int64 a,k,m;
    while(~scanf("%I64d%I64d%I64d",&a,&k,&m)){
        printf("%I64d\n",quickpow(a,k,m));
    }
    return 0;
}
2012-03-19 01:20
cb_1212
Rank: 1
等 级:新手上路
帖 子:126
专家分:5
注 册:2011-4-28
收藏
得分:0 
回复 12楼 czz5242199
嗯嗯,计算过程中溢出了吧。谢谢了!
2012-03-19 12:34
cb_1212
Rank: 1
等 级:新手上路
帖 子:126
专家分:5
注 册:2011-4-28
收藏
得分:0 
照czz改了下,怎么一直错呢,奇怪。
http://218.194.91.48/acmhome/problemdetail.do?&method=showdetail&id=1323
2012-03-19 23:05
cb_1212
Rank: 1
等 级:新手上路
帖 子:126
专家分:5
注 册:2011-4-28
收藏
得分:0 
a=(a%mod)*(a%mod)%mod;//mod也可能取到10^18,a<mod时(a%mod)*(a%mod)爆__int64了,算法还要优化啊。
2012-03-20 00:15
beyondyf
Rank: 19Rank: 19Rank: 19Rank: 19Rank: 19Rank: 19
等 级:贵宾
威 望:103
帖 子:3282
专家分:12654
注 册:2008-1-21
收藏
得分:6 
再拆分一个乘数就可以了

重剑无锋,大巧不工
2012-03-20 11:50
laoyang103
Rank: 19Rank: 19Rank: 19Rank: 19Rank: 19Rank: 19
来 自:内蒙古包头
等 级:贵宾
威 望:19
帖 子:3082
专家分:11056
注 册:2010-5-22
收藏
得分:0 
回复 16楼 beyondyf
这个a貌似太大了。。。。二分要溢出

求杨大哥的代码

                                         
===========深入<----------------->浅出============
2012-03-20 15:47
beyondyf
Rank: 19Rank: 19Rank: 19Rank: 19Rank: 19Rank: 19
等 级:贵宾
威 望:103
帖 子:3282
专家分:12654
注 册:2008-1-21
收藏
得分:0 
离开一下,一会儿写给你

重剑无锋,大巧不工
2012-03-20 15:56
beyondyf
Rank: 19Rank: 19Rank: 19Rank: 19Rank: 19Rank: 19
等 级:贵宾
威 望:103
帖 子:3282
专家分:12654
注 册:2008-1-21
收藏
得分:0 
久等了,这段代码虽然AC了,但时间很不理想,近2s,很想知道是哪组数据如此耗时。有时间我再重新构建一个算法。

程序代码:
#include<stdio.h>
long long AxBmodM(long long a, long long b, long long m)
{
    long long r;
    for(r = 0; b; b >>= 3, a = (a << 3) % m)
        if(b & 7) r = (r + a * (b & 7)) % m;
    return r;
}

long long quickPow(long long a, long long k, long long m)
{
    long long r;
    if(a == 0) return 0;
    if(k == 0) return 1;
    if(k == 1) return a % m;
    r = quickPow(a, k >> 1, m);
    r = AxBmodM(r, r, m);
    if(k & 1) r = AxBmodM(r, a, m);
    return r;
}

int main()
{
    long long a, k, m;
    while(scanf("%I64d%I64d%I64d", &a, &k, &m) != EOF)
        printf("%I64d\n", quickPow(a, k, m));
    return 0;
}

重剑无锋,大巧不工
2012-03-20 17:41
cb_1212
Rank: 1
等 级:新手上路
帖 子:126
专家分:5
注 册:2011-4-28
收藏
得分:0 
拆分乘数算的。
程序代码:
#include<cstdio>
__int64 quickpow(__int64 a,__int64 k,__int64 mod)
{
    __int64 aa,sum,temp=1;
    while(k){
        if(k&1)
            temp=a*temp%mod;
        if(a<1000000000) a=a*a%mod;
        else
        {
            aa=a;
            sum=0;
            while(aa){
                sum=(sum+a*(aa%10)%mod)%mod;
                aa/=10;
                a=a*2%mod;
                a=a*5%mod;
            }
            a=sum;            
        }
        k>>=1;
    }
    return temp;
}
int main()
{
    __int64 a,k,m;
    while(~scanf("%I64d%I64d%I64d",&a,&k,&m)){
        printf("%I64d\n",quickpow(a%m,k,m));
    }
    return 0;
}
2012-03-20 20:13
快速回复:快速幂,哪错了?
数据加载中...
 
   



关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.013901 second(s), 7 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved