鱼C论坛

 找回密码
 立即注册
查看: 2522|回复: 0

一道作业题,有关低级处理,效率,Karatsuba算法

[复制链接]
发表于 2014-7-4 17:37:18 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能^_^

您需要 登录 才可以下载或查看,没有账号?立即注册

x
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#define SWAP(a,b) { a ^= b; b ^= a; a ^= b; }

unsigned long int karatsuba(unsigned long int a, unsigned long int b);
int main () {
    unsigned long int a, b;
    srand(time(NULL));
    
    
#ifndef TEST
    a=rand(); b=rand();
    printf("%ld*%ld=%ld\n",a,b,karatsuba(a,b));
#endif
    
#ifdef TEST
    for(int i=0; i< 1000000; i++) {
        a=rand(); b=rand();
        if(karatsuba(a,b)!=a*b) {
            fprintf(stderr,"Error: a=%ld, b=%ld\n",a,b);
            exit(-1);
        }
    }
    
#endif
}

unsigned long int karatsuba(unsigned long int a, unsigned long int b) {
    int i, n, N;
    unsigned long int x0, y0;
    unsigned long int z0, z1, z2;
    if(a<b) SWAP(a,b);
    if(b==0) return 0;
    for(n=-1, i = 1; i <= b; i<<=1, n++);  /* not optimal */
    for(N=n; i <= a; i<<=1, N++);
    y0=b&((1<<n)-1);
    x0=a&((1<<N)-1);
    z2=1<<(N-n);
    z1=karatsuba(z2,y0)+x0;
    z0=karatsuba(x0,y0);
    return (z2<<(2*n))+(z1<<n)+z0;
}
                                                                                                                                Explainwhat specific adjustments were made to the algorithm inorder to improve the efficiency.



                                                                                                Write an efficient function to replace the for loops marked as “notoptimal”. Hint: you can use some assembler instructions (keyword asm), some build-incompiler functions (dependent on your compiler), or the operators &, |, << and >>
                               
                       
               
                               
                       
               


题目大致意思是说这段代码用karatsuba算法在算两个数乘积,但是用了low level的处理,让解释一下它的高效在哪.
谁能解释一下那个karatsuba函数里的代码是在干些什么啊?顺便,如何改进/* not optimal */ 的地方?
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

小黑屋|手机版|Archiver|鱼C工作室 ( 粤ICP备18085999号-1 | 粤公网安备 44051102000585号)

GMT+8, 2024-11-24 14:50

Powered by Discuz! X3.4

© 2001-2023 Discuz! Team.

快速回复 返回顶部 返回列表