编程珠玑2中的习题

/*

 * Filename: 

 * Auther: mlj

 *第一章 第四题  产生0-n-1 之间k个随机顺序的随机整数

 * Create date: 2013/ 3/12

 */

 #include <stdlib.h> 

#include <stdio.h> 

#include<stdio.h>

#include <time.h>  

#define N 10000000

#define K 1000000



void swap(int *pm,int *pn)      /*必须用指针进行交换*/

{

   int temp;

   temp=*pm;

   *pm=*pn;

   *pn=temp;

}



int randint(l,u)

{

    return rand()%(u-l+1)+l;

}

int a[N];

int main(void)

{

    int   i;

    FILE *fp =NULL;

    fp = fopen("test.dat","wb"); 

    srand( (unsigned)time( NULL ) ); /*定义这个可以产生不同的随机数*/

    for(i=0;   i<N;   i++) {a[i]=i;}

    for(i=0;   i<K;   i++)

    {

     /* pa=&a[i]; pb=&a[rand()%i+1];*/

      swap(&a[i], &a[rand()%(N-i)+i]);     

      fwrite(&a[i],sizeof(int),1,fp);

      //printf("%d ",a[i] );

    }

    fclose(fp);     

}

int a[N]; 不能定义成局部变量，局部变量在栈区，大小有限。

/*

 * Filename: 

 * Auther: mlj

 *第一章 第3题  产生的k个 随机整数 排序

 * Create date: 2013/ 3/12

 */

 #include <stdlib.h> 

#include <stdio.h> 

#include<stdio.h>

#include <time.h>  

#define N 10000000

#define K 1000000



#define BITSPERWORD 32

#define SHIFT 5

#define MASK 0x1f

int a[1+N/BITSPERWORD];



void set(int i) {        a[i>>SHIFT] |=  (1<<(i & MASK)); }

void clr(int i) {        a[i>>SHIFT] &= ~(1<<(i & MASK)); }

int test(int i) {return  a[i>>SHIFT] &   (1<<(i & MASK)); }



int intcomp(const void *a,const void *b)

{

    return *((int*)a) - *((int*)b);

}



int main(void)

{

    int   i;

    FILE *fp =NULL;

    float starttime = 0.0f;

    int *stream =(int*) malloc( sizeof(int)*K);

    int *p=NULL;

    fp = fopen("test.dat","rb"); 

    

    fread(stream,sizeof(int),K,fp) ;

    //printf("读文件耗时 %0.2f",clock()*0.001f-starttime);

    starttime = clock()*0.001f;

    for(i=0;   i<N;   i++) 

        clr(i);  

    for(i=0;   i<K;   i++)      

         set(*(stream+i)); 

    p = stream;

    for(i=0;   i<N;   i++)

    {

     /* pa=&a[i]; pb=&a[rand()%i+1];*/

      if(test(i))     

         *(p++) = i;

      //printf("%d ",a[i] );

    }

    printf("总耗时 %0.2f",clock()*0.001f-starttime);



    fseek(fp,0,SEEK_SET);

    fread(stream,sizeof(int),K,fp) ;

    starttime = clock()*0.001f;

    qsort(stream,K,sizeof(int),intcomp);

    printf("总耗时 %0.2f",clock()*0.001f-starttime);

    fclose(fp);     

    free(stream);

    return 0;

}

总耗时 0.10总耗时 0.17

C/位图0.1s  ,qsort 也才0.17s .