设为首页收藏本站
切换到宽版

 找回密码
 注册
搜索
批处理之家»论坛 › 批处理专区: BAT CMD DOS › 开源原创工具 › 数学计算工具 i 的重构版new i
[新手上路]批处理新手入门导读[视频教程]批处理基础视频教程[视频教程]VBS基础视频教程[批处理精品]批处理版照片整理器
[批处理精品]纯批处理备份&还原驱动[批处理精品]CMD命令50条不能说的秘密[在线下载]第三方命令行工具[在线帮助]VBScript / JScript 在线参考
返回列表 发新帖
查看: 21125|回复: 6

数学计算工具 i 的重构版new i

[复制链接]
happy886rr
发表于 2017-2-23 20:56:13 | 显示全部楼层 |阅读模式
NEW I.EXE  版本1.0

摘要：
=============================================================================
数学计算工具I.EXE的重构版，更换了新的内核和计算核心，改用自主研发的静逆波兰
解析器，多处if改用case实现，抛弃低端的strcmp函数， 重写了程序主体架构，完善
了之前不安全的代码。采用双引擎素数筛法，可以快捷获取一万亿以内的素数。
new i 工具的用法与旧 i 已经不同，旧版本宣布作废，今后将只对new i 做升级维护。
=============================================================================

下载链接: https://pan.baidu.com/s/1tPFRm3Mck7XYBXWME2pD4w?pwd=9w8g

用法：(开关支持近似输入法，比如 i /get-permutation 可以直接输入成 i /get-pe )
-----------------------------------------------------------------------------
  i [options] {arguments}...
-----------------------------------------------------------------------------
    /help Show help information

    /get-fib  [n]
    /get-fact [n]
    /get-pi   [Accuracy]
    /get-e    [Accuracy]
    /get-combination [m] [n]
    /get-permutation [m] [n]
    /get-datedifference [start_day] [end_day]
    /get-hash [string]
    /get-result [expression]
    /get-statistics [arrays]
    /get-solution [coefficients]
    /get-elements
    /get-mytax [principal] [interest rate] [number of compound interest]
    /get-interest [annuity] [interest rate] [interest rate period]
    /get-annuity [income] [three insurance a gold]

    /pri-list  [m] [n]
    /pri-count [max]
    /pri-index [m] [n]
    /pri-check [prime]
    /pri-sum   [max]
-----------------------------------------------------------------------------

示例：
-----------------------------------------------------------------------------
REM 计算字符串哈希值
i /get-hash "i"
Crc32 : e66c3671
Md5   : 865c0c0b4ab0e063e5caa3387c1a8741
Sha1  : 042dc4512fa3d391c5170cf3aa61e6a638f84342
Sha256: de7d1b721a1e0632b7cf04edf5032c8ecffa9f9a08492152b926f1a5a7e765d7
-----------------------------------------------------------------------------
REM 计算组合数 C(33,6)
i /get-comb 33 6
=1107568
-----------------------------------------------------------------------------
REM 计算排列数 P(33,6)
i /get-perm 33 3
=32736
-----------------------------------------------------------------------------
REM 科学计算  5*(arctan(cos(sin(ln(lg(2.71828))))))
i /get-result 5*(arctan(cos(sin(ln(lg(2.71828))))))
=3.178987410777
-----------------------------------------------------------------------------
REM 素数切片
i /pri-list 100000000 100000500
┌────────────────────────┐
100000007 100000037 100000039 100000049 100000073
100000081 100000123 100000127 100000193 100000213
100000217 100000223 100000231 100000237 100000259
100000267 100000279 100000357 100000379 100000393
100000399 100000421 100000429 100000463 100000469
100000471 100000493
└────────────────────────┘
-----------------------------------------------------------------------------
REM 第一千万个素数
i /pri-index 10000000
=179424673
-----------------------------------------------------------------------------
REM 素性检测
i /pri-check 100000357
100000357 是素数
-----------------------------------------------------------------------------
REM 计算30亿内的素数和
i /pri-sum 3000000000
=211169315361937456
-----------------------------------------------------------------------------
REM 统计300亿内的素数个数
i /pri-count 30000000000
在30000000000以内共有1300005926个素数
-----------------------------------------------------------------------------
RME 解高次方程
i /get-solu 3 0 5 -6 -9
一元四次方程3x^4+0x^3+5x^2-6x-9=0的解为：
┌────────────────────────┐
  共轭复根：-0.252438186547+1.696390660707i
            -0.252438186547-1.696390660707i
  实数解为： 1.293411063722
            -0.788534690627
