数学研发网设为首页收藏本站
切换到窄版

 找回密码
 欢迎注册
数学研发论坛»论坛 › 【算法研发】 › 编程擂台 › 用过程语言实现group by cube的代码
返回列表 发新帖
查看: 11061|回复: 9

[擂台] 用过程语言实现group by cube的代码

[复制链接]
〇〇
发表于 2010-2-25 10:10:12 | 显示全部楼层 |阅读模式

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要 登录 才可以下载或查看，没有账号？欢迎注册

×
有100万行数据，要求对c1,c2,c3,c4作各种组合的分组汇总，计算行数和v1的合计,即
  2. 0( c1, c2, c3, c4);
  3. 1( c1, c2, c3,NULL);
  4. 2( c1, c2,NULL, c4);
  5. 3( c1, c2,NULL,NULL);
  6. 4( c1,NULL, c3, c4);
  7. 5( c1,NULL, c3,NULL);
  8. 6( c1,NULL,NULL, c4);
  9. 7( c1,NULL,NULL,NULL);
  10. 8(NULL, c2, c3, c4);
  11. 9(NULL, c2, c3,NULL);
  12. 10(NULL, c2,NULL, c4);
  13. 11(NULL, c2,NULL,NULL);
  14. 12(NULL,NULL, c3, c4);
  15. 13(NULL,NULL, c3,NULL);
  16. 14(NULL,NULL,NULL, c4);
  17. 15(NULL,NULL,NULL,NULL);
  18. /*源数据*/
  19. int main()
  20. {
  21. char* c1=(char*)malloc(1E6L);
  22. char* c2=(char*)malloc(1E6L);
  23. char* c3=(char*)malloc(1E6L);
  24. char* c4=(char*)malloc(1E6L);
  25. int* v1=(int* )malloc(4E6L);
  26. __int64* v2;//store sum(v1);
  27. long t=clock();
  28. for (int i=0;i<1E6L;i++)
  29. {
  30. int j=i+1;
  31. c1[i]=j%4;
  32. c2[i]=j%5;
  33. c3[i]=j%67;
  34. c4[i]=(j%13?j%113:127); //用127模拟NULL
  35. v1[i]=j;
  36. }
复制代码
用尽量少的空间和时间
毋因群疑而阻独见　　毋任己意而废人言
毋私小惠而伤大体　　毋借公论以快私情
回复

使用道具 举报

qianyb
发表于 2010-2-25 16:44:11 | 显示全部楼层
这种计算空间占用不大，时间很难再优化了（先把相同的行数统计出来，再汇总）
毋因群疑而阻独见　　毋任己意而废人言
毋私小惠而伤大体　　毋借公论以快私情
回复 支持 反对

使用道具 举报

〇〇
 楼主| 发表于 2010-2-25 20:38:16 | 显示全部楼层
很多地方没有优化，应该还能提高100毫秒
  2. lenovo t61
  3. 453ms,recn=152760
  4. 547ms,recn=234600,c1=-1,c2=127,c3=127,c4=127,count=1000000,sum=500000500000
  6. hp dl360
  7. 483ms,recn=152760
  8. 577ms,recn=234600,c1=-1,c2=127,c3=127,c4=127,count=1000000,sum=500000500000
  10. /* compiled by cl of vc2003*/
  11. #include <time.h>
  12. #include <stdio.h>
  13. #include <stdlib.h>
  14. #include <map>
  15. #include <iostream>
  16. #define add_key {if((key_recn.find(key))== key_recn.end())key_recn[key]=c++;}
  17. #define addkey {if((key_recn2.find(key))== key_recn2.end())key_recn2[key]=c++;}
  18. using namespace std;
  19. struct ltstr
  20. {
  21. bool operator()(const char* s1, const char* s2) const
  22. {
  23. return strcmp(s1, s2) < 0;
  24. }
  25. };
  26. struct ltint
  27. {
  28. bool operator()(const int s1, const int s2) const
  29. {
  30. return (s1<s2) ;
  31. }
  32. };
  33. int main()
  34. {
  35. char* c1=(char*)malloc(1E6L);
  36. char* c2=(char*)malloc(1E6L);
  37. char* c3=(char*)malloc(1E6L);
  38. char* c4=(char*)malloc(1E6L);
  39. int* v1=(int* )malloc(4E6L);
  40. __int64* v2;//store sum(v1);
  41. long t=clock();
  42. for (int i=0;i<1E6L;i++)
  43. {
  44. int j=i+1;
  45. c1[i]=j%4;
  46. c2[i]=j%5;
  47. c3[i]=j%67;
  48. c4[i]=(j%13?j%113:127);
  49. v1[i]=j;
  50. }
  53. map<const int, int, ltint> key_recn;
  54. map<const int, int, ltint> key_recn2;
  55. int c=0;
  56. for (int i=0;i<1E6L;i++)
  57. {
  58. int key=(c1[i]<<24)+(c2[i]<<16)+(c3[i]<<8)+c4[i];
  59. //printf("key=%d\n",key);
  60. add_key;
  62. }
  65. int *C=(int*)malloc(c*4);//store count(c1)
  66. int *k=(int*)malloc(c*4);//store key
  68. v2=(__int64*)malloc(c*8);
  70. memset(C,0,c*4); //need?
  71. memset(v2,0,c*8);
  74. for (int i=0;i<1E6L;i++)
  75. {
  76. int recn=key_recn[(c1[i]<<24)+(c2[i]<<16)+(c3[i]<<8)+c4[i]];
  77. k[recn]=(c1[i]<<24)+(c2[i]<<16)+(c3[i]<<8)+c4[i];
  78. C[recn]++;
  79. v2[recn]=v2[recn]+v1[i];
  81. }
  82. printf("%dms,recn=%d\n",clock()-t,c);
  83. t=clock();
  84. int c_2=c;
  85. for (int i=0;i<c_2;i++)
  86. {
  87. int key=(k[i]&0xFFFFFF)+((0x80+15)<<24);//set c1 null ,gid=8
  88. addkey;
  90. key=(k[i]&0xFF00FFFF)+((0x40)<<24)+(127<<16);//set c2 null,gid=4
  91. addkey;
  93. key=(k[i]&0xFFFF00FF)+((0x20)<<24)+(127<<8);//set c3 null,gid=2
  94. addkey;
  96. key=(k[i]&0xFFFFFF00)+((0x10)<<24)+(127);//set c4 null,gid=1
  97. addkey;
  99. key=(k[i]&0xFFFF00)+((0x90+15)<<24)+(127);//set c1,c4 null ,gid=9
  100. addkey;
  102. key=(k[i]&0xFF00FF00)+((0x50)<<24)+(127<<16)+(127);//set c2,c4 null,gid=5
  103. addkey;
  105. key=(k[i]&0xFFFF0000)+((0x30)<<24)+(127<<8)+(127);//set c3,c4 null,gid=3
  106. addkey;
  108. key=(k[i]&0xFF0000FF)+((0x60)<<24)+(127<<16)+(127<<8);//set c2,c3 null,gid=6
  109. addkey;
  111. key=(k[i]&0x00FF00FF)+((0xA0+15)<<24)+(127<<8);//set c1,c3 null,gid=10
  112. addkey;
  114. key=(k[i]&0xFFFF)+((0xC0+15)<<24)+(127<<16);//set c1,c2 null,gid=12
  115. addkey;
  117. key=(k[i]&0xFF00)+((0xD0+15)<<24)+(127<<16)+(127);//set c1,c2,c4 null,gid=13
  118. addkey;
  120. key=(k[i]&0xFF)+((0xE0+15)<<24)+(127<<16)+(127<<8);//set c1,c2,c3 null,gid=14
  121. addkey;
  123. key=(k[i]&0xFF0000)+((0xB0+15)<<24)+(127<<8)+(127);//set c1,c3,c4 null,gid=11
  124. addkey;
  126. key=(k[i]&0xFF000000)+((0x70)<<24)+(127<<16)+(127<<8)+(127);//set c2,c3,c4 null,gid=7
  127. addkey;
  129. key=((0xF0+15)<<24)+(127<<16)+(127<<8)+(127);//set c1,c2,c3,c4 null,gid=15
  130. addkey;
  131. }
  133. int mask[15]={
  134. 0x00FFFFFF,
  135. 0xFF00FFFF,
  136. 0xFFFF00FF,
  137. 0xFFFFFF00,
  138. 0x00FFFF00,
  139. 0xFF00FF00,
  140. 0xFFFF0000,
  141. 0xFF0000FF,
  142. 0x00FF00FF,
  143. 0x0000FFFF,
  144. 0x0000FF00,
  145. 0x000000FF,
  146. 0x00FF0000,
  147. 0xFF000000,
  148. 0x00000000};
  150. int null[15]={
  151. ((0x80+15)<<24),
  152. ((0x40)<<24)+(127<<16),
  153. ((0x20)<<24)+(127<<8),
  154. ((0x10)<<24)+(127),
  155. ((0x90+15)<<24)+(127),
  156. ((0x50)<<24)+(127<<16)+(127),
  157. ((0x30)<<24)+(127<<8)+(127),
  158. ((0x60)<<24)+(127<<16)+(127<<8),
  159. ((0xA0+15)<<24)+(127<<8),
  160. ((0xC0+15)<<24)+(127<<16),
  161. ((0xD0+15)<<24)+(127<<16)+(127),
  162. ((0xE0+15)<<24)+(127<<16)+(127<<8),
  163. ((0xB0+15)<<24)+(127<<8)+(127),
  164. ((0x70)<<24)+(127<<16)+(127<<8)+(127),
  165. ((0xF0+15)<<24)+(127<<16)+(127<<8)+(127)
  166. };
  169. int *C2=(int*)malloc((c-c_2)*4);//store count(c1)
  170. int *k2=(int*)malloc((c-c_2)*4);//store key
  171. __int64* v3=(__int64*)malloc((c-c_2)*8);
  173. memset(C2,0,(c-c_2)*4); //need?
  174. memset(v3,0,(c-c_2)*8);
  176. int y=0;
  177. for (int i=0;i<c_2;i++)
  178. {
  179. for(int j=0;j<15;j++)
  180. {
  181. int recn=key_recn2[(k[i]&mask[j])+null[j]]-c_2;
  182. k2[recn]=(k[i]&mask[j])+null[j];
  183. C2[recn]+=C[i];
  184. v3[recn]=v3[recn]+v2[i];
  186. y=recn;
  187. }
  188. }
  189. printf("%dms,recn=%d,c1=%d,c2=%d,c3=%d,c4=%d,count=%d,sum=%I64d",clock()-t,c,
  190. (k2[y]>>24),((k2[y]>>16)&255),((k2[y]>>8)&255),(k2[y]&255),C2[y],v3[y]);
  192. free(c1);
  193. free(c2);
  194. free(c3);
  195. free(c4);
  196. free(v1);
  197. free(v2);
  198. free(C);
  199. free(k);
  200. free(C2);
  201. free(k2);
  202. free(v3);
  204. return 0;
  206. }
复制代码
毋因群疑而阻独见　　毋任己意而废人言
毋私小惠而伤大体　　毋借公论以快私情
回复 支持 反对

使用道具 举报

〇〇
 楼主| 发表于 2010-3-1 01:53:19 | 显示全部楼层
  2. /* compiled by cl of vc2003*/
  3. #include <time.h>
  4. #include <stdio.h>
  5. #include <stdlib.h>
  6. #include <map>
  7. #include <iostream>
  8. #define add_key {if((key_recn.find(key))== key_recn.end())key_recn[key]=c++;}
  9. #define addkey {if((key_recn2.find(key))== key_recn2.end())key_recn2[key]=c++;}
  10. #define adkey(i) {if((key_recn2[i].find(key))== key_recn2[i].end())key_recn2[i][key]=c++;}
  11. #define sumkey(i) {f[i]++;int recn=key_recn2[i][key]-c_2;k2[recn]=key;C2[recn]+=C2[(iter->second)-c_2];v3[recn]+=v3[(iter->second)-c_2];}
  12. #define __int64 long long
  13. using namespace std;
  14. struct ltstr
  15. {
  16. bool operator()(const char* s1, const char* s2) const
  17. {
  18. return strcmp(s1, s2) < 0;
  19. }
  20. };
  21. struct ltint
  22. {
  23. bool operator()(const int s1, const int s2) const
  24. {
  25. return (s1<s2) ;
  26. }
  27. };
  28. int main()
  29. {
  30. int s[3][3]={ {4,5,67},{9,19,109},{9,97,119}};
  31. //printf("%d,%d",t[0][1],t[1][2]);
  32. //return 0;
  33. for(int r=2;r<3;r++)
  34. {
  35. long t=clock();
  36. printf("%d,%d,%d,114,%d\n",s[r][0],s[r][1],s[r][2],s[r][0]*s[r][1]*s[r][2]*114);
  37. char* c1=(char*)malloc(1E6L);
  38. char* c2=(char*)malloc(1E6L);
  39. char* c3=(char*)malloc(1E6L);
  40. char* c4=(char*)malloc(1E6L);
  41. int* v1=(int* )malloc(4E6L);
  43. printf("申请源表内存%dms\n",clock()-t);
  44. t=clock();
  46. __int64* v2;//store sum(v1);
  47. //FILE *fp=fopen("in.txt","wb");
  48. //FILE *fp=fopen("in.txt","rb");
  49. for (int i=0;i<1E6L;i++)
  50. {
  51. int j=i+1;
  52. c1[i]=j%s[r][0];
  53. c2[i]=j%s[r][1];
  54. c3[i]=j%s[r][2];
  55. c4[i]=(j%13?j%113:127);
  56. v1[i]=j;
  57. //fscanf(fp,"%d,%d,%d,%d,%d\n",&c1[i],&c2[i],&c3[i],&c4[i],&v1[i]);
  58. //fprintf(fp,"%d,%d,%d,%d,%d\n",c1[i],c2[i],c3[i],c4[i],v1[i]);
  60. }
  61. //fclose(fp);
  62. //return 1;
  63. map<const int, int, ltint> key_recn;
  64. map<const int, int, ltint> key_recn2[16];
  65. int c=0;
  66. printf("填充源表数据%dms,recn=%d,每毫秒%d行\n",clock()-t,c,(long)1E6/(clock()-t));
  67. t=clock();
  69. for (int i=0;i<1E6L;i++)
  70. {
  71. int key=(c1[i]<<24)+(c2[i]<<16)+(c3[i]<<8)+c4[i];
  72. //printf("key=%d\n",key);
  73. add_key;
  75. }
  77. printf("新建一层节点%dms,recn=%d,每毫秒读%d行,每毫秒写%d节点\n",clock()-t,c,(long)1E6/(clock()-t),c/(clock()-t));
  78. t=clock();
  79. int *C=(int*)malloc(c*4);//store count(c1)
  80. int *k=(int*)malloc(c*4);//store key
  82. v2=(__int64*)malloc(c*8);
  84. memset(C,0,c*4); //need?
  85. memset(v2,0,c*8);
  86. printf("申请一层内存%dms\n",clock()-t);
  87. t=clock();
  90. for (int i=0;i<1E6L;i++)
  91. {
  92. int recn=key_recn[(c1[i]<<24)+(c2[i]<<16)+(c3[i]<<8)+c4[i]];
  93. k[recn]=(c1[i]<<24)+(c2[i]<<16)+(c3[i]<<8)+c4[i];
  94. C[recn]++;
  95. v2[recn]=v2[recn]+v1[i];
  97. }
  98. printf("写入一层数据%dms,recn=%d,每毫秒读%d行，每毫秒写%d行\n",clock()-t,c,(long)1E6/(clock()-t),c/(clock()-t));
  99. t=clock();
  100. int c_2=c;
  101. for (int i=0;i<c_2;i++)
  102. {
  103. int key=(k[i]&0xFFFFFF)+((0x80+15)<<24);//set c1 null ,gid=8
  104. adkey(8);
  106. key=(k[i]&0xFF00FFFF)+((0x40)<<24)+(127<<16);//set c2 null,gid=4
  107. adkey(4);
  109. key=(k[i]&0xFFFF00FF)+((0x20)<<24)+(127<<8);//set c3 null,gid=2
  110. adkey(2);
  112. key=(k[i]&0xFFFFFF00)+((0x10)<<24)+(127);//set c4 null,gid=1
  113. adkey(1);
  115. /*
  116. key=(k[i]&0xFFFF00)+((0x90+15)<<24)+(127);//set c1,c4 null ,gid=9
  117. adkey(9);
  119. key=(k[i]&0xFF00FF00)+((0x50)<<24)+(127<<16)+(127);//set c2,c4 null,gid=5
  120. adkey(5);
  122. key=(k[i]&0xFFFF0000)+((0x30)<<24)+(127<<8)+(127);//set c3,c4 null,gid=3
  123. adkey(3);
  125. key=(k[i]&0xFF0000FF)+((0x60)<<24)+(127<<16)+(127<<8);//set c2,c3 null,gid=6
  126. adkey(6);
  128. key=(k[i]&0x00FF00FF)+((0xA0+15)<<24)+(127<<8);//set c1,c3 null,gid=10
  129. adkey(10);
  131. key=(k[i]&0xFFFF)+((0xC0+15)<<24)+(127<<16);//set c1,c2 null,gid=12
  132. adkey(12);
  134. key=(k[i]&0xFF00)+((0xD0+15)<<24)+(127<<16)+(127);//set c1,c2,c4 null,gid=13
  135. adkey(13);
  137. key=(k[i]&0xFF)+((0xE0+15)<<24)+(127<<16)+(127<<8);//set c1,c2,c3 null,gid=14
  138. adkey(14);
  140. key=(k[i]&0xFF0000)+((0xB0+15)<<24)+(127<<8)+(127);//set c1,c3,c4 null,gid=11
  141. adkey(11);
  143. key=(k[i]&0xFF000000)+((0x70)<<24)+(127<<16)+(127<<8)+(127);//set c2,c3,c4 null,gid=7
  144. adkey(7);
  146. key=((0xF0+15)<<24)+(127<<16)+(127<<8)+(127);//set c1,c2,c3,c4 null,gid=15
  147. adkey(15);
  148. */
  149. }
  150. map<const int, int,ltint>::iterator iter;
  151. for(int i=1;i<7;i++)
  152. {
  153. if (i==1)
  154. for(iter=key_recn2[i].begin();iter != key_recn2[i].end(); iter++)
  155. {
  156. int key=((iter->first)&0x0FFF0000)+((0x30)<<24)+(127<<8)+(127);//set c3,c4 null,gid=3
  157. adkey(3);
  159. key=((iter->first)&0xFFFF00)+((0x90+15)<<24)+(127);//set c1,c4 null ,gid=9
  160. adkey(9);
  162. key=((iter->first)&0x0F00FF00)+((0x50)<<24)+(127<<16)+(127);//set c2,c4 null,gid=5
  163. adkey(5);
  165. }
  166. if(i==2)
  167. for(iter=key_recn2[i].begin();iter != key_recn2[i].end(); iter++)
  168. {
  169. int key=((iter->first)&0x00FF00FF)+((0xA0+15)<<24)+(127<<8);//set c1,c3 null,gid=10
  170. adkey(10);
  171. }
  172. if(i==4)
  173. for(iter=key_recn2[i].begin();iter != key_recn2[i].end(); iter++)
  174. {
  175. int key=((iter->first)&0x0F0000FF)+((0x60)<<24)+(127<<16)+(127<<8);//set c2,c3 null,gid=6
  176. adkey(6);
  179. key=((iter->first)&0xFFFF)+((0xC0+15)<<24)+(127<<16);//set c1,c2 null,gid=12
  180. adkey(12);
  181. }
  182. if(i==3)
  183. for(iter=key_recn2[i].begin();iter != key_recn2[i].end(); iter++)
  184. {
  185. int key=((iter->first)&0xFF0000)+((0xB0+15)<<24)+(127<<8)+(127);//set c1,c3,c4 null,gid=11
  186. adkey(11);
  188. key=((iter->first)&0x0F000000)+((0x70)<<24)+(127<<16)+(127<<8)+(127);//set c2,c3,c4 null,gid=7
  189. adkey(7);
  190. }
  191. if(i==5)
  192. for(iter=key_recn2[i].begin();iter != key_recn2[i].end(); iter++)
  193. {
  194. int key=((iter->first)&0xFF00)+((0xD0+15)<<24)+(127<<16)+(127);//set c1,c2,c4 null,gid=13
  195. adkey(13);
  196. }
  197. if(i==6)
  198. for(iter=key_recn2[i].begin();iter != key_recn2[i].end(); iter++)
  199. {
  200. int key=((iter->first)&0xFF)+((0xE0+15)<<24)+(127<<16)+(127<<8);//set c1,c2,c3 null,gid=14
  201. adkey(14);
  202. }
  203. }
  204. int key=((0xF0+15)<<24)+(127<<16)+(127<<8)+(127);//set c1,c2,c3,c4 null,gid=15
  205. adkey(15);
  207. printf("新建后面节点%dms,recn=%d,每毫秒读%d行key,每毫秒写%d节点\n",clock()-t,c,c_2/(clock()-t),(c-c_2)/(clock()-t));
  208. for(int i=0;i<16;i++)
  209. {
  210. printf("%d,first=%0X,renc=%d\n",i,key_recn2[i].begin()->first,key_recn2[i][key_recn2[i].begin()->first]);
  211. }
  212. t=clock();
  213. int o[15]={8,4,2,1,9,5,3,6,10,12,13,14,11,7,15};
  214. int mask[15]={
  215. 0x00FFFFFF,
  216. 0xFF00FFFF,
  217. 0xFFFF00FF,
  218. 0xFFFFFF00,
  219. 0x00FFFF00,
  220. 0xFF00FF00,
  221. 0xFFFF0000,
  222. 0xFF0000FF,
  223. 0x00FF00FF,
  224. 0x0000FFFF,
  225. 0x0000FF00,
  226. 0x000000FF,
  227. 0x00FF0000,
  228. 0xFF000000,
  229. 0x00000000};
  231. int null[15]={
  232. ((0x80+15)<<24),
  233. ((0x40)<<24)+(127<<16),
  234. ((0x20)<<24)+(127<<8),
  235. ((0x10)<<24)+(127),
  236. ((0x90+15)<<24)+(127),
  237. ((0x50)<<24)+(127<<16)+(127),
  238. ((0x30)<<24)+(127<<8)+(127),
  239. ((0x60)<<24)+(127<<16)+(127<<8),
  240. ((0xA0+15)<<24)+(127<<8),
  241. ((0xC0+15)<<24)+(127<<16),
  242. ((0xD0+15)<<24)+(127<<16)+(127),
  243. ((0xE0+15)<<24)+(127<<16)+(127<<8),
  244. ((0xB0+15)<<24)+(127<<8)+(127),
  245. ((0x70)<<24)+(127<<16)+(127<<8)+(127),
  246. ((0xF0+15)<<24)+(127<<16)+(127<<8)+(127)
  247. };
  250. int *C2=(int*)malloc((c-c_2)*4);//store count(c1)
  251. int *k2=(int*)malloc((c-c_2)*4);//store key
  252. __int64* v3=(__int64*)malloc((c-c_2)*8);
  254. memset(C2,0,(c-c_2)*4); //need?
  255. memset(v3,0,(c-c_2)*8);
  256. printf("申请后面内存%dms,%d\n",clock()-t,c);
  257. t=clock();
  258. //goto f;
  259. int y=0;
  260. int f[16];
  261. memset(f,0,15*4);
  263. for (int i=0;i<c_2;i++)
  264. {
  265. for(int j=0;j<4;j++)
  266. {
  267. //int key=(k[i]&mask[j])+null[j];
  268. //if(key_recn2[o[j]].find(key)!= key_recn2[o[j]].end())
  269. {f[o[j]]++;
  270. int recn=key_recn2[o[j]][(k[i]&mask[j])+null[j]]-c_2;
  271. k2[recn]=(k[i]&mask[j])+null[j];
  272. C2[recn]+=C[i];
  273. v3[recn]+=v2[i];
  275. y=recn;
  276. }
  277. }
  278. }
  280. for(int i=1;i<7;i++)
  281. {
  282. if (i==1)
  283. for(iter=key_recn2[i].begin();iter != key_recn2[i].end(); iter++)
  284. {
  285. int key=((iter->first)&0x0FFF0000)+((0x30)<<24)+(127<<8)+(127);//set c3,c4 null,gid=3
  286. /* {f[3]++;
  287. int recn=key_recn2[3][key]-c_2;
  288. k2[recn]=key;
  289. C2[recn]+=C2[(iter->second)-c_2];
  290. v3[recn]+=v3[(iter->second)-c_2];
  291. }*/
  292. sumkey(3);
  294. key=((iter->first)&0xFFFF00)+((0x90+15)<<24)+(127);//set c1,c4 null ,gid=9
  295. sumkey(9);
  297. key=((iter->first)&0x0F00FF00)+((0x50)<<24)+(127<<16)+(127);//set c2,c4 null,gid=5
  298. sumkey(5);
  300. }
  301. if(i==2)
  302. for(iter=key_recn2[i].begin();iter != key_recn2[i].end(); iter++)
  303. {
  304. int key=((iter->first)&0x00FF00FF)+((0xA0+15)<<24)+(127<<8);//set c1,c3 null,gid=10
  305. sumkey(10);
  306. }
  308. if(i==4)
  309. for(iter=key_recn2[i].begin();iter != key_recn2[i].end(); iter++)
  310. {
  311. int key=((iter->first)&0x0F0000FF)+((0x60)<<24)+(127<<16)+(127<<8);//set c2,c3 null,gid=6
  312. sumkey(6);
  315. key=((iter->first)&0xFFFF)+((0xC0+15)<<24)+(127<<16);//set c1,c2 null,gid=12
  316. sumkey(12);
  317. }
  318. if(i==3)
  319. for(iter=key_recn2[i].begin();iter != key_recn2[i].end(); iter++)
  320. {
  321. int key=((iter->first)&0xFF0000)+((0xB0+15)<<24)+(127<<8)+(127);//set c1,c3,c4 null,gid=11
  322. sumkey(11);
  324. key=((iter->first)&0x0F000000)+((0x70)<<24)+(127<<16)+(127<<8)+(127);//set c2,c3,c4 null,gid=7
  325. sumkey(7);
  327. key=((0xF0+15)<<24)+(127<<16)+(127<<8)+(127);//set c1,c2,c3,c4 null,gid=15
  328. sumkey(15);
  329. }
  330. if(i==5)
  331. for(iter=key_recn2[i].begin();iter != key_recn2[i].end(); iter++)
  332. {
  333. int key=((iter->first)&0xFF00)+((0xD0+15)<<24)+(127<<16)+(127);//set c1,c2,c4 null,gid=13
  334. sumkey(13);
  335. }
  336. if(i==6)
  337. for(iter=key_recn2[i].begin();iter != key_recn2[i].end(); iter++)
  338. {
  339. int key=((iter->first)&0xFF)+((0xE0+15)<<24)+(127<<16)+(127<<8);//set c1,c2,c3 null,gid=14
  340. sumkey(14);
  341. }
  342. }
  344. printf("汇总后面节点%dms,recn=%d,每毫秒读%d行，每毫秒写%d行\n",clock()-t,c,c_2/(clock()-t),(c-c_2)/(clock()-t));
  345. t=clock();
  346. /*
  347. printf("%fms,recn=%d,c1=%d,c2=%d,c3=%d,c4=%d,count=%d,sum=%I64d\n",(float)(clock()-t)/CLOCKS_PER_SEC ,c,
  348. (k2[y]>>24),((k2[y]>>16)&255),((k2[y]>>8)&255),(k2[y]&255),C2[y],v3[y]);
  350. */
  351. FILE *fp=fopen("out.txt","wb");
  353. //fp=fopen("out.txt","wb");
  354. fprintf(fp,"gid,c1,c2,c3,c4,count,sum\n");
  355. for (int i=0;i<c_2;i++)
  356. {
  357. y=i;
  358. fprintf(fp,"%d,%d,%d,%d,%d,%d,%I64d\n",
  359. 0,(k[y]>>24),((k[y]>>16)&255),((k[y]>>8)&255),(k[y]&255),C[y],v2[y]);
  361. }
  363. for (int i=0;i<c-c_2;i++)
  364. {
  365. y=i;
  366. fprintf(fp,"%d,%d,%d,%d,%d,%d,%I64d\n",
  367. (16+((k2[y]>>24)>>4))&15,(k2[y]>>24)&15,((k2[y]>>16)&255),((k2[y]>>8)&255),(k2[y]&255),C2[y],v3[y]);
  369. }
  370. fclose(fp);
  371. printf("保存到文件%dms,recn=%d,每毫秒写%d行\n",clock()-t,c,c/(clock()-t));
  372. t=clock();
  374. f:
  375. free(c1);
  376. free(c2);
  377. free(c3);
  378. free(c4);
  379. free(v1);
  380. free(v2);
  381. free(C);
  382. free(k);
  383. free(C2);
  384. free(k2);
  385. free(v3);
  386. printf("释放所有动态内存%dms,%d\n",clock()-t,c);
  387. for(int i=0;i<16;i++)
  388. printf("%d,%d\n",i,f[i]);
  392. }
  393. return 0;
  394. }
复制代码
毋因群疑而阻独见　　毋任己意而废人言
毋私小惠而伤大体　　毋借公论以快私情
回复 支持 反对

使用道具 举报

〇〇
 楼主| 发表于 2010-3-2 12:48:04 | 显示全部楼层
rollup的选择还是有讲究的
毋因群疑而阻独见　　毋任己意而废人言
毋私小惠而伤大体　　毋借公论以快私情
回复 支持 反对

使用道具 举报

〇〇
 楼主| 发表于 2010-3-2 21:34:47 | 显示全部楼层
  2. /* compiled by cl of vc2003*/
  3. #include <time.h>
  4. #include <stdio.h>
  5. #include <stdlib.h>
  6. #include <map>
  7. #include <iostream>
  8. #define add_key {if((key_recn.find(key))== key_recn.end())key_recn[key]=c++;}
  9. #define addkey {if((key_recn2.find(key))== key_recn2.end())key_recn2[key]=c++;}
  10. #define addskey {if((iter=key_recn.find(key))== key_recn.end()){key_recn[key]=c++;int recn=c;k[recn]=key;C[recn]++;v2[recn]+=v1[i];}\
  11. else{int recn=iter->second;k[recn]=key;C[recn]++;v2[recn]+=v1[i];}}
  12. /*int recn=key_recn[(c1[i]<<24)+(c2[i]<<16)+(c3[i]<<8)+c4[i]];
  13. k[recn]=(c1[i]<<24)+(c2[i]<<16)+(c3[i]<<8)+c4[i];
  14. C[recn]++;
  15. v2[recn]=v2[recn]+v1[i];*/
  17. #define adkey(i) {if((key_recn2[i].find(key))== key_recn2[i].end())key_recn2[i][key]=c++;}
  18. #define sumkey(i) {f[i]++;int recn=key_recn2[i][key]-c_2;k2[recn]=key;C2[recn]+=C2[(iter->second)-c_2];v3[recn]+=v3[(iter->second)-c_2];}
  19. #define adsmkey(j) {if((iter=key_recn2[j].find(key))== key_recn2[j].end()){key_recn2[j][key]=c++;f[j]++;int recn=c-c_2;k2[recn]=key;C2[recn]+=C[i];v3[recn]+=v2[i];}\
  20. else{f[j]++;int recn=/*key_recn2[j][key]-c_2*/iter->second-c_2;k2[recn]=key;C2[recn]+=C[i];v3[recn]+=v2[i];}}
  22. #define __int64 long long
  23. using namespace std;
  24. struct ltstr
  25. {
  26. bool operator()(const char* s1, const char* s2) const
  27. {
  28. return strcmp(s1, s2) < 0;
  29. }
  30. };
  31. struct ltint
  32. {
  33. bool operator()(const int s1, const int s2) const
  34. {
  35. return (s1<s2) ;
  36. }
  37. };
  38. int main()
  39. {
  40. int dsize=0;
  41. int s[3][3]={ {4,5,67},{9,109,19},{9,97,119}};
  42. //printf("%d,%d",t[0][1],t[1][2]);
  43. //return 0;
  44. for(int it=1;it<2;it++)
  45. for(int r=2;r<3;r++)
  46. {
  47. long t=clock();
  48. long t1=t;
  49. printf("%d,%d,%d,114,%d\n",s[r][0],s[r][1],s[r][2],s[r][0]*s[r][1]*s[r][2]*114);
  50. char* c1=(char*)malloc(1E6L);
  51. char* c2=(char*)malloc(1E6L);
  52. char* c3=(char*)malloc(1E6L);
  53. char* c4=(char*)malloc(1E6L);
  54. int* v1=(int* )malloc(4E6L);
  56. printf("申请源表内存%dms\n",clock()-t);
  57. t=clock();
  59. __int64* v2;//store sum(v1);
  60. //FILE *fp=fopen("in.txt","wb");
  61. //FILE *fp=fopen("in.txt","rb");
  62. for (int i=0;i<1E6L;i++)
  63. {
  64. int j=i+1;
  65. c1[i]=j%s[r][0];
  66. c2[i]=j%s[r][1];
  67. c3[i]=j%s[r][2];
  68. c4[i]=(j%13?j%113:127);
  69. v1[i]=j;
  70. //fscanf(fp,"%d,%d,%d,%d,%d\n",&c1[i],&c2[i],&c3[i],&c4[i],&v1[i]);
  71. //fprintf(fp,"%d,%d,%d,%d,%d\n",c1[i],c2[i],c3[i],c4[i],v1[i]);
  73. }
  74. //fclose(fp);
  75. //return 1;
  76. map<const int, int, ltint> key_recn;
  77. map<const int, int, ltint> key_recn2[16];
  78. map<const int, int,ltint>::iterator iter;
  80. int c=1E6L;
  81. printf("填充源表数据%dms,recn=%d,每毫秒%d行\n",clock()-t,c,(long)1E6/(clock()-t));
  82. t=clock();
  84. int *C=(int*)malloc(c*4);//store count(c1)
  85. int *k=(int*)malloc(c*4);//store key
  87. v2=(__int64*)malloc(c*8);
  89. memset(C,0,c*4); //need?
  90. memset(v2,0,c*8);
  91. printf("申请一层内存%dms\n",clock()-t);
  92. t=clock();
  93. c=0;
  94. for (int i=0;i<1E6L;i++)
  95. {
  96. int key=(c1[i]<<24)+(c2[i]<<16)+(c3[i]<<8)+c4[i];
  97. //printf("key=%d\n",key);
  98. addskey;
  100. }
  102. printf("新建一层节点%dms,recn=%d,每毫秒读%d行,每毫秒写%d节点\n",clock()-t,c,(long)1E6/(clock()-t),c/(clock()-t));
  103. t=clock();
  105. /*
  106. for (int i=0;i<1E6L;i++)
  107. {
  108. int recn=key_recn[(c1[i]<<24)+(c2[i]<<16)+(c3[i]<<8)+c4[i]];
  109. k[recn]=(c1[i]<<24)+(c2[i]<<16)+(c3[i]<<8)+c4[i];
  110. C[recn]++;
  111. v2[recn]=v2[recn]+v1[i];
  113. }
  114. */
  115. //printf("写入一层数据%dms,recn=%d,每毫秒读%d行，每毫秒写%d行\n",clock()-t,c,(long)1E6/(clock()-t),c/(clock()-t));
  116. t=clock();
  117. int c_2=c;
  118. int o[15]={8,4,2,1,9,5,3,6,10,12,13,14,11,7,15};
  119. int mask[15]={
  120. 0x00FFFFFF,
  121. 0xFF00FFFF,
  122. 0xFFFF00FF,
  123. 0xFFFFFF00,
  124. 0x00FFFF00,
  125. 0xFF00FF00,
  126. 0xFFFF0000,
  127. 0xFF0000FF,
  128. 0x00FF00FF,
  129. 0x0000FFFF,
  130. 0x0000FF00,
  131. 0x000000FF,
  132. 0x00FF0000,
  133. 0xFF000000,
  134. 0x00000000};
  136. int null[15]={
  137. ((0x80+15)<<24),
  138. ((0x40)<<24)+(127<<16),
  139. ((0x20)<<24)+(127<<8),
  140. ((0x10)<<24)+(127),
  141. ((0x90+15)<<24)+(127),
  142. ((0x50)<<24)+(127<<16)+(127),
  143. ((0x30)<<24)+(127<<8)+(127),
  144. ((0x60)<<24)+(127<<16)+(127<<8),
  145. ((0xA0+15)<<24)+(127<<8),
  146. ((0xC0+15)<<24)+(127<<16),
  147. ((0xD0+15)<<24)+(127<<16)+(127),
  148. ((0xE0+15)<<24)+(127<<16)+(127<<8),
  149. ((0xB0+15)<<24)+(127<<8)+(127),
  150. ((0x70)<<24)+(127<<16)+(127<<8)+(127),
  151. ((0xF0+15)<<24)+(127<<16)+(127<<8)+(127)
  152. };
  154. dsize=c_2;
  155. int *C2=(int*)malloc((c_2)*4);//store count(c1)
  156. int *k2=(int*)malloc((c_2)*4);//store key
  157. __int64* v3=(__int64*)malloc((c_2)*8);
  159. memset(C2,0,(c_2)*4); //need?
  160. memset(v3,0,(c_2)*8);
  161. printf("申请后面内存%dms,%d\n",clock()-t,c_2);
  162. t=clock();
  163. //goto f;
  164. int y=0;
  165. int f[16];
  166. memset(f,0,16*4);
  167. for (int i=0;i<c_2;i++)
  168. {
  169. if(dsize<c-c_2+4)
  170. {
  171. dsize+=dsize;
  172. C2=(int*)realloc(C2,(dsize)*4);//store count(c1)
  173. k2=(int*)realloc(k2,(dsize)*4);//store key
  174. v3=(__int64*)realloc(v3,(dsize)*8);
  176. //memset(C2,0,(c_2)*4); //need?
  177. //memset(v3,0,(c_2)*8);
  178. }
  179. int key=(k[i]&0xFFFFFF)+((0x80+15)<<24);//set c1 null ,gid=8
  180. adsmkey(8);
  182. key=(k[i]&0xFF00FFFF)+((0x40)<<24)+(127<<16);//set c2 null,gid=4
  183. adsmkey(4);
  185. key=(k[i]&0xFFFF00FF)+((0x20)<<24)+(127<<8);//set c3 null,gid=2
  186. adsmkey(2);
  188. key=(k[i]&0xFFFFFF00)+((0x10)<<24)+(127);//set c4 null,gid=1
  189. adsmkey(1);
  191. /*
  192. for(int j=0;j<4;j++)
  193. {
  194. //int key=(k[i]&mask[j])+null[j];
  195. //if(key_recn2[o[j]].find(key)!= key_recn2[o[j]].end())
  196. {f[o[j]]++;
  197. int recn=key_recn2[o[j]][(k[i]&mask[j])+null[j]]-c_2;
  198. k2[recn]=(k[i]&mask[j])+null[j];
  199. C2[recn]+=C[i];
  200. v3[recn]+=v2[i];
  202. y=recn;
  203. }
  205. }
  207. key=(k[i]&0xFFFF00)+((0x90+15)<<24)+(127);//set c1,c4 null ,gid=9
  208. adkey(9);
  210. key=(k[i]&0xFF00FF00)+((0x50)<<24)+(127<<16)+(127);//set c2,c4 null,gid=5
  211. adkey(5);
  213. key=(k[i]&0xFFFF0000)+((0x30)<<24)+(127<<8)+(127);//set c3,c4 null,gid=3
  214. adkey(3);
  216. key=(k[i]&0xFF0000FF)+((0x60)<<24)+(127<<16)+(127<<8);//set c2,c3 null,gid=6
  217. adkey(6);
  219. key=(k[i]&0x00FF00FF)+((0xA0+15)<<24)+(127<<8);//set c1,c3 null,gid=10
  220. adkey(10);
  222. key=(k[i]&0xFFFF)+((0xC0+15)<<24)+(127<<16);//set c1,c2 null,gid=12
  223. adkey(12);
  225. key=(k[i]&0xFF00)+((0xD0+15)<<24)+(127<<16)+(127);//set c1,c2,c4 null,gid=13
  226. adkey(13);
  228. key=(k[i]&0xFF)+((0xE0+15)<<24)+(127<<16)+(127<<8);//set c1,c2,c3 null,gid=14
  229. adkey(14);
  231. key=(k[i]&0xFF0000)+((0xB0+15)<<24)+(127<<8)+(127);//set c1,c3,c4 null,gid=11
  232. adkey(11);
  234. key=(k[i]&0xFF000000)+((0x70)<<24)+(127<<16)+(127<<8)+(127);//set c2,c3,c4 null,gid=7
  235. adkey(7);
  237. key=((0xF0+15)<<24)+(127<<16)+(127<<8)+(127);//set c1,c2,c3,c4 null,gid=15
  238. adkey(15);
  239. */
  240. }
  243. //map<const int, int,ltint>::iterator iter;
  244. switch(it)
  245. {
  246. case 1:
  247. for(int i=1;i<7;i++)
  248. {
  249. if (i==1)
  250. for(iter=key_recn2[i].begin();iter != key_recn2[i].end(); iter++)
  251. {
  252. int key=((iter->first)&0x0FFF0000)+((0x30)<<24)+(127<<8)+(127);//set c3,c4 null,gid=3
  253. adkey(3);
  255. key=((iter->first)&0xFFFF00)+((0x90+15)<<24)+(127);//set c1,c4 null ,gid=9
  256. adkey(9);
  258. key=((iter->first)&0x0F00FF00)+((0x50)<<24)+(127<<16)+(127);//set c2,c4 null,gid=5
  259. adkey(5);
  261. }
  262. if(i==2)
  263. for(iter=key_recn2[i].begin();iter != key_recn2[i].end(); iter++)
  264. {
  265. int key=((iter->first)&0x00FF00FF)+((0xA0+15)<<24)+(127<<8);//set c1,c3 null,gid=10
  266. adkey(10);
  267. }
  268. if(i==4)
  269. for(iter=key_recn2[i].begin();iter != key_recn2[i].end(); iter++)
  270. {
  271. int key=((iter->first)&0x0F0000FF)+((0x60)<<24)+(127<<16)+(127<<8);//set c2,c3 null,gid=6
  272. adkey(6);
  275. key=((iter->first)&0xFFFF)+((0xC0+15)<<24)+(127<<16);//set c1,c2 null,gid=12
  276. adkey(12);
  277. }
  278. if(i==3)
  279. for(iter=key_recn2[i].begin();iter != key_recn2[i].end(); iter++)
  280. {
  281. int key=((iter->first)&0xFF0000)+((0xB0+15)<<24)+(127<<8)+(127);//set c1,c3,c4 null,gid=11
  282. adkey(11);
  284. key=((iter->first)&0x0F000000)+((0x70)<<24)+(127<<16)+(127<<8)+(127);//set c2,c3,c4 null,gid=7
  285. adkey(7);
  286. }
  287. if(i==5)
  288. for(iter=key_recn2[i].begin();iter != key_recn2[i].end(); iter++)
  289. {
  290. int key=((iter->first)&0xFF00)+((0xD0+15)<<24)+(127<<16)+(127);//set c1,c2,c4 null,gid=13
  291. adkey(13);
  292. }
  293. if(i==6)
  294. for(iter=key_recn2[i].begin();iter != key_recn2[i].end(); iter++)
  295. {
  296. int key=((iter->first)&0xFF)+((0xE0+15)<<24)+(127<<16)+(127<<8);//set c1,c2,c3 null,gid=14
  297. adkey(14);
  298. }
  299. }
  300. break;
  301. case 2:
  302. for(int i=1;i<7;i++)
  303. {
  304. if (i==1)
  305. for(iter=key_recn2[i].begin();iter != key_recn2[i].end(); iter++)
  306. {
  308. }
  309. if(i==2)
  310. for(iter=key_recn2[i].begin();iter != key_recn2[i].end(); iter++)
  311. {
  312. int key=((iter->first)&0x0FFF0000)+((0x30)<<24)+(127<<8)+(127);//set c3,c4 null,gid=3
  313. adkey(3);
  315. key=((iter->first)&0xFFFF00)+((0x90+15)<<24)+(127);//set c1,c4 null ,gid=9
  316. adkey(9);
  318. key=((iter->first)&0x0F00FF00)+((0x50)<<24)+(127<<16)+(127);//set c2,c4 null,gid=5
  319. adkey(5);
  321. key=((iter->first)&0x00FF00FF)+((0xA0+15)<<24)+(127<<8);//set c1,c3 null,gid=10
  322. adkey(10);
  323. }
  324. if(i==4)
  325. for(iter=key_recn2[i].begin();iter != key_recn2[i].end(); iter++)
  326. {
  327. int key=((iter->first)&0x0F0000FF)+((0x60)<<24)+(127<<16)+(127<<8);//set c2,c3 null,gid=6
  328. adkey(6);
  331. key=((iter->first)&0xFFFF)+((0xC0+15)<<24)+(127<<16);//set c1,c2 null,gid=12
  332. adkey(12);
  333. }
  334. if(i==3)
  335. for(iter=key_recn2[i].begin();iter != key_recn2[i].end(); iter++)
  336. {
  337. int key=((iter->first)&0xFF0000)+((0xB0+15)<<24)+(127<<8)+(127);//set c1,c3,c4 null,gid=11
  338. adkey(11);
  340. key=((iter->first)&0x0F000000)+((0x70)<<24)+(127<<16)+(127<<8)+(127);//set c2,c3,c4 null,gid=7
  341. adkey(7);
  342. }
  343. if(i==5)
  344. for(iter=key_recn2[i].begin();iter != key_recn2[i].end(); iter++)
  345. {
  346. int key=((iter->first)&0xFF00)+((0xD0+15)<<24)+(127<<16)+(127);//set c1,c2,c4 null,gid=13
  347. adkey(13);
  348. }
  349. if(i==6)
  350. for(iter=key_recn2[i].begin();iter != key_recn2[i].end(); iter++)
  351. {
  352. int key=((iter->first)&0xFF)+((0xE0+15)<<24)+(127<<16)+(127<<8);//set c1,c2,c3 null,gid=14
  353. adkey(14);
  354. }
  355. }
  356. break;
  357. default:break;
  358. }
  359. int key=((0xF0+15)<<24)+(127<<16)+(127<<8)+(127);//set c1,c2,c3,c4 null,gid=15
  360. adkey(15);
  362. printf("新建后面节点%dms,recn=%d,每毫秒读%d行key,每毫秒写%d节点\n",clock()-t,c,c_2/(clock()-t),(c-c_2)/(clock()-t));
  363. for(int i=0;i<16;i++)
  364. {
  365. printf("%d,first=%0X,renc=%d\n",i,key_recn2[i].begin()->first,key_recn2[i][key_recn2[i].begin()->first]);
  366. }
  367. t=clock();
  369. /*
  370. for (int i=0;i<c_2;i++)
  371. {
  372. for(int j=0;j<4;j++)
  373. {
  374. //int key=(k[i]&mask[j])+null[j];
  375. //if(key_recn2[o[j]].find(key)!= key_recn2[o[j]].end())
  376. {f[o[j]]++;
  377. int recn=key_recn2[o[j]][(k[i]&mask[j])+null[j]]-c_2;
  378. k2[recn]=(k[i]&mask[j])+null[j];
  379. C2[recn]+=C[i];
  380. v3[recn]+=v2[i];
  382. y=recn;
  383. }
  384. }
  385. }
  386. */
  387. for(int i=1;i<7;i++)
  388. {
  389. if (i==1)
  390. for(iter=key_recn2[i].begin();iter != key_recn2[i].end(); iter++)
  391. {
  392. int key=((iter->first)&0x0FFF0000)+((0x30)<<24)+(127<<8)+(127);//set c3,c4 null,gid=3
  393. /* {f[3]++;
  394. int recn=key_recn2[3][key]-c_2;
  395. k2[recn]=key;
  396. C2[recn]+=C2[(iter->second)-c_2];
  397. v3[recn]+=v3[(iter->second)-c_2];
  398. }*/
  399. sumkey(3);
  401. key=((iter->first)&0xFFFF00)+((0x90+15)<<24)+(127);//set c1,c4 null ,gid=9
  402. sumkey(9);
  404. key=((iter->first)&0x0F00FF00)+((0x50)<<24)+(127<<16)+(127);//set c2,c4 null,gid=5
  405. sumkey(5);
  407. }
  408. if(i==2)
  409. for(iter=key_recn2[i].begin();iter != key_recn2[i].end(); iter++)
  410. {
  411. int key=((iter->first)&0x00FF00FF)+((0xA0+15)<<24)+(127<<8);//set c1,c3 null,gid=10
  412. sumkey(10);
  413. }
  415. if(i==4)
  416. for(iter=key_recn2[i].begin();iter != key_recn2[i].end(); iter++)
  417. {
  418. int key=((iter->first)&0x0F0000FF)+((0x60)<<24)+(127<<16)+(127<<8);//set c2,c3 null,gid=6
  419. sumkey(6);
  422. key=((iter->first)&0xFFFF)+((0xC0+15)<<24)+(127<<16);//set c1,c2 null,gid=12
  423. sumkey(12);
  424. }
  425. if(i==3)
  426. for(iter=key_recn2[i].begin();iter != key_recn2[i].end(); iter++)
  427. {
  428. int key=((iter->first)&0xFF0000)+((0xB0+15)<<24)+(127<<8)+(127);//set c1,c3,c4 null,gid=11
  429. sumkey(11);
  431. key=((iter->first)&0x0F000000)+((0x70)<<24)+(127<<16)+(127<<8)+(127);//set c2,c3,c4 null,gid=7
  432. sumkey(7);
  434. key=((0xF0+15)<<24)+(127<<16)+(127<<8)+(127);//set c1,c2,c3,c4 null,gid=15
  435. sumkey(15);
  436. }
  437. if(i==5)
  438. for(iter=key_recn2[i].begin();iter != key_recn2[i].end(); iter++)
  439. {
  440. int key=((iter->first)&0xFF00)+((0xD0+15)<<24)+(127<<16)+(127);//set c1,c2,c4 null,gid=13
  441. sumkey(13);
  442. }
  443. if(i==6)
  444. for(iter=key_recn2[i].begin();iter != key_recn2[i].end(); iter++)
  445. {
  446. int key=((iter->first)&0xFF)+((0xE0+15)<<24)+(127<<16)+(127<<8);//set c1,c2,c3 null,gid=14
  447. sumkey(14);
  448. }
  449. }
  451. printf("汇总后面节点%dms,recn=%d,每毫秒读%d行，每毫秒写%d行\n",clock()-t,c,c_2/(clock()-t),(c-c_2)/(clock()-t));
  452. t=clock();
  453. /*
  454. printf("%fms,recn=%d,c1=%d,c2=%d,c3=%d,c4=%d,count=%d,sum=%I64d\n",(float)(clock()-t)/CLOCKS_PER_SEC ,c,
  455. (k2[y]>>24),((k2[y]>>16)&255),((k2[y]>>8)&255),(k2[y]&255),C2[y],v3[y]);
  458. FILE *fp=fopen("out.txt","wb");
  460. //fp=fopen("out.txt","wb");
  461. fprintf(fp,"gid,c1,c2,c3,c4,count,sum\n");
  462. for (int i=0;i<c_2;i++)
  463. {
  464. y=i;
  465. fprintf(fp,"%d,%d,%d,%d,%d,%d,%I64d\n",
  466. 0,(k[y]>>24),((k[y]>>16)&255),((k[y]>>8)&255),(k[y]&255),C[y],v2[y]);
  468. }
  470. for (int i=0;i<c-c_2;i++)
  471. {
  472. y=i;
  473. fprintf(fp,"%d,%d,%d,%d,%d,%d,%I64d\n",
  474. (16+((k2[y]>>24)>>4))&15,(k2[y]>>24)&15,((k2[y]>>16)&255),((k2[y]>>8)&255),(k2[y]&255),C2[y],v3[y]);
  476. }
  477. fclose(fp);
  478. printf("保存到文件%dms,recn=%d,每毫秒写%d行\n",clock()-t,c,c/(clock()-t));
  479. t=clock();
  480. */
  481. f:
  482. free(c1);
  483. free(c2);
  484. free(c3);
  485. free(c4);
  486. free(v1);
  487. free(v2);
  488. free(C);
  489. free(k);
  490. free(C2);
  491. free(k2);
  492. free(v3);
  493. printf("释放所有动态内存%dms,%d\n",clock()-t,c);
  494. for(int i=0;i<16;i++)
  495. printf("%d,%d\n",i,f[i]);
  497. printf("总时间%dms,%d\n",clock()-t1,c);
  500. }
  501. return 0;
  502. }
