SELECT  1
USE D:\连续整数方程\相同记录数1至4.DBF ALIAS 相同记录1至4
SELECT  2
USE D:\连续整数方程\相同记录数5加6.DBF ALIAS 相同记录56
SELECT  3
USE D:\连续整数方程\自然数二元运算结果表.DBF ALIAS 自然数二元运算
kssj=SECONDS()
 

          SELECT  2
          GO 1
          For j=1 to 4300
          @ 15,22 say j
          A=自然数
          B=相同记录数
              SELECT 3
             INSERT INTO 自然数二元运算 (自然数,方法数) SELECT 自然数+A ,方法数*B FROM 相同记录1至4
          SELECT 2
          SKIP
          ENDFOR
    =MESSAGEBOX("运行时间："+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
制作方阵对方阵二元运算结果表

SELECT  1
USE D:\连续整数方程\三生素数逆1.DBF ALIAS 三生逆1
SELECT  2
USE D:\连续整数方程\三生素数逆2.DBF ALIAS 三生逆2
SELECT  3
USE D:\连续整数方程\三生素数二元运算表逆.DBF ALIAS 三生二元表
kssj=SECONDS()
 
          SELECT  2
          GO 1
          For j=1 to 1009
          @ 15,22 say j
          B=三中42
         
             SELECT 3
             INSERT INTO 三生二元表 (合成数) SELECT 三中42+B FROM 三生逆1
          SELECT  2
          SKIP
          ENDFOR
    =MESSAGEBOX("运行时间："+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
制作三生素数连续追加方式的二元运算（本程序实际上是完成一个二维表的加法运算并储存）

SELECT  1
USE D:\连续整数方程\相同记录数1至4.DBF ALIAS 相同记录1至4
SELECT  2
USE D:\连续整数方程\相同记录数5加6.DBF ALIAS 相同记录56
SELECT  3
USE D:\连续整数方程\自然数二元运算结果表.DBF ALIAS 自然数二元运算
kssj=SECONDS()
 

          SELECT  2
          GO 1
          For j=1 to 4300
          @ 15,22 say j
          A=自然数
          B=相同记录数
              SELECT 3
             INSERT INTO 自然数二元运算 (自然数,方法数) SELECT 自然数+A ,方法数*B FROM 相同记录1至4
          SELECT 2
          SKIP
          ENDFOR
    =MESSAGEBOX("运行时间："+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
连续整数方程\制作方阵对方阵二元运算结果表
两个方阵的二元同时运算，一个方阵计算类名，另一个方阵计算合成方法数，形成两个字段的vfp表，类名与合成方法数关联成一个记录条，然后继续操作，汇总同类名的合成方法数。

SELECT 1
USE D:\素数式裂解\素数式8.DBF ALIAS 素数式8
SELECT 2
USE D:\素数式裂解\素数式8横.DBF ALIAS 素数式8横
kssj=SECONDS()
    SELECT 1
    ss1=0
    ss8=42
    For i=1 to 38760
      ss2=素式8
      SKIP
      ss3=素式8
      SKIP
      ss4=素式8
      SKIP
      ss5=素式8
       SKIP
      ss6=素式8
      skip
      ss7=素式8
       SELECT  2  &&打开盛放素数式的表
       APPEND BLANK     &&增加一条空记录
       REPLACE 素1 WITH ss1 &&将B值付给素数式
       REPLACE 素2 WITH ss2 &&将B值付给素数式
       REPLACE 素3 WITH ss3 &&将B值付给素数式
       REPLACE 素4 WITH ss4 &&将B值付给素数式
       REPLACE 素5 WITH ss5 &&将B值付给素数式
       REPLACE 素6 WITH ss6 &&将B值付给素数式
       REPLACE 素7 WITH ss7 &&将B值付给素数式
       REPLACE 素8 WITH ss8 &&将B值付给素数式
       SELECT 1
       skip
   endfor
    =MESSAGEBOX("运行时间："+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
先用程序制作各种情况的素数式（对开头的“0”及末尾的“42”不做处理），只产生2至40的一个抽取6个的组合数。（制作一个仅含一个字段的vfp表）。
     上述完成后，用本程序制作一个含8个字段的vfp表，记录条是从0到42的一个组合（每条记录的第一个字段的值为0，末字段（即第八个字段的）值是42，说白了就是在2-40之间抽取6个数的组合，再加个0和42构成含8个值的记录条，只是偶数，奇数不参与）。

SELECT 1
USE D:\三角递增法\四生素数表中.DBF ALIAS 四生素数中
SELECT 2
USE D:\三角递增法\各段四生素数个数细分.DBF ALIAS 四生个数细分
kssj=SECONDS()
SELECT 1
     For j=1 to 500
      @ 3,6 say j
      fw=j*10^7
            SELECT 1
            GO 1
            COUNT ALL FOR fw>四中 TO tj
              SELECT  2  &&打开盛放素数式的表
              APPEND BLANK     &&增加一条空记录
              REPLACE 范围 WITH fw
              REPLACE 总个数 WITH tj
     ENDFOR
   
    =MESSAGEBOX("运行时间："+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
这是把四生素数表分成500个段落（每段落1000万，总区间跨度50亿），为后续工作打基础，我们知道任何软件都有尺寸（容量）限制，所以好多计算需要分步完成。

在这个问题上laowan001版主给予很大帮助。在我开始学习vfp编程时，得到那宝吉先生的鼎力帮助，从一个菜鸟到独立完成作业，没有那宝吉先生的帮助是无法自己完成的。

算法来源于，问题处理方式，和最后解决问题，前提是：你必须对问题有深刻的理解。

CLOSE DATA
USE d:\连续整数方程\三生素数三元运算结果表逆.dbf
sele 合成数,sum(方法数) as 方法数 from 三生素数三元运算结果表逆 group by 合成数  into dbf d:\连续整数方程\三生三元运算统计表逆
USE IN 三生三元运算统计表逆
时间久了，自己把自己写出来的东西都给忘个差不多了。

编程有时也需要迂回战术，如果不能直达，不妨先解决相关，近似问题，然后在反过头来从新考虑，是否可以在那种基础上发展下去。

单位元矩阵算法：
模5矩阵    0.2    0.4    0.6    0.8
0.2    0.4    0.6    0.8    1
0.4    0.6    0.8    1    1.2
0.6    0.8    1    1.2    1.4
0.8    1    1.2    1.4    1.6
5的剩余类    统计2
0.4    1
0.6    2
0.8    3
1    4
1.2    3
1.4    2
1.6    1
合计    16
周期    0    1    2    3    4    5    6    7    8    9
0    0    1    2    3    4    5    6    7    8    9
1    1    2    3    4    5    6    7    8    9    10
2    2    3    4    5    6    7    8    9    10    11
3    3    4    5    6    7    8    9    10    11    12
4    4    5    6    7    8    9    10    11    12    13
5    5    6    7    8    9    10    11    12    13    14
6    6    7    8    9    10    11    12    13    14    15
7    7    8    9    10    11    12    13    14    15    16
8    8    9    10    11    12    13    14    15    16    17
9    9    10    11    12    13    14    15    16    17    18
周期    统计
0    1
1    2
2    3
3    4
4    5
5    6
6    7
7    8
8    9
9    10
10    9
11    8
12    7
13    6
14    5
15    4
16    3
17    2
18    1
合计    100
这种算法，是耦合的实际应用，为什么，把一个算法要发在这里，一个非常重要的客观原因就是，算法来源于，解决实际问题，我的合成方法论是解决线性不定方程的满足条件的正整数解组数。