└────────────────────────┘
-----------------------------------------------------------------------------
RME 计算日期差
i /get-dat 19951207 20170223
=7749
-----------------------------------------------------------------------------
REM 其他功能请自行体验，不再详叙。
-----------------------------------------------------------------------------

备注：
-----------------------------------------------------------------------------
常数类
        pi    3.1415926535897932
        e     2.7182818284590452       

通用类
        rand  随机数
        round 四舍五入
        int   取整
        ceil  向上舍入
        floor 向下舍入
        abs   绝对值
        sqrt  开方
        lg    常用对数，以10为底
        ln    自然对数
        exp   e的次幂
        gamma 伽玛函数
        torad 度转弧度
        +     加
        -     减
        *     乘
        /     除
        %     取余数
        ^     次方
        !     阶乘

三角函数类
        sin、cos、tan   
        arcsin、arccos、arctan

双曲函数类
        sinh、cosh、tanh
        arcsinh、arccosh、arctanh
-----------------------------------------------------------------------------

英译：
-----------------------------------------------------------------------------
CONSOLE CALCULATE TOOL, COPYRIGHT@2017~2019 BY HAPPY, VERSION 1.0
-----------------------------------------------------------------------------
  i [options] {arguments}...
-----------------------------------------------------------------------------
    /help Show help information

    /get-fib  [n]
    /get-fact [n]
    /get-pi   [Accuracy]
    /get-e    [Accuracy]
    /get-combination [m] [n]
    /get-permutation [m] [n]
    /get-datedifference [start_day] [end_day]
    /get-hash [string]
    /get-result [expression]
    /get-statistics [arrays]
    /get-solution [coefficients]
    /get-elements
    /get-mytax [principal] [interest rate] [number of compound interest]
    /get-interest [annuity] [interest rate] [interest rate period]
    /get-annuity [income] [three insurance a gold]

    /pri-list  [m] [n]
    /pri-count [max]
    /pri-index [m] [n]
    /pri-check [prime]
    /pri-sum   [max]
-----------------------------------------------------------------------------
2017-02-23


