| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 1058 人关注过本帖
标题:如何简单地求两个很大的数(10^10以上)的最大公约数
只看楼主 加入收藏
李西甲
Rank: 1
等 级:等待验证会员
帖 子:16
专家分:7
注 册:2011-12-5
结帖率:100%
收藏
已结贴  问题点数:10 回复次数:11 
如何简单地求两个很大的数(10^10以上)的最大公约数
高手们,给点意见
搜索更多相关主题的帖子: 公约数 如何 
2011-12-05 13:36
czz5242199
Rank: 11Rank: 11Rank: 11Rank: 11
等 级:小飞侠
威 望:4
帖 子:660
专家分:2400
注 册:2011-10-26
收藏
得分:3 
还是辗转相除,如果数据超过了long long,就自己写高精度运算
2011-12-05 14:47
原味好
Rank: 4
来 自:西安
等 级:业余侠客
帖 子:59
专家分:250
注 册:2011-11-29
收藏
得分:3 
干什么都得现实点
2011-12-05 14:52
李西甲
Rank: 1
等 级:等待验证会员
帖 子:16
专家分:7
注 册:2011-12-5
收藏
得分:0 
回复 2楼 czz5242199
那样太慢了,呵呵
2011-12-05 22:54
李西甲
Rank: 1
等 级:等待验证会员
帖 子:16
专家分:7
注 册:2011-12-5
收藏
得分:0 
回复 4楼 李西甲
试试看而已
2011-12-05 22:54
laoyang103
Rank: 19Rank: 19Rank: 19Rank: 19Rank: 19Rank: 19
来 自:内蒙古包头
等 级:贵宾
威 望:19
帖 子:3082
专家分:11056
注 册:2010-5-22
收藏
得分:3 
超级最小公倍数
时间限制:8000 ms  |  内存限制:16384 KB
描述
给2个正整数a,b(1<=a,b<=10^100),求a和b的最小公倍数。

输入
输入包含多组数据,每组数据一行,包含两个正整数a和b,中间以一个空格隔开。输入以0 0结束。

输出
每组数据输出一行,为a,b的最小公倍数。

样例输入
123 321
123456789 987654321
0 0
样例输出
13161
13548070123626141
http://www.

[ 本帖最后由 laoyang103 于 2011-12-5 23:20 编辑 ]

                                         
===========深入<----------------->浅出============
2011-12-05 23:17
laoyang103
Rank: 19Rank: 19Rank: 19Rank: 19Rank: 19Rank: 19
来 自:内蒙古包头
等 级:贵宾
威 望:19
帖 子:3082
专家分:11056
注 册:2010-5-22
收藏
得分:0 
程序代码:
#include<stdio.h>
#include<string.h>
void output(char * a, int an)
{
    int i;
    for(i = an - 1; i >= 0; i--) putchar(a[i] + '0');
    putchar('\n');
}
void format(char * a, int an)
{
    int i, j;
    char t;
    for(i = 0, j = an - 1; i < j; i++, j--)
    {
        a[i] -= '0';
        a[j] -= '0';
        t = a[i];
        a[i] = a[j];
        a[j] = t;
    }
    if(i == j) a[i] -= '0';
}
void mul(char * a, int * an, char * b, int bn)
{
    char t[256], s[256] = {0};
    int i, j, k, f;
    for(i = 0; i < bn; i++)
    {
        if(b[i] == 0) continue;
        for(j = 0; j < i; t[j++] = 0);
        for(f = 0, k = 0; k < *an; k++)
        {
            t[j] = b[i] * a[k] + f;
            f = t[j] / 10;
            t[j++] %= 10;
        }
        if(f) t[j++] = f;
        for(f = 0, k = 0; k < j; k++)
        {
            s[k] += t[k] + f;
            if(s[k] >= 10){ s[k] -= 10; f = 1;} else f = 0;
        }
        if(f)s[j++] = 1;
    }
    for(i = 0; i < j; a[i] = s[i++]);
    *an = j;
}
void div(char * a, int * an, char * b, int bn)
{
    char t[256], r[256] = {0};
    int f, i, j;
    if(*an < bn)
    {
        a[0] = 0;
        a[1] = 0;
        *an = 1;
        return;
    }
    for(i = *an - bn; i >= 0; i--)
    {
        for(f = 0, j = 0; j < bn; j++)
        {
            t[j] = a[i + j] - b[j] - f;
            if(t[j] < 0){ t[j] += 10; f = 1;} else f = 0;
        }
        if(a[bn + i] - f >= 0)
        {
            a[bn + i] -= f;
            for(j = 0; j < bn; a[i + j] = t[j++]);
            r[i]++;
            i++;
        }
    }
    for(i = *an - bn; i && !r[i]; i--);
    for(j = i; j >= 0; a[j] = r[j--]);
    *an = i + 1;
    a[*an] = 0;
}

void mod(char * a, int * an, char * b, int bn)
{
    if(*an<bn)
        return ;
    char c[101] = {0};
    int i,j,k,low;
    for(i = *an-bn;i>=0;--i)
    {
        low = 0;
        for(j = 0;j<bn;++j)
        {
            c[j] = a[i+j] - b[j] - low;
            if(c[j]<0)
            {c[j] += 10;low = 1;}
            else
                low = 0;
        }
        if(a[i+bn] - low >= 0)
        {
            a[i+bn] -= low;
            for(k = 0;k<bn;++k)
                a[i+k] = c[k];
            i++;
        }
    }
    for(i = bn-1;!a[i]&&i>=0;i--);
    *an = i+1;
}

void lcm(char *a, int * an, char *b, int bn)
{
    char ta[256], tb[256], *pa, *pb, *pt;
    int pan, pbn, ptn, i;
    for(i = 0; i <= *an; ta[i] = a[i++]);
    for(i = 0; i <= bn; tb[i] = b[i++]);
    pa = ta;
    pb = tb;
    pan = *an;
    pbn = bn;
    mod(pa, &pan, pb, pbn);
    while(pan > 1 || pa[0])
    {
        pt = pa;
        pa = pb;
        pb = pt;
        ptn = pan;
        pan = pbn;
        pbn = ptn;
        mod(pa, &pan, pb, pbn);
    }
    div(a, an, pb, pbn);
    mul(a, an, b, bn);
}
int main()
{
    char a[256], b[256];
    int an, bn;
    while(scanf("%s %s", a, b), a[0] > '0')
    {
        an = strlen(a);
        bn = strlen(b);
        format(a, an);
        format(b, bn);
        //mod(a, &an, b, bn);
        lcm(a, &an, b, bn);
        output(a, an);
    }
    return 0;
}
AC代码  留着当大数模板用吧

                                         
===========深入<----------------->浅出============
2011-12-05 23:19
LOVE_IT123
Rank: 1
等 级:新手上路
帖 子:2
专家分:2
注 册:2011-12-5
收藏
得分:3 
试试看
2011-12-05 23:28
czz5242199
Rank: 11Rank: 11Rank: 11Rank: 11
等 级:小飞侠
威 望:4
帖 子:660
专家分:2400
注 册:2011-10-26
收藏
得分:0 
回复 4楼 李西甲
你确定太慢了?你能分析出他的时间复杂度吗?还是你就是觉得太慢了?
2011-12-06 09:06
李西甲
Rank: 1
等 级:等待验证会员
帖 子:16
专家分:7
注 册:2011-12-5
收藏
得分:0 
回复 9楼 czz5242199
我是个新手,,感觉而已,莫见怪!
2011-12-06 22:13
快速回复:如何简单地求两个很大的数(10^10以上)的最大公约数
数据加载中...
 
   



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

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