复制代码
毋因群疑而阻独见　　毋任己意而废人言
毋私小惠而伤大体　　毋借公论以快私情
回复 支持 反对

使用道具 举报

〇〇
 楼主| 发表于 2010-3-3 13:12:14 | 显示全部楼层
源代码 mapp.rar (2.78 KB, 下载次数: 0) 比较结果 comp.gif
毋因群疑而阻独见　　毋任己意而废人言
毋私小惠而伤大体　　毋借公论以快私情
回复 支持 反对

使用道具 举报

wayne
发表于 2010-3-3 13:14:13 | 显示全部楼层
酒向知己饮 可惜我不会饮酒，
毋因群疑而阻独见　　毋任己意而废人言
毋私小惠而伤大体　　毋借公论以快私情
回复 支持 反对

使用道具 举报

〇〇
 楼主| 发表于 2010-3-3 16:40:52 | 显示全部楼层
所有的代码

maps.cpp.rar

193.61 KB, 下载次数: 0, 下载积分: 金币 -1 枚, 经验 1 点, 下载 1 次
毋因群疑而阻独见　　毋任己意而废人言
毋私小惠而伤大体　　毋借公论以快私情
回复 支持 反对

使用道具 举报

〇〇
 楼主| 发表于 2010-3-9 16:15:48 | 显示全部楼层
map9d

map9d.rar

8.53 KB, 下载次数: 0, 下载积分: 金币 -1 枚, 经验 1 点, 下载 1 次
毋因群疑而阻独见　　毋任己意而废人言
毋私小惠而伤大体　　毋借公论以快私情
回复 支持 反对

使用道具 举报

返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

小黑屋|手机版|数学研发网 ( 苏ICP备07505100号 )

GMT+8, 2024-11-22 02:05 , Processed in 0.036402 second(s), 20 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

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