核心源码：
  1. /*
  2.         CONSOLE CALCULATE TOOL, COPYRIGHT@2017~2019 BY HAPPY, VERSION 1.0
  3.         THE NEW I.EXE
  4.         LIN: EQUATION.LIB  HASH.LIB  REVPOLISH.LIB SMALLFUN.LIB PRIME.LIB
  5. */

  7. #include    <math.h>
  8. #include   <stdio.h>
  9. #include <stdbool.h>
  10. #include  <stdlib.h>
  11. #include  <string.h>

  13. //定义帮助说明//////////////////////////////////////////////////////
  14. #define HELP_INFORMATION "\
  15. CONSOLE CACLUTE TOOL, COPYRIGHT@2017~2018 BY HAPPY, VERSION 1.0  \n\
  16. -----------------------------------------------------------------\n\
  17.   i [option] {parameters}\n\
  18. -----------------------------------------------------------------\n\
  19.     /get-fib  [n]\n\
  20.     /get-fact [n]\n\
  21.     /get-pi   [Accuracy]\n\
  22.     /get-e    [Accuracy]\n\
  23.     /get-combination [m] [n]\n\
  24.     /get-permutation [m] [n]\n\
  25.     /get-datedifference [start_day] [end_day]\n\
  26.     /get-hash [string]\n\
  27.     /get-result [expression]\n\
  28.     /get-statistics [arrays]\n\
  29.     /get-solution [coefficients]\n\
  30.     /get-elements\n\
  31.     /get-mytax [principal] [interest rate] [number of compound interest]\n\
  32.     /get-interest [annuity] [interest rate] [interest rate period]\n\
  33.     /get-annuity [income] [three insurance a gold]\n\
  34. \n\
  35.     /pri-list  [m] [n]\n\
  36.     /pri-count [max]\n\
  37.     /pri-index [m] [n]\n\
  38.     /pri-check [prime]\n\
  39.     /pri-sum   [max]\n\
  40. -----------------------------------------------------------------\n\
  41. 02-23-2017\n"
  42. ////////////////////////////////////////////////////////////////////
  43. #define CHEMICAL_ELEMENTS "------Chemical_elements------\n\
  44.     ⅠA                                                                                                     0   \n\
  45.   ┌──┐                                                                                              ┌──┐\n\
  46.   │1　H│                                                                                              │2 He│\n\
  47. 1│ 氢 │                                                                                              │ 氦 │\n\
  48.   │1.01│ ⅡA                                                               ⅢA   ⅣA   ⅤA   ⅥA   ⅦA│4.00│\n\
  49.   ├──┼──┐                                                          ┌──┬──┬──┬──┬──┼──┤\n\
  50.   │3 Li│4 Be│                                                          │5  B│6  C│7  N│8  O│9  F│10Ne│\n\
  51. 2│ 锂 │ 铍 │                                                          │ 硼 │ 碳 │ 氮 │ 氧 │ 氟 │ 氖 │\n\
  52.   │6.94│9.01│                                                          │10.8│12.0│14.0│16.0│19.0│20.2│\n\
  53.   ├──┼──┤                                                          ├──┼──┼──┼──┼──┼──┤\n\
  54.   │11Li│12Be│                                                          │13Al│14Si│15 P│16 S│17Cl│18Ar│\n\
  55. 3│ 钠 │ 镁 │                                                          │ 铝 │ 硅 │ 磷 │ 硫 │ 氯 │ 氩 │\n\
  56.   │23.0│24.3│ ⅢB   ⅣB   ⅤB   ⅥB   ⅦB         Ⅷ         ⅠB   ⅡB │27.0│28.1│31.0│32.1│35.5│40.0│\n\
  57.   ├──┼──┼──┬──┬──┬──┬──┬──┬──┬──┬──┬──┼──┼──┼──┼──┼──┼──┤\n\
  58.   │19 K│20Ca│21Sc│22Ti│23 V│24Cr│25Mn│26Fe│27Co│28Ni│29Cu│30Zn│31Ga│32Ge│33As│34Se│35Br│36Kr│\n\
  59. 4│ 钾 │ 钙 │ 钪 │ 钛 │ 钒 │ 铬 │ 锰 │ 铁 │ 钴 │ 镍 │ 铜 │ 锌 │ 镓 │ 锗 │ 砷 │ 硒 │ 溴 │ 氪 │\n\
  60.   │39.1│40.1│45.0│47.8│51.0│52.0│54.9│55.9│58.9│58.7│63.6│65.4│69.8│72.6│74.9│79.0│79.9│83.8│\n\
  61.   ├──┼──┼──┼──┼──┼──┼──┼──┼──┼──┼──┼──┼──┼──┼──┼──┼──┼──┤\n\
  62.   │37Rb│38Sr│39 Y│40Zr│41Nb│42Mo│43Tc│44Ru│45Rh│46Pd│47Ag│48Cd│49In│50Sn│51Sb│52Te│53 I│54Xe│\n\
  63. 5│ 铷 │ 锶 │ 钇 │ 锆 │ 铌 │ 钼 │ 锝 │ 钌 │ 铑 │ 钯 │ 银 │ 镉 │ 铟 │ 锡 │ 锑 │ 碲 │ 碘 │ 氙 │\n\
  64.   │85.5│87.6│88.9│91.2│92.9│95.9│[99]│101.│103.│106.│108.│112.│115.│119.│122.│128.│127.│131.│\n\
  65.   ├──┼──┼──┼──┼──┼──┼──┼──┼──┼──┼──┼──┼──┼──┼──┼──┼──┼──┤\n\
  66.   │55Cs│56Ba│57- │72Hf│73Ta│74 W│75Re│76Os│77Ir│78Pt│79Au│80Hg│81Tl│82Pb│83Bi│84Po│85At│86Rn│\n\
  67. 6│ 铯 │ 钡 │  71│ 铪 │ 钽 │ 钨 │ 铼 │ 锇 │ 铱 │ 铂 │ 金 │ 汞 │ 铊 │ 铅 │ 铋 │ 钋 │ 砹 │ 氡 │\n\
  68.   │133.│137.│镧系│179.│181.│184.│186.│190.│192.│195.│197.│201.│204.│207.│209.│209*│210*│222*│\n\
  69.   ├──┼──┼──┼──┼──┼──┼──┼──┼──┼──┼──┼──┼──┼──┼──┼──┼──┼──┤\n\
  70.   │87Fr│88Ra│89- │104Rf 105Db 106Sg 107Bh 108Hs 109Mt 110Ds 111Rg 112Cn 113 │114 │115 │116 │117 │118 │\n\
  71. 7│ 钫 │ 镭 │ 103│钅卢│钅杜│钅喜│钅波│钅黑│钅麦│钅达│钅仑│ 鎶 │ Uut│ Fl │ Uup│ Lv │ Uus│ Uuo│\n\
  72.   │223*│226.│锕系│265*│268*│271*│270*│277*│276*│281*│280*│285*│284*│289*│289*│283*│294*│294*│\n\
  73.   └──┴──┴──┴──┴──┴──┴──┴──┴──┴──┴──┴──┴──┴──┴──┴──┴──┴──┘\n\
  74. \n\
  75.   ┌──┬──┬──┬──┬──┬──┬──┬──┬──┬──┬──┬──┬──┬──┬──┐\n\
  76. 镧│57La│58Ce│59Pr│60Nd│61Pm│62Sm│63Eu│64Gd│65Tb│66Dy│67Ho│68Er│69Tm│70Yb│71Lu│\n\
  77. 系│ 镧 │ 铈 │ 镨 │ 钕 │ 钷 │ 钐 │ 铕 │ 钆 │ 铽 │ 镝 │ 钬 │ 铒 │ 铥 │ 镱 │ 镥 │\n\
  78.   │139.│140.│141.│144.│147*│150.│152.│157.│159.│163.│165.│167.│169.│173.│175.│\n\
  79.   ├──┼──┼──┼──┼──┼──┼──┼──┼──┼──┼──┼──┼──┼──┼──┤\n\
  80. 锕│89Ac│90Th│91Pa│92 U│93Np│94Pu│95Am│96Cm│97Bk│98Cf│99Es│100Fm 101Md 102No 103Lr \n\
  81. 系│ 锕 │ 钍 │ 镤 │ 铀 │ 镎 │ 钚 │ 镅 │ 锔 │ 锫 │ 锎 │ 锿 │ 镄 │ 钔 │ 锘 │ 铹 │\n\
  82.   │227.│232.│231.│238.│237.│244*│243*│247*│247*│251*│252*│257*│258*│259*│260*│\n\
  83.   └──┴──┴──┴──┴──┴──┴──┴──┴──┴──┴──┴──┴──┴──┴──┘\n"
  84. ////////////////////////////////////////////////////////////////////
  85. #define SUPER_PRIMES "     The larger number of prime numbers is as follows:\n\
  86.   ┌─────────┬──────────────────┐\n\
  87.   │     range N      │   Number of prime numbers PI(N)    │\n\
  88.   ├─────────┼──────────────────┤\n\
  89.   │     10**1        │                                  4 │\n\
  90.   │     10**2        │                                 25 │\n\
  91.   │     10**3        │                                168 │\n\
  92.   │     10**4        │                              1,229 │\n\
  93.   │     10**5        │                              9,592 │\n\
  94.   │     10**6        │                             78,498 │\n\
  95.   │     10**7        │                            664,579 │\n\
  96.   │     10**8        │                          5,761,455 │\n\
  97.   │     10**9        │                         50,847,534 │\n\
  98.   │     10**10       │                        455,052,511 │\n\
  99.   │     10**11       │                      4,118,054,813 │\n\
  100.   │     10**12       │                     37,607,912,018 │\n\
  101.   │     10**13       │                    346,065,536,839 │\n\
  102.   │     10**14       │                  3,204,941,750,802 │\n\
  103.   │     10**15       │                 29,844,570,422,669 │\n\
  104.   │     10**16       │                279,238,341,033,925 │\n\
  105.   │     10**17       │              2,623,557,157,654,233 │\n\
  106.   │     10**18       │             24,739,954,287,740,860 │\n\
  107.   │     10**19       │            234,057,667,276,344,607 │\n\
  108.   │     10**20       │          2,220,819,602,560,918,840 │\n\
  109.   │     10**21       │         21,127,269,486,018,731,928 │\n\
  110.   │     10**22       │        201,467,286,689,315,906,290 │\n\
  111.   │     10**23       │      1,925,320,391,606,803,968,923 │\n\
  112.   │     10**24       │     18,435,599,767,349,200,867,866 │\n\
  113.   │     10**25       │    176,846,309,399,143,769,411,680 │\n\
  114.   │     10**26       │  1,699,246,750,872,437,141,327,603 │\n\
  115.   └─────────┴──────────────────┘"
  116. ////////////////////////////////////////////////////////////////////

  118. /***************定义全局变量*************/
  119. //定义开关词目
  120. #define            SENSITIVE_NUM  20
  121. typedef  enum{                          GET_FIB ,   GET_FACT ,   GET_PI ,   GET_E ,   GET_COMBINATION ,   GET_PERMUTATION ,   GET_DATEDIFFERENCE ,   GET_HASH ,   GET_RESULT ,   GET_STATISTICS ,   GET_SOLUTION ,   GET_ELEMENTS ,   GET_MYTAX ,   GET_INTEREST ,   GET_ANNUITY ,   PRI_LIST ,   PRI_COUNT ,   PRI_INDEX ,   PRI_CHECK ,   PRI_SUM }senum;
  122. static const char* SENSITIVE_WORDS[]={"/GET-FIB", "/GET-FACT", "/GET-PI", "/GET-E", "/GET-COMBINATION", "/GET-PERMUTATION", "/GET-DATEDIFFERENCE", "/GET-HASH", "/GET-RESULT", "/GET-STATISTICS", "/GET-SOLUTION", "/GET-ELEMENTS", "/GET-MYTAX", "/GET-INTEREST", "/GET-ANNUITY", "/PRI-LIST", "/PRI-COUNT", "/PRI-INDEX", "/PRI-CHECK", "/PRI-SUM"};
  123. static const char  SENSITIVE_ARGVS[]={         1,           1,         1,        1,                  2,                  2,                     2,           1,             1,                 1,               1,               0,            2,               3,              3,           1,            1,            1,            1,          1};

  125. /***************库函数声明***************/
  126. //字符串哈希
  127. void GetStringHash(const char*);
  128. //解高次方程
  129. void Univariate_Equation(int, char**);
  130. //逆波兰核心
  131. double RevPolishCore(const char*);
  132. //阶乘函数
  133. void SmallFact(int);
  134. //π函数
  135. int SmallPI(int, long long, long long, long long, long long);
  136. //e函数
  137. void SmallE(int);
  138. //个税计算
  139. double SmallTax(int, int);
  140. //计算日期差
  141. void SmallDateminus(int, int);

  143. /***************功能函数群***************/
  144. //辨关键词
  145. int itifyWORDS(const char* strARGV)
  146. {
  147.         int SN;
  148.         for(SN=0; SN<SENSITIVE_NUM; SN++){
  149.                 char *op=(char*)strARGV, *kp=(char*)SENSITIVE_WORDS[SN];
  150.                 while(*kp!='\0'){
  151.                         if( (('a'<= *op && *op<='z')?*op-32:*op) != (('a'<= *kp && *kp<='z')?*kp-32:*kp) ){break;}
  152.                         op++;kp++;
  153.                 }
  154.                  if(
  155.                         ((*kp=='\0') && (*op==' '||*op=='\t'||*op=='\r'||*op=='\n'||*op=='\0')) ||
  156.                         (int)kp-(int)SENSITIVE_WORDS[SN] >6
  157.                 ){
  158.                         return SN;
  159.                 }
  160.         }                               
  161.         return -1;
  162. }
  163. //计算组合数C(n,m)
  164. long long Combination(long long  n,long long  m)
  165. {
  166.         if(m==1){
  167.                 return n;
  168.         }
  169.         else{
  170.                 if(m>n/2){m=n-m;}
  171.                 long long  sum=1, *r=(long long*)calloc(1000, sizeof(long long)), i, j, count=1;
  172.                 for(j=n;j>n-m;j--){
  173.                         sum*=j,i=1;
  174.                         for(i=2;i<=m;i++){
  175.                                 if(r[i]==0 && sum%i==0){
  176.                                         sum/=i, r[i]=1, count+=1;
  177.                                 }
  178.                         }
  179.                 }
  180.                 if(count==m && sum>0){free(r);return sum;}
  181.         }
  182.         return -1;
  183. }
  184. //计算排列数P(n,m)
  185. long long Permutation(long long n,long long m)
  186. {
  187.         long long  sum=1, i, j;
  188.         for(i=1,j=n; i<=m; i++,j--){
  189.                 sum*=j;
  190.                 if(sum<0){
  191.                           return -1;
  192.                         break;
  193.                 }
  194.        
  195.         }
  196.         return sum;
  197. }
  198. //斐波那契函数
  199. void Fibonacci(int n, int k)
  200. {
  201.         int i,j,tmp=0,add;
  202.         int *Fa=(int*)calloc(110000, sizeof(int)), *Fb=(int*)calloc(110000, sizeof(int));

  204.         Fa[1]=Fb[1]=1;
  205.         if(n==0 &&k==0){
  206.                 printf("fib(0) = 0");
  207.         }else if(n<=2 &&k==0){
  208.                 printf("fib(%d) = 1",n);
  209.         }else{
  210.                 int N=ceil(n*log10((1+sqrt(5))/2)-log10(5)/2);
  211.                 for(i=3;i<=n;i++){
  212.                         add=0;
  213.                         for(j=1;j<=N/9+1;j++){
  214.                                 tmp=Fa[j]+Fb[j]+add, Fa[j]=Fb[j], add=tmp/1000000000, Fb[j]=tmp%1000000000;
  215.                         }
  216.                 }
  217.                 if(n<45){
  218.                         printf("fib(%d) = %d", n, Fb[1]);
  219.                 }else if(n>=45){
  220.                         printf("fib(%d) = %d", n, Fb[N/9+1]);
  221.                         for(j=N/9;j>1;j--){
  222.                                 printf("%09d",Fb[j]);
  223.                         }
  224.                         printf("%09d",Fb[1]);
  225.                 }
  226.         }
  227.         free(Fa);
  228.         free(Fb);
  229. }
  230. //统计函数
  231. void Statistics(int sargc, char** sargv)
  232. {
  233.         double s=0,q=0;
  234.         int i;
  235.         for(i=2;i<sargc;i++){
  236.                 s+=atof(sargv[i]);
  237.         }
  238.         double e=s/(sargc-2);
  239.         for(i=2;i<sargc;i++){
  240.                 q+=(atof(sargv[i])-e)*(atof(sargv[i])-e);
  241.         }
  242.         fprintf(stdout,
  243.                 "[Statistics]\n"
  244.                 "  Sample expectation  E(x)=%16.12lf\n"
  245.                 "  Sample variance     D(x)=%16.12lf\n"
  246.                 "  Free deviation      Q(x)=%16.12lf\n"
  247.                 "  Standard deviation  V(x)=%16.12lf\n"
  248.                 "  And values          S(x)=%16.12lf\n"
  249.                 ,
  250.                 e, q/(sargc-2), sqrt(q/(sargc-3)), sqrt(q/(sargc-2)), s
  251.         );
  252. }



  256. /*************MAIN主函数入口*************/
  257. int main(int argc, char** argv)
  258. {
  259.         //参数过少则抛出使用说明
  260.         if( argc==1 || (argc>1 && argv[1][0]=='/' && (argv[1][1]=='?'||argv[1][1]=='h'||argv[1][1]=='H')) ){
  261.                 fputs(HELP_INFORMATION, stdout);
  262.                 exit(0);
  263.         }

  265.         //解析开关命令
  266.         //setlocale(LC_ALL, "chs");
  267.         int SN=itifyWORDS(argv[1]);

  269.         if(!(0<=SN &&SN<=19)){
  270.                 //无法识别的开关
  271.                 fprintf(stderr, "Unrecognized switch '%s'\n", argv[1]);
  272.                 exit(1);
  273.         }

  275.         if(argc< 1+1+SENSITIVE_ARGVS[SN]){
  276.                 //缺少必要参数
  277.                 fprintf(stderr, "The switch '%s' needs %d parameters\n", argv[1], SENSITIVE_ARGVS[SN]);
  278.                 exit(1);
  279.         }

  281.         //执行开关命令
  282.         switch((senum)SN)
  283.         {
  284.         //GET系列开关
  285.         case GET_FIB:
  286.                 Fibonacci(atoi(argv[2]), 0);
  287.                 break;
  288.         case GET_FACT:
  289.                 SmallFact(atoi(argv[2]));
  290.                 break;
  291.         case GET_PI:
  292.                 {
  293.                         int n=((int)(atoll(argv[2])/4+1)*4)*3.5;
  294.                         if(n>126000){
  295.                                 printf("超出设定计算范围");
  296.                         }else{
  297.                                 SmallPI(n,0,0,0,0);
  298.                         }
  299.                 }
  300.                 break;
  301.         case GET_E:
  302.                 SmallE(atoi(argv[2]));
  303.                 break;
  304.         case GET_COMBINATION:
  305.                 fprintf(stdout, "%lld\n", Combination(atoll(argv[2]),  atoll(argv[3])) );
  306.                 break;
  307.         case GET_PERMUTATION:
  308.                 fprintf(stdout, "%lld\n", Permutation(atoll(argv[2]),  atoll(argv[3])) );
  309.                 break;
  310.         case GET_DATEDIFFERENCE:
  311.                 SmallDateminus(atoi(argv[2]),  atoi(argv[3]));
  312.                 break;
  313.         case GET_HASH:
  314.                 GetStringHash((const char*)argv[2]);
  315.                 break;
  316.         case GET_RESULT:
  317.                 fprintf(stdout, "%.12lf\n", RevPolishCore(argv[2]));
  318.                 break;
  319.         case GET_STATISTICS:
  320.                 Statistics(argc, argv);
  321.                 break;
  322.         case GET_SOLUTION:
  323.                 Univariate_Equation(argc, argv);
  324.                 break;
  325.         case GET_ELEMENTS:
  326.                 system("mode con cols=120 lines=30");
  327.                 fputs(CHEMICAL_ELEMENTS, stdout);
  328.                 break;
  329.         case GET_MYTAX:
  330.                 SmallTax(atoi(argv[2]), atoi(argv[3]));
  331.                 break;
  332.         case GET_INTEREST:
  333.                 fprintf(stdout, "%.12lf\n", atof(argv[2])*pow(atof(argv[3])+1.0,atof(argv[4])) );
  334.                 break;
  335.         case GET_ANNUITY:
  336.                 fprintf(stdout, "%.12lf\n", atof(argv[2])*(pow(atof(argv[3])+1,atoi(argv[4]))-1)/atof(argv[3]));
  337.                 break;

  339.         //素数系列开关
  340.         case PRI_LIST:
  341.                 if(argc==3){
  342.                         long long n =atoll(argv[2]);
  343.                         if(n<500){
  344.                                 preprime(n);
  345.                         }else if(n<=30000000){
  346.                                 prime((int)(n));
  347.                         }else if(n<999900000000){
  348.                                 class_prime(1000000);
  349.                                 prif(n);
  350.                         }else{
  351.                                 fputs(SUPER_PRIMES, stdout);       
  352.                         }
  353.                 }else if(argc==4){  
  354.                         long long N =atoll(argv[3]);          
  355.                         if(N<=30000000){
  356.                                 listprimes(atoll(argv[2]),atoll(argv[3]));
  357.                         }else if(N<999900000000){
  358.                                 class_prime(1000000);
  359.                                 long long M=atoll(argv[2]);
  360.                                 M=M-(M&1),N=N+(N&1);
  361.                                 if(N-M<=10000000 && N-M>0){
  362.                                         int i, t=FFJ(M, N), T=0, LP;
  363.                                         bool* bm=get_boolp();
  364.                                         printf("%s 到 %s 的素数切片，共计%d个素数\n", argv[2], argv[3], t);
  365.                                         printf("┌────────────────────────────────┐\n");
  366.                                         for(i=1; i<=10000000; i+=2){
  367.                                                 if(!bm[i]){
  368.                                                         T++;
  369.                                                         printf("%13lld", i+M);
  370.                                                         if(T%5==0){printf("\n");}
  371.                                                         if(T==t){LP=i;break;}
  372.                                                 }
  373.                                         }
  374.                                         printf( "\n└────────────────────────────────┘"
  375.                                                 "\n该范围内最大的素数为：%lld"
  376.                                                 ,M+LP
  377.                                         );               
  378.                                 }

  380.                         }else{
  381.                                 printf("超出切片范围");
  382.                         }
  383.                 }
  384.                 break;
  385.         case PRI_COUNT:
  386.                 {
  387.                         long long n=atoll(argv[2]);
  388.                         if(n<=50000000000){
  389.                                 countprimes(n);
  390.                         }else{
  391.                                 fputs(SUPER_PRIMES, stdout);
  392.                         }
  393.                 }
  394.                 break;
  395.         case PRI_INDEX:
  396.                 if(argc==3){
  397.                         indexprime(atoll(argv[2]));
  398.                        
  399.                 }else{
  400.                         index2prime(atoll(argv[2]),atoll(argv[3]));
  401.                 }
  402.                 break;
  403.         case PRI_CHECK:
  404.             if(isPrime(atoll(argv[2]),5)){
  405.             fprintf(stdout, "%s是素数",  argv[2]);
  406.         }
  407.         else{
  408.             fprintf(stdout, "%s不是素数",argv[2]);
  409.         }
  410.                 break;
  411.         case PRI_SUM:
  412.                 {
  413.                         long long n =atoll(argv[2]);
  414.                         if(n<=4200000000){
  415.                                 fprintf(stdout, "%lld", primesum(n));
  416.                         }
  417.                 }
  418.                 break;
  419.         }

  421.         return 0;
  422. }
复制代码

评分

参与人数 2技术 +2 收起 理由
523066680 + 1 回来点赞！
老刘1号 + 1 666

查看全部评分

回复

举报

codegay
发表于 2017-2-24 09:49:46 | 显示全部楼层
i new i new new i
回复

举报

happy886rr
 楼主| 发表于 2017-2-24 10:17:16 | 显示全部楼层
回复 2# codegay
，不错，就叫这个名字吧。inew
回复

举报

niudd
发表于 2017-2-28 11:38:41 | 显示全部楼层
很牛X的代码，学习了
回复

举报

niudd
发表于 2017-2-28 11:38:57 | 显示全部楼层
很牛X的代码，学习了
回复

举报

老刘1号
发表于 2017-4-19 22:46:57 | 显示全部楼层
突然看到了周期表……吃鲸……
回复

举报

wu198726
发表于 2025-6-30 11:03:07 | 显示全部楼层
不错，，学习了，，
回复

举报

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

本版积分规则

Archiver|手机版|小黑屋|批处理之家 ( 渝ICP备10000708号 )

GMT+8, 2026-3-16 22:03 , Processed in 0.027438 second(s), 11 queries , File On.

Powered by Discuz! X3.5

© 2001-2026 Discuz! Team.

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