注册 登录
编程论坛 VFP论坛

通过单位矩阵与周期矩阵耦合统计方程解组数

独木星空 发布于 2022-10-29 15:24, 4911 次点击
我构造了一种算法,可以解决限制条件下,线性不定方程x+y+z+u+v+m=N满足条件的解组数,条件是未知数不能取某一个大于1的整数的倍数,如未知数不能取2t(t为正整数),这是一种限制;有或者不能取3t的数;不能取4t的数;5t的数;....,等等。这里的Pt(t为任意的正整数,P为大于1的正整数)中,P为单位矩阵的单位(相当于周期值),t为周期,我们限定t=10,这样对于任一个P来说,我们就不从新计算周期矩阵了,而只计算单位矩阵。
周期矩阵计算:(第一步)
周期10    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

周期数    统计2
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

第二步
周期数    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
10    10    11    12    13    14    15    16    17    18    19
11    11    12    13    14    15    16    17    18    19    20
12    12    13    14    15    16    17    18    19    20    21
13    13    14    15    16    17    18    19    20    21    22
14    14    15    16    17    18    19    20    21    22    23
15    15    16    17    18    19    20    21    22    23    24
16    16    17    18    19    20    21    22    23    24    25
17    17    18    19    20    21    22    23    24    25    26
18    18    19    20    21    22    23    24    25    26    27

统计2    1    2    3    4    5    6    7    8    9    10
1    1    2    3    4    5    6    7    8    9    10
2    2    4    6    8    10    12    14    16    18    20
3    3    6    9    12    15    18    21    24    27    30
4    4    8    12    16    20    24    28    32    36    40
5    5    10    15    20    25    30    35    40    45    50
6    6    12    18    24    30    36    42    48    54    60
7    7    14    21    28    35    42    49    56    63    70
8    8    16    24    32    40    48    56    64    72    80
9    9    18    27    36    45    54    63    72    81    90
10    10    20    30    40    50    60    70    80    90    100
9    9    18    27    36    45    54    63    72    81    90
8    8    16    24    32    40    48    56    64    72    80
7    7    14    21    28    35    42    49    56    63    70
6    6    12    18    24    30    36    42    48    54    60
5    5    10    15    20    25    30    35    40    45    50
4    4    8    12    16    20    24    28    32    36    40
3    3    6    9    12    15    18    21    24    27    30
2    2    4    6    8    10    12    14    16    18    20
1    1    2    3    4    5    6    7    8    9    10

周期数    统计4
0    1
1    4
2    10
3    20
4    35
5    56
6    84
7    120
8    165
9    220
10    282
11    348
12    415
13    480
14    540
15    592
16    633
17    660
18    670
19    660
20    633
21    592
22    540
23    480
24    415
25    348
26    282
27    220
28    165
29    120
30    84
31    56
32    35
33    20
34    10
35    4
36    1
合计    10000
这第二步仅列出了10列数据,后9列数据未列出,补全即可。周期数二元运算符号“+”;统计二元运算符号“*”;统计是把周期数一样的,对应位置的统计数加在一起:SUMIF(B$24:T$42,AQ24,W$24:AO$42),Excel中的函数。
57 回复
#2
sam_jiang2022-10-29 15:25
干嘛用的?
#3
独木星空2022-10-29 15:38
第三步
周期数    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
10    10    11    12    13    14    15    16    17    18    19
11    11    12    13    14    15    16    17    18    19    20
12    12    13    14    15    16    17    18    19    20    21
13    13    14    15    16    17    18    19    20    21    22
14    14    15    16    17    18    19    20    21    22    23
15    15    16    17    18    19    20    21    22    23    24
16    16    17    18    19    20    21    22    23    24    25
17    17    18    19    20    21    22    23    24    25    26
18    18    19    20    21    22    23    24    25    26    27
19    19    20    21    22    23    24    25    26    27    28
20    20    21    22    23    24    25    26    27    28    29
21    21    22    23    24    25    26    27    28    29    30
22    22    23    24    25    26    27    28    29    30    31
23    23    24    25    26    27    28    29    30    31    32
24    24    25    26    27    28    29    30    31    32    33
25    25    26    27    28    29    30    31    32    33    34
26    26    27    28    29    30    31    32    33    34    35
27    27    28    29    30    31    32    33    34    35    36
28    28    29    30    31    32    33    34    35    36    37
29    29    30    31    32    33    34    35    36    37    38
30    30    31    32    33    34    35    36    37    38    39
31    31    32    33    34    35    36    37    38    39    40
32    32    33    34    35    36    37    38    39    40    41
33    33    34    35    36    37    38    39    40    41    42
34    34    35    36    37    38    39    40    41    42    43
35    35    36    37    38    39    40    41    42    43    44
36    36    37    38    39    40    41    42    43    44    45

统计4/2    1    2    3    4    5    6    7    8    9    10
1    1    2    3    4    5    6    7    8    9    10
4    4    8    12    16    20    24    28    32    36    40
10    10    20    30    40    50    60    70    80    90    100
20    20    40    60    80    100    120    140    160    180    200
35    35    70    105    140    175    210    245    280    315    350
56    56    112    168    224    280    336    392    448    504    560
84    84    168    252    336    420    504    588    672    756    840
120    120    240    360    480    600    720    840    960    1080    1200
165    165    330    495    660    825    990    1155    1320    1485    1650
220    220    440    660    880    1100    1320    1540    1760    1980    2200
282    282    564    846    1128    1410    1692    1974    2256    2538    2820
348    348    696    1044    1392    1740    2088    2436    2784    3132    3480
415    415    830    1245    1660    2075    2490    2905    3320    3735    4150
480    480    960    1440    1920    2400    2880    3360    3840    4320    4800
540    540    1080    1620    2160    2700    3240    3780    4320    4860    5400
592    592    1184    1776    2368    2960    3552    4144    4736    5328    5920
633    633    1266    1899    2532    3165    3798    4431    5064    5697    6330
660    660    1320    1980    2640    3300    3960    4620    5280    5940    6600
670    670    1340    2010    2680    3350    4020    4690    5360    6030    6700
660    660    1320    1980    2640    3300    3960    4620    5280    5940    6600
633    633    1266    1899    2532    3165    3798    4431    5064    5697    6330
592    592    1184    1776    2368    2960    3552    4144    4736    5328    5920
540    540    1080    1620    2160    2700    3240    3780    4320    4860    5400
480    480    960    1440    1920    2400    2880    3360    3840    4320    4800
415    415    830    1245    1660    2075    2490    2905    3320    3735    4150
348    348    696    1044    1392    1740    2088    2436    2784    3132    3480
282    282    564    846    1128    1410    1692    1974    2256    2538    2820
220    220    440    660    880    1100    1320    1540    1760    1980    2200
165    165    330    495    660    825    990    1155    1320    1485    1650
120    120    240    360    480    600    720    840    960    1080    1200
84    84    168    252    336    420    504    588    672    756    840
56    56    112    168    224    280    336    392    448    504    560
35    35    70    105    140    175    210    245    280    315    350
20    20    40    60    80    100    120    140    160    180    200
10    10    20    30    40    50    60    70    80    90    100
4    4    8    12    16    20    24    28    32    36    40
1    1    2    3    4    5    6    7    8    9    10
同样只列出了10列数据,后9列数据未列出,补全即可。
周期数    统计6
0    1
1    6
2    21
3    56
4    126
5    252
6    462
7    792
8    1287
9    2002
10    2997
11    4332
12    6062
13    8232
14    10872
15    13992
16    17577
17    21582
18    25927
19    30492
20    35127
21    39662
22    43917
23    47712
24    50877
25    53262
26    54747
27    55252
28    54747
29    53262
30    50877
31    47712
32    43917
33    39662
34    35127
35    30492
36    25927
37    21582
38    17577
39    13992
40    10872
41    8232
42    6062
43    4332
44    2997
45    2002
46    1287
47    792
48    462
49    252
50    126
51    56
52    21
53    6
54    1
合计    1000000
0周期代表第一个周期,10个周期用0至9表示,最大值是9,所以最大周期数为9*6=54,所以六个未知数的和最大值是54,最小值是0,它们的和全部分布到0-54上,所以,统计6的合计是10^6=100万。
这对于10个周期已经处理完毕。   周期矩阵算法步骤
#4
独木星空2022-10-29 15:41
模7矩阵    1    2    3    4    5    6
1    2    3    4    5    6    7
2    3    4    5    6    7    8
3    4    5    6    7    8    9
4    5    6    7    8    9    10
5    6    7    8    9    10    11
6    7    8    9    10    11    12

7的剩余类    统计2
2    1
3    2
4    3
5    4
6    5
7    6
8    5
9    4
10    3
11    2
12    1
合计    36
周期值是7,未知数不能取7的倍数,单位矩阵第一步二元运算,运算符“+”,统计合成值相同的个数。
#5
独木星空2022-10-29 15:44
第二步
7的剩余类    2    3    4    5    6    7    8    9    10    11    12
2    4    5    6    7    8    9    10    11    12    13    14
3    5    6    7    8    9    10    11    12    13    14    15
4    6    7    8    9    10    11    12    13    14    15    16
5    7    8    9    10    11    12    13    14    15    16    17
6    8    9    10    11    12    13    14    15    16    17    18
7    9    10    11    12    13    14    15    16    17    18    19
8    10    11    12    13    14    15    16    17    18    19    20
9    11    12    13    14    15    16    17    18    19    20    21
10    12    13    14    15    16    17    18    19    20    21    22
11    13    14    15    16    17    18    19    20    21    22    23
12    14    15    16    17    18    19    20    21    22    23    24

统计2    1    2    3    4    5    6    5    4    3    2    1
1    1    2    3    4    5    6    5    4    3    2    1
2    2    4    6    8    10    12    10    8    6    4    2
3    3    6    9    12    15    18    15    12    9    6    3
4    4    8    12    16    20    24    20    16    12    8    4
5    5    10    15    20    25    30    25    20    15    10    5
6    6    12    18    24    30    36    30    24    18    12    6
5    5    10    15    20    25    30    25    20    15    10    5
4    4    8    12    16    20    24    20    16    12    8    4
3    3    6    9    12    15    18    15    12    9    6    3
2    2    4    6    8    10    12    10    8    6    4    2
1    1    2    3    4    5    6    5    4    3    2    1

7的剩余类    统计4
4    1
5    4
6    10
7    20
8    35
9    56
10    80
11    104
12    125
13    140
14    146
15    140
16    125
17    104
18    80
19    56
20    35
21    20
22    10
23    4
24    1
合计    1296
剩余类数据区域用加法,统计区域数据用乘法,首列+首行,首列*首行,统计函数:SUMIF(B$16:L$26,AA16,O$16:Y$26)
#6
独木星空2022-10-29 15:46
第三步
7的剩余类    2    3    4    5    6    7    8
4    6    7    8    9    10    11    12
5    7    8    9    10    11    12    13
6    8    9    10    11    12    13    14
7    9    10    11    12    13    14    15
8    10    11    12    13    14    15    16
9    11    12    13    14    15    16    17
10    12    13    14    15    16    17    18
11    13    14    15    16    17    18    19
12    14    15    16    17    18    19    20
13    15    16    17    18    19    20    21
14    16    17    18    19    20    21    22
15    17    18    19    20    21    22    23
16    18    19    20    21    22    23    24
17    19    20    21    22    23    24    25
18    20    21    22    23    24    25    26
19    21    22    23    24    25    26    27
20    22    23    24    25    26    27    28
21    23    24    25    26    27    28    29
22    24    25    26    27    28    29    30
23    25    26    27    28    29    30    31
24    26    27    28    29    30    31    32

统计4    1    2    3    4    5    6    5
1    1    2    3    4    5    6    5
4    4    8    12    16    20    24    20
10    10    20    30    40    50    60    50
20    20    40    60    80    100    120    100
35    35    70    105    140    175    210    175
56    56    112    168    224    280    336    280
80    80    160    240    320    400    480    400
104    104    208    312    416    520    624    520
125    125    250    375    500    625    750    625
140    140    280    420    560    700    840    700
146    146    292    438    584    730    876    730
140    140    280    420    560    700    840    700
125    125    250    375    500    625    750    625
104    104    208    312    416    520    624    520
80    80    160    240    320    400    480    400
56    56    112    168    224    280    336    280
35    35    70    105    140    175    210    175
20    20    40    60    80    100    120    100
10    10    20    30    40    50    60    50
4    4    8    12    16    20    24    20
1    1    2    3    4    5    6    5

7的剩余类    统计6
6    1
7    6
8    21
9    56
10    126
11    252
12    456
13    756
14    1161
15    1666
16    2247
17    2856
18    3431
19    3906
20    4221
21    4332
22    4221
23    3906
24    3431
25    2856
26    2247
27    1666
28    1161
29    756
30    456
31    252
32    126
33    56
34    21
35    6
36    1
合计    46656
数据区域后4列未列出,补全即可,运算同上楼。
#7
独木星空2022-10-29 16:01
算到上一步,最后就是:单位矩阵与周期矩阵耦合了,即把周期矩阵的周期数*周期值,作为运算元素,而单位矩阵的值直接参与加法运算。
     统计数据是相乘关系,最后做一下统计,10个周期的N就出来了,7*10=70,这样计算出来的数,N=70以前正确,以后的不确保。
周期/分值    6    7    8    9    10    11    12
0    6    7    8    9    10    11    12
7    13    14    15    16    17    18    19
14    20    21    22    23    24    25    26
21    27    28    29    30    31    32    33
28    34    35    36    37    38    39    40
35    41    42    43    44    45    46    47
42    48    49    50    51    52    53    54
49    55    56    57    58    59    60    61
56    62    63    64    65    66    67    68
63    69    70    71    72    73    74    75
70    76    77    78    79    80    81    82
77    83    84    85    86    87    88    89
84    90    91    92    93    94    95    96
91    97    98    99    100    101    102    103
98    104    105    106    107    108    109    110
105    111    112    113    114    115    116    117
112    118    119    120    121    122    123    124
119    125    126    127    128    129    130    131
126    132    133    134    135    136    137    138
133    139    140    141    142    143    144    145
140    146    147    148    149    150    151    152
147    153    154    155    156    157    158    159
154    160    161    162    163    164    165    166
161    167    168    169    170    171    172    173
168    174    175    176    177    178    179    180
175    181    182    183    184    185    186    187
182    188    189    190    191    192    193    194
189    195    196    197    198    199    200    201
196    202    203    204    205    206    207    208
203    209    210    211    212    213    214    215
210    216    217    218    219    220    221    222
217    223    224    225    226    227    228    229
224    230    231    232    233    234    235    236
231    237    238    239    240    241    242    243
238    244    245    246    247    248    249    250
245    251    252    253    254    255    256    257
252    258    259    260    261    262    263    264
259    265    266    267    268    269    270    271
266    272    273    274    275    276    277    278
273    279    280    281    282    283    284    285
280    286    287    288    289    290    291    292
287    293    294    295    296    297    298    299
294    300    301    302    303    304    305    306
301    307    308    309    310    311    312    313
308    314    315    316    317    318    319    320
315    321    322    323    324    325    326    327
322    328    329    330    331    332    333    334
329    335    336    337    338    339    340    341
336    342    343    344    345    346    347    348
343    349    350    351    352    353    354    355
350    356    357    358    359    360    361    362
357    363    364    365    366    367    368    369
364    370    371    372    373    374    375    376
371    377    378    379    380    381    382    383
378    384    385    386    387    388    389    390

统计6    1    6    21    56    126    252    456
1    1    6    21    56    126    252    456
6    6    36    126    336    756    1512    2736
21    21    126    441    1176    2646    5292    9576
56    56    336    1176    3136    7056    14112    25536
126    126    756    2646    7056    15876    31752    57456
252    252    1512    5292    14112    31752    63504    114912
462    462    2772    9702    25872    58212    116424    210672
792    792    4752    16632    44352    99792    199584    361152
1287    1287    7722    27027    72072    162162    324324    586872
2002    2002    12012    42042    112112    252252    504504    912912
2997    2997    17982    62937    167832    377622    755244    1366632
4332    4332    25992    90972    242592    545832    1091664    1975392
6062    6062    36372    127302    339472    763812    1527624    2764272
8232    8232    49392    172872    460992    1037232    2074464    3753792
10872    10872    65232    228312    608832    1369872    2739744    4957632
13992    13992    83952    293832    783552    1762992    3525984    6380352
17577    17577    105462    369117    984312    2214702    4429404    8015112
21582    21582    129492    453222    1208592    2719332    5438664    9841392
25927    25927    155562    544467    1451912    3266802    6533604    11822712
30492    30492    182952    640332    1707552    3841992    7683984    13904352
35127    35127    210762    737667    1967112    4426002    8852004    16017912
39662    39662    237972    832902    2221072    4997412    9994824    18085872
43917    43917    263502    922257    2459352    5533542    11067084    20026152
47712    47712    286272    1001952    2671872    6011712    12023424    21756672
50877    50877    305262    1068417    2849112    6410502    12821004    23199912
53262    53262    319572    1118502    2982672    6711012    13422024    24287472
54747    54747    328482    1149687    3065832    6898122    13796244    24964632
55252    55252    331512    1160292    3094112    6961752    13923504    25194912
54747    54747    328482    1149687    3065832    6898122    13796244    24964632
53262    53262    319572    1118502    2982672    6711012    13422024    24287472
50877    50877    305262    1068417    2849112    6410502    12821004    23199912
47712    47712    286272    1001952    2671872    6011712    12023424    21756672
43917    43917    263502    922257    2459352    5533542    11067084    20026152
39662    39662    237972    832902    2221072    4997412    9994824    18085872
35127    35127    210762    737667    1967112    4426002    8852004    16017912
30492    30492    182952    640332    1707552    3841992    7683984    13904352
25927    25927    155562    544467    1451912    3266802    6533604    11822712
21582    21582    129492    453222    1208592    2719332    5438664    9841392
17577    17577    105462    369117    984312    2214702    4429404    8015112
13992    13992    83952    293832    783552    1762992    3525984    6380352
10872    10872    65232    228312    608832    1369872    2739744    4957632
8232    8232    49392    172872    460992    1037232    2074464    3753792
6062    6062    36372    127302    339472    763812    1527624    2764272
4332    4332    25992    90972    242592    545832    1091664    1975392
2997    2997    17982    62937    167832    377622    755244    1366632
2002    2002    12012    42042    112112    252252    504504    912912
1287    1287    7722    27027    72072    162162    324324    586872
792    792    4752    16632    44352    99792    199584    361152
462    462    2772    9702    25872    58212    116424    210672
252    252    1512    5292    14112    31752    63504    114912
126    126    756    2646    7056    15876    31752    57456
56    56    336    1176    3136    7056    14112    25536
21    21    126    441    1176    2646    5292    9576
6    6    36    126    336    756    1512    2736
1    1    6    21    56    126    252    456

合成值    统计6
6    1
7    6
8    21
9    56
10    126
11    252
12    456
13    762
14    1197
15    1792
16    2583
17    3612
18    4943
19    6642
20    8778
21    11424
22    14658
23    18564
24    23213
25    28734
26    35259
27    42924
28    51870
29    62244
30    74215
31    87870
32    103425
33    121100
34    141120
35    163716
36    189120
37    217650
38    249375
39    284620
40    323715
41    366996
42    414807
43    467466
44    525546
45    589036
46    658392
47    734076
48    816557
49    906318
50    1003737
51    1109808
52    1224342
53    1347948
54    1481242
55    1624848
56    1779420
57    1945328
58    2124234
59    2315628
60    2520294
61    2739024
62    2972619
63    3221946
64    3487239
65    3771180
66    4072740
67    4392900
68    4732650
69    5092990
70    5475057
71    5878740
72    6308225
73    6761700
74    7240365
75    7745430
76    8278110
77    8839872
78    9429910
79    10054500
80    10710615
81    11399570
82    12122685
83    12881280
84    13677132
85    14508180
86    15383535
87    16298450
88    17254245
89    18252240
90    19293750
91    20380872
92    21509550
93    22692670
94    23923215
95    25202400
96    26531435
97    27911520
98    29345127
99    30825270
100    32369730
101    33968640
102    35623000
103    37333800
104    39102015
105    40930602
106    42808515
107    44759700
108    46770850
109    48842640
110    50975730
111    53170760
112    55431342
113    57741060
114    60131400
115    62585040
116    65102220
117    67683160
118    70328055
119    73041402
120    75799955
121    78644160
122    81552150
123    84523620
124    87558240
125    90655650
126    93821517
127    97024200
128    100314495
129    103665540
130    107076375
131    110546010
132    114073420
133    117665772
134    121281420
135    124982780
136    128737665
137    132544350
138    136401075
139    140306040
140    144268272
141    148234540
142    152279925
143    156366750
144    160492415
145    164654280
146    168849675
147    173089902
148    177308775
149    181594890
150    185905355
151    190236900
152    194586225
153    198950010
154    203342502
155    207683610
156    212075905
157    216471780
158    220867515
159    225259370
160    229643595
161    234038022
162    238348095
163    242690280
164    247012935
165    251312110
166    255583845
167    259824180
168    264055092
169    268167480
170    272291055
171    276370990
172    280403325
173    284384100
174    288309365
175    292205682
176    295949865
177    299683650
178    303350025
179    306945240
180    310465555
181    313907250
182    317301762
183    320512850
184    323692515
185    326782680
186    329780025
187    332681250
188    335483085
189    338221962
190    340750185
191    343227780
192    345596805
193    347854590
194    349998495
195    352025920
196    353978232
197    355698420
198    357351905
199    358882110
200    360287235
201    361565520
202    362715255
203    363782502
204    364603555
205    365346330
206    365956815
207    366434300
208    366778125
209    366987690
210    367113342
211    366987690
212    366778125
213    366434300
214    365956815
215    365346330
216    364603555
217    363782502
218    362715255
219    361565520
220    360287235
221    358882110
222    357351905
223    355698420
224    353978232
225    352025920
226    349998495
227    347854590
228    345596805
229    343227780
230    340750185
231    338221962
232    335483085
233    332681250
234    329780025
235    326782680
236    323692515
237    320512850
238    317301762
239    313907250
240    310465555
241    306945240
242    303350025
243    299683650
244    295949865
245    292205682
246    288309365
247    284384100
248    280403325
249    276370990
250    272291055
251    268167480
252    264055092
253    259824180
254    255583845
255    251312110
256    247012935
257    242690280
258    238348095
259    234038022
260    229643595
261    225259370
262    220867515
263    216471780
264    212075905
265    207683610
266    203342502
267    198950010
268    194586225
269    190236900
270    185905355
271    181594890
272    177308775
273    173089902
274    168849675
275    164654280
276    160492415
277    156366750
278    152279925
279    148234540
280    144268272
281    140306040
282    136401075
283    132544350
284    128737665
285    124982780
286    121281420
287    117665772
288    114073420
289    110546010
290    107076375
291    103665540
292    100314495
293    97024200
294    93821517
295    90655650
296    87558240
297    84523620
298    81552150
299    78644160
300    75799955
301    73041402
302    70328055
303    67683160
304    65102220
305    62585040
306    60131400
307    57741060
308    55431342
309    53170760
310    50975730
311    48842640
312    46770850
313    44759700
314    42808515
315    40930602
316    39102015
317    37333800
318    35623000
319    33968640
320    32369730
321    30825270
322    29345127
323    27911520
324    26531435
325    25202400
326    23923215
327    22692670
328    21509550
329    20380872
330    19293750
331    18252240
332    17254245
333    16298450
334    15383535
335    14508180
336    13677132
337    12881280
338    12122685
339    11399570
340    10710615
341    10054500
342    9429910
343    8839872
344    8278110
345    7745430
346    7240365
347    6761700
348    6308225
349    5878740
350    5475057
351    5092990
352    4732650
353    4392900
354    4072740
355    3771180
356    3487239
357    3221946
358    2972619
359    2739024
360    2520294
361    2315628
362    2124234
363    1945328
364    1779420
365    1624848
366    1481242
367    1347948
368    1224342
369    1109808
370    1003737
371    906318
372    816557
373    734076
374    658392
375    589036
376    525546
377    467466
378    414807
379    366996
380    323715
381    284620
382    249375
383    217650
384    189120
385    163716
386    141120
387    121100
388    103425
389    87870
390    74215
391    62244
392    51870
393    42924
394    35259
395    28734
396    23213
397    18564
398    14658
399    11424
400    8778
401    6642
402    4943
403    3612
404    2583
405    1792
406    1197
407    762
408    456
409    252
410    126
411    56
412    21
413    6
414    1
合计    46656000000
数据区域仅列出了7列,后边24列数据未列出。
#8
独木星空2022-10-29 16:06
这是在Excel上的计算过程,我希望把它翻译成VFP编程语言,而且模P可以任意指定那种,P为周期值,未知数不能取周期值P的倍数,求线性不定方程满足条件的解组数,有其他方法也行,我需要6个周期的数据。
#9
独木星空2022-10-29 20:40
实际上就是解决这么一类问题:x+y+z+u+v+m=N,未知数取值不是5的倍数,求线性不定方程满足条件的正整数解组数。用vfp程序解决。
#10
独木星空2022-10-30 07:19
vfp的汇总统计功能是强悍的。就是不知道你是否敢去触碰它。
#11
独木星空2022-10-30 19:38
一般情况下,对于问题的产生比较熟悉的,肯定是有办法把它转化成流程图。
#12
gs25367856782022-10-31 07:46
看过,朋友是在展示你的【独孤九剑】
#13
my23182022-10-31 08:21
解法不知道,不可能写出算法,解法不熟悉,很难写出精炼的算法
#14
独木星空2022-10-31 08:58
laowan001先生也不来试一试,在我的印象里,他还是有办法解决的。
另外,吹水佬版主更不在话下,不知道为什么,就是不出手,一展其锋芒。
#15
laowan0012022-10-31 09:17
回复 14楼 独木星空
你这个问题太专业,一时半会儿搞不懂啊
#16
吹水佬2022-10-31 14:51
以下是引用独木星空在2022-10-31 08:58:41的发言:

laowan001先生也不来试一试,在我的印象里,他还是有办法解决的。
另外,吹水佬版主更不在话下,不知道为什么,就是不出手,一展其锋芒。

老实讲,你那么多东东,我还没看明一点。
你可以先说说:
1、这个问题的概念(是什么名堂)
2、相关的指导理论(理论基础是什么)
人脑有点认知后再看电脑能不能认知


#17
独木星空2022-10-31 15:45
回复 16楼 吹水佬
问题是这样的:x+y+z+u+v+m=N,如果未知数不加限制,则用排列组合中的挡板法可以直接求出,现在的问题是限制未知数取某数的倍数,这个某数是大于等于2即可,我们也可以先研究个具体的,比如不能取5的倍数。
     我们把未知数x(y,z,u,v,m)=Pt+r的形式,与数论中的剩余类(剩余系)有点相似,0<r≤P-1,
这样每个未知数就是两部分组成的,而且t与r各自独立,根据乘法原理x=Pt+r是分步关系,如果t有i种方法,r有j种方法,则对于一个x来说,共有ij种方法。
      x+y+z+u+v+m=N→(x+y)+(z+u+v+m)=N→((x+y)+[(z+u)+(v+m)]=N
我们先打分,充分利用对称原理,把(x+y)分成一个区块,(z+u+v+m)分成一个区块,假说x+y=N1,z+u+v+m=N2,N1+N2=N的情况下,得到N1有j1种方法,得到N2有j2种方法,则共有j1j2种方法,在整个运算过程中,始终是把未知数分区,2个变量和+(2个变量和+2个变量和)=2个变量和+4个变量和=N,(最开始一步是1个变量+1个变量=N2),然后,用2+2=4(数字都是代表的变量个数),2+4=6(同前注解),第一步是1+1=2.
数据区域用“+”,统计区域用“*”,第一步例外,直接统计各种和值的个数,第二步在第一步的基础上进行,把“和值”用在数据区域,首列,顶行都是用它,形成二维数据结构;把统计2的值照样处理,排列成首行,首列二维数据(它的二元运算用“*”,前边的数据区域用“+”),然后统计数据区域和值一样的有多少,即找到统计对应位置上的数求和。
     实际上,这些叙述不一定起多大作用,需要看一下整个帖子实例,每步干什么,最后得到什么数值。
#18
独木星空2022-10-31 16:12
回复 16楼 吹水佬
可以先研究周期矩阵,它们不受限制,先不管单位矩阵,计算10个周期的即可,周期一定从0开始,因为单位矩阵的缘故,也就是X+Y+Z+V+U+M=N中,未知数可以取0-9的数字,看一看它们的分布,肯定之和是0--54,有10^6=100万种方法,落到0--54上,数据分布是对称的。
     当研究整个问题时,只需要把0--54乘限制条件即可(实际上,乘的是周期值,比方不能取5的倍数,那就是0---54*5了)。
      实际上,你研究完周期矩阵,下来研究单位矩阵时,也就好办了,比如不能取5的倍数,那么,未知数取值就是1--4之间,最小值是6,最大值是4*6=24,也就是说,4^6个值一定落到1至24上。
      最后,耦合一步,也简单,用周期倍数值+单位矩阵的和值,作为结果,合成方法数=周期倍数值的统计个数*单位和值的统计个数,说白了,周期矩阵,单位矩阵原来是各自独立,后来它们混合了一下。
#19
独木星空2022-10-31 20:14
解6元一次不定方程的解组数,条件变量不能取5的倍数,给出1至100之间的解组数(即N=1到100)
#20
吹水佬2022-10-31 23:28
以下是引用独木星空在2022-10-31 20:14:27的发言:

解6元一次不定方程的解组数,条件变量不能取5的倍数,给出1至100之间的解组数(即N=1到100)

算出的结果做什么用的?
#21
独木星空2022-11-01 03:54
回复 20楼 吹水佬
只要算出6周的数据,就可以用待定系数法求出多项式公式表达式。一般这类单条件的线性不定方程的解组数公式是一个比未知数个数少一次的多项式形式。它是六元的,所以,解组数表达式是一个最高5次的多项式,它有6个系数要求,所以,最少需要6个周期的数据。
#22
独木星空2022-11-02 07:22
最终还是没有人,一试身手。
#23
吹水佬2022-11-02 07:36
以下是引用独木星空在2022-11-2 07:22:32的发言:

最终还是没有人,一试身手。

你不是算出来了吗? 那一大堆的数是什么东东。
电算有限,手算无限。
#24
xinjie2022-11-02 11:23
先搞定线性代数再来吧。。。。哎
#25
独木星空2022-11-03 14:51
回复 24楼 xinjie
这些表与它的关系,第一步,相当于x+y=N
                  第二步,相当于(x+y)+(z+u)=N
                  第三步,相当于[(x+y)+(z+u)]+(v+m)=N
     这三步,周期矩阵,与单位矩阵,操作步骤一样,只是起点,与参与的元素不同。
     
     单位矩阵,用1到P-1之间的几个元素,不能用0,也不能用P(限制条件值,比如,未知数,或者说变量,不能取7的倍数,则P=7,所用值就是1,2,3,4,5,6了)

      但是周期矩阵是从0开始,到某一个值即可,如果,想获得10个周期的数据,则此时,是从0到9,即0代表第一个周期,9实际上是第10个周期的。
      如N=36时,可以写成7*5+1的形式,7是限制条件(不能取它的倍数),5是周期值,5个周期,实际求出来的是第六个周期上的数据,在周期矩阵与单位矩阵耦合时,周期矩阵的值(周期数),需要乘7(周期值,或不能取得某数的倍数中,某数,这里不能取7的倍数,某数即7了);但是,单位矩阵的合成值不变,直接参与运算,统计值,是得到这个合成值的方法数,即得到对应N的不定方程的解组数。
       第一步,统计用的是Excel中的COUNTIF(B$2:G$7,I2)计数函数
       以后各步,统计用的是Excel中的SUMIF(B$16:L$26,AA16,O$16:Y$26)计数函数

你的建议让我不知所措(让我先学代数,后学编程?无所适从)
#26
独木星空2022-11-03 21:21
好的编程都是建立在对问题的熟练程度上,对问题的来龙去脉不清楚,你就不可能有好的算法。
#27
吹水佬2022-11-03 22:50
以下是引用独木星空在2022-11-3 21:21:37的发言:

好的编程都是建立在对问题的熟练程度上,对问题的来龙去脉不清楚,你就不可能有好的算法。

电算与手算有时不是基于同一思路过程。
就如上次“开方”之说,模拟手算的套路可能要迭代或递归,但有人用一个常数系数算一句就搞定。
所以,电算有时不一定是按纯数理出牌,虽然没有纯数理的精确,但只要能高效达到实际要求的就是好牌。
从电算的角度看,算法有时不定是用数理推导出来,数理只能作指导。
程序 = 数据结构 + 算法

[此贴子已经被作者于2022-11-3 22:52编辑过]

#28
jiazhefish2022-11-04 08:18
以下是引用吹水佬在2022-11-3 22:50:15的发言:


电算与手算有时不是基于同一思路过程。
就如上次“开方”之说,模拟手算的套路可能要迭代或递归,但有人用一个常数系数算一句就搞定。
所以,电算有时不一定是按纯数理出牌,虽然没有纯数理的精确,但只要能高效达到实际要求的就是好牌。
从电算的角度看,算法有时不定是用数理推导出来,数理只能作指导。
程序 = 数据结构 + 算法

#29
独木星空2022-11-06 15:11
单位矩阵9        1        2        3        4        5        6        7        8
1        2        3        4        5        6        7        8        9
2        3        4        5        6        7        8        9        10
3        4        5        6        7        8        9        10        11
4        5        6        7        8        9        10        11        12
5        6        7        8        9        10        11        12        13
6        7        8        9        10        11        12        13        14
7        8        9        10        11        12        13        14        15
8        9        10        11        12        13        14        15        16

模9        统计2
2        1
3        2
4        3
5        4
6        5
7        6
8        7
9        8
10        7
11        6
12        5
13        4
14        3
15        2
16        1
合计        64
单位矩阵,第一步计算。
#30
独木星空2022-11-06 15:11
模9        1        2        3        4        5        6        7        8
2        3        4        5        6        7        8        9        10
3        4        5        6        7        8        9        10        11
4        5        6        7        8        9        10        11        12
5        6        7        8        9        10        11        12        13
6        7        8        9        10        11        12        13        14
7        8        9        10        11        12        13        14        15
8        9        10        11        12        13        14        15        16
9        10        11        12        13        14        15        16        17
10        11        12        13        14        15        16        17        18
11        12        13        14        15        16        17        18        19
12        13        14        15        16        17        18        19        20
13        14        15        16        17        18        19        20        21
14        15        16        17        18        19        20        21        22
15        16        17        18        19        20        21        22        23
16        17        18        19        20        21        22        23        24

统计2/1        1        1        1        1        1        1        1        1
1        1        1        1        1        1        1        1        1
2        2        2        2        2        2        2        2        2
3        3        3        3        3        3        3        3        3
4        4        4        4        4        4        4        4        4
5        5        5        5        5        5        5        5        5
6        6        6        6        6        6        6        6        6
7        7        7        7        7        7        7        7        7
8        8        8        8        8        8        8        8        8
7        7        7        7        7        7        7        7        7
6        6        6        6        6        6        6        6        6
5        5        5        5        5        5        5        5        5
4        4        4        4        4        4        4        4        4
3        3        3        3        3        3        3        3        3
2        2        2        2        2        2        2        2        2
1        1        1        1        1        1        1        1        1

模9        统计3
3        1
4        3
5        6
6        10
7        15
8        21
9        28
10        36
11        42
12        46
13        48
14        48
15        46
16        42
17        36
18        28
19        21
20        15
21        10
22        6
23        3
24        1
合计        512
单位矩阵,第二步计算2+1(2是列数据,1是行数据,交叉点是实际:数据区域,有行*列或得)
#31
独木星空2022-11-06 15:12
模9        3        4        5        6        7        8        9
3        6        7        8        9        10        11        12
4        7        8        9        10        11        12        13
5        8        9        10        11        12        13        14
6        9        10        11        12        13        14        15
7        10        11        12        13        14        15        16
8        11        12        13        14        15        16        17
9        12        13        14        15        16        17        18
10        13        14        15        16        17        18        19
11        14        15        16        17        18        19        20
12        15        16        17        18        19        20        21
13        16        17        18        19        20        21        22
14        17        18        19        20        21        22        23
15        18        19        20        21        22        23        24
16        19        20        21        22        23        24        25
17        20        21        22        23        24        25        26
18        21        22        23        24        25        26        27
19        22        23        24        25        26        27        28
20        23        24        25        26        27        28        29
21        24        25        26        27        28        29        30
22        25        26        27        28        29        30        31
23        26        27        28        29        30        31        32
24        27        28        29        30        31        32        33

统计3/3        1        3        6        10        15        21        28
1        1        3        6        10        15        21        28
3        3        9        18        30        45        63        84
6        6        18        36        60        90        126        168
10        10        30        60        100        150        210        280
15        15        45        90        150        225        315        420
21        21        63        126        210        315        441        588
28        28        84        168        280        420        588        784
36        36        108        216        360        540        756        1008
42        42        126        252        420        630        882        1176
46        46        138        276        460        690        966        1288
48        48        144        288        480        720        1008        1344
48        48        144        288        480        720        1008        1344
46        46        138        276        460        690        966        1288
42        42        126        252        420        630        882        1176
36        36        108        216        360        540        756        1008
28        28        84        168        280        420        588        784
21        21        63        126        210        315        441        588
15        15        45        90        150        225        315        420
10        10        30        60        100        150        210        280
6        6        18        36        60        90        126        168
3        3        9        18        30        45        63        84
1        1        3        6        10        15        21        28

模9        统计6
6        1
7        6
8        21
9        56
10        126
11        252
12        462
13        792
14        1281
15        1966
16        2877
17        4032
18        5432
19        7056
20        8856
21        10752
22        12642
23        14412
24        15946
25        17136
26        17892
27        18152
28        17892
29        17136
30        15946
31        14412
32        12642
33        10752
34        8856
35        7056
36        5432
37        4032
38        2877
39        1966
40        1281
41        792
42        462
43        252
44        126
45        56
46        21
47        6
48        1
合计        262144
单位矩阵,第三步计算,3+3;数据区域仅仅显示了7列数据,还有15列数据未显示。
#32
独木星空2022-11-06 15:12
模9        3        4        5        6        7        8        9
6        9        10        11        12        13        14        15
7        10        11        12        13        14        15        16
8        11        12        13        14        15        16        17
9        12        13        14        15        16        17        18
10        13        14        15        16        17        18        19
11        14        15        16        17        18        19        20
12        15        16        17        18        19        20        21
13        16        17        18        19        20        21        22
14        17        18        19        20        21        22        23
15        18        19        20        21        22        23        24
16        19        20        21        22        23        24        25
17        20        21        22        23        24        25        26
18        21        22        23        24        25        26        27
19        22        23        24        25        26        27        28
20        23        24        25        26        27        28        29
21        24        25        26        27        28        29        30
22        25        26        27        28        29        30        31
23        26        27        28        29        30        31        32
24        27        28        29        30        31        32        33
25        28        29        30        31        32        33        34
26        29        30        31        32        33        34        35
27        30        31        32        33        34        35        36
28        31        32        33        34        35        36        37
29        32        33        34        35        36        37        38
30        33        34        35        36        37        38        39
31        34        35        36        37        38        39        40
32        35        36        37        38        39        40        41
33        36        37        38        39        40        41        42
34        37        38        39        40        41        42        43
35        38        39        40        41        42        43        44
36        39        40        41        42        43        44        45
37        40        41        42        43        44        45        46
38        41        42        43        44        45        46        47
39        42        43        44        45        46        47        48
40        43        44        45        46        47        48        49
41        44        45        46        47        48        49        50
42        45        46        47        48        49        50        51
43        46        47        48        49        50        51        52
44        47        48        49        50        51        52        53
45        48        49        50        51        52        53        54
46        49        50        51        52        53        54        55
47        50        51        52        53        54        55        56
48        51        52        53        54        55        56        57

统计6/3        1        3        6        10        15        21        28
1        1        3        6        10        15        21        28
6        6        18        36        60        90        126        168
21        21        63        126        210        315        441        588
56        56        168        336        560        840        1176        1568
126        126        378        756        1260        1890        2646        3528
252        252        756        1512        2520        3780        5292        7056
462        462        1386        2772        4620        6930        9702        12936
792        792        2376        4752        7920        11880        16632        22176
1281        1281        3843        7686        12810        19215        26901        35868
1966        1966        5898        11796        19660        29490        41286        55048
2877        2877        8631        17262        28770        43155        60417        80556
4032        4032        12096        24192        40320        60480        84672        112896
5432        5432        16296        32592        54320        81480        114072        152096
7056        7056        21168        42336        70560        105840        148176        197568
8856        8856        26568        53136        88560        132840        185976        247968
10752        10752        32256        64512        107520        161280        225792        301056
12642        12642        37926        75852        126420        189630        265482        353976
14412        14412        43236        86472        144120        216180        302652        403536
15946        15946        47838        95676        159460        239190        334866        446488
17136        17136        51408        102816        171360        257040        359856        479808
17892        17892        53676        107352        178920        268380        375732        500976
18152        18152        54456        108912        181520        272280        381192        508256
17892        17892        53676        107352        178920        268380        375732        500976
17136        17136        51408        102816        171360        257040        359856        479808
15946        15946        47838        95676        159460        239190        334866        446488
14412        14412        43236        86472        144120        216180        302652        403536
12642        12642        37926        75852        126420        189630        265482        353976
10752        10752        32256        64512        107520        161280        225792        301056
8856        8856        26568        53136        88560        132840        185976        247968
7056        7056        21168        42336        70560        105840        148176        197568
5432        5432        16296        32592        54320        81480        114072        152096
4032        4032        12096        24192        40320        60480        84672        112896
2877        2877        8631        17262        28770        43155        60417        80556
1966        1966        5898        11796        19660        29490        41286        55048
1281        1281        3843        7686        12810        19215        26901        35868
792        792        2376        4752        7920        11880        16632        22176
462        462        1386        2772        4620        6930        9702        12936
252        252        756        1512        2520        3780        5292        7056
126        126        378        756        1260        1890        2646        3528
56        56        168        336        560        840        1176        1568
21        21        63        126        210        315        441        588
6        6        18        36        60        90        126        168
1        1        3        6        10        15        21        28

模9        统计9
9        1
10        9
11        45
12        165
13        495
14        1287
15        3003
16        6435
17        12861
18        24229
19        43353
20        74097
21        121515
22        191907
23        292743
24        432399
25        619677
26        863109
27        1170073
28        1545777
29        1992195
30        2507067
31        3083103
32        3707559
33        4362297
34        5024385
35        5667237
36        6262237
37        6780735
38        7196247
39        7486635
40        7635987
41        7635987
42        7486635
43        7196247
44        6780735
45        6262237
46        5667237
47        5024385
48        4362297
49        3707559
50        3083103
51        2507067
52        1992195
53        1545777
54        1170073
55        863109
56        619677
57        432399
58        292743
59        191907
60        121515
61        74097
62        43353
63        24229
64        12861
65        6435
66        3003
67        1287
68        495
69        165
70        45
71        9
72        1
合计        134217728
单位矩阵,第四步,6+3,数据区域仅仅显示了7列数据,还有15列数据未显示。
单位矩阵计算完成。  周期矩阵可以重复上述步骤,只不过从0周开始,算10个周期的吧,即参与元素是0至9的10个数。
#33
独木星空2022-11-06 15:17
9个未知数(或称,9个自变量)的线性不定方程,限制条件,自变量不能取9的倍数,求N对应的解组数。上边是用Excel手工计算的。单位矩阵计算步骤,周期矩阵计算步骤未列出,最后一步,单位矩阵与周期矩阵耦合也未列出。
#34
独木星空2022-11-07 13:30
在Excel上能实现的,在vfp上也一定可以实现,还有vb6上都可以实现。
#35
独木星空2022-11-19 16:29
在vb6版块,给出了多种解决方案。其中一人,还在Excel中建了宏进行了计算。
#36
独木星空2023-06-27 20:56
CLEAR ALL
SELECT 1
USE  d:\孪生素数\孪素19.DBF ALIAS 孪19表
  SELECT 2
USE  d:\孪生素数\孪素基表.DBF ALIAS 孪基表
    kssj=SECONDS()                      &&取出开始时间
   For N=1 to 19
    @ 3,6 say N
       SELECT 孪基表
       go 1
       for m=1 to 22275
           && @ 4,8 say 孪素
           X=孪素+510510*(N-1)
           MZ=X%19
                IF  MZ=0 OR MZ=2
                  else
                  SELECT   孪19表      &&打开盛放素数的表
                  APPEND BLANK                &&增加一条空记录
                  REPLACE 孪19 WITH X     &&将N值付给素
                ENDIF
                SELECT 孪基表
                 SKIP
       endfor
    ENDFOR
   =MESSAGEBOX("运行时间:"+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
这是用vfp程序通过孪素17基表制作孪素19生成元表。
发在此标题下并不合适,无奈找不到自己的贴切主题,就这样吧,主题太少了。
#37
shizi02023-06-27 22:20
如果纯讨论数学,把数学模型转变成程序,那这个程序的应用价值何在?因为旁观者要看明白、听明白楼主所述,耗时一定很久,更何况还要演算、编程、测试,所以在不明白价值的情况下,难有“闲人”深入。
比如说,楼主解决这问题是用于“围棋机器人”?“复杂排课算法”?.....
#38
独木星空2023-06-28 13:01
CLEAR ALL
SELECT 1
USE d:\孪生素数\素数表万.DBF ALIAS 万素
SELECT 2
USE d:\孪生素数\孪素19.DBF ALIAS 孪素19基表
SELECT 3
USE d:\孪生素数\孪生素数表.DBF ALIAS 孪素表
     kssj=SECONDS()                      &&取出开始时间
  For i=1 to 10
         @ 3,6 say i
         SELECT 孪素19基表
       go 1
       SELECT 孪素19基表
    for m=1 to 378675
       k=孪19
       N=k+9699690*(i-1)
        Kf=INT(SQRT(N))       &&求出奇数的开方根
        SELECT  万素             &&打开素数表
        LOCATE FOR 素数>=kf      &&根据开方根,查找最大素数
        DO CASE
                CASE EOF()          &&如果超出素数表最后一条记录
                        GO BOTTOM
                CASE 素数>kf     &&如果找到的素数比开方根值大,就将记录指针退回一个,保证使用的最大素数在开方根内
                        SKIP -1
        ENDCASE
        SSS=RECNO()    &&记录下最大素数(开方根内)记录位置(如果找到的素数等于开方根,则直接读取这个素数的位置)
        GO 9           &&从第一条记录开始读取素数(3)
        FOR j=1 TO SSS-8      &&内循环开始。这个循环实质上是从小到大顺序,依次读取素数。循环值是记录序号
                qmz=MOD(N,素数)   &&以读取的素数为条件,对被判断值求模 【qmz=MOD(bpz,素数)】文本号括着的原在此行语句体部分.特别需要注意的是在用标点符号时,一定改成英文模式及半角状态。
                IF qmz=0 OR qmz=2
                        EXIT         
                ENDIF
                SKIP            &&素数表指针向下移动一个
        ENDFOR
        IF j>sss-8
                SELECT  孪素表  &&打开盛放素数的表
                APPEND BLANK                        &&增加一条空记录
                REPLACE 孪素 WITH N             &&将N值付给素数
        ENDIF
        SELECT 2
        skip
      ENDFOR
   ENDFOR
        =MESSAGEBOX("运行时间:"+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
这是筛选10个周期中的孪生素数程序(周期值9699690).
有专业知识,又有vfp编程,那一方面都会影响到读者。
#39
独木星空2023-07-03 21:50
SELECT 1
USE H:\素数式系数\素数表十亿.dbf ALIAS 素数表
kssj=SECONDS()  &&取出开始时间
s=6.000000000000000000
SELECT 1
  GO 3
  A=素数
  s=s*(1+8/(A-2)^3)
  SKIP 1
        FOR  j=1 TO 50876310
        SELECT 1
        B=素数
        s=s*(1+8/(B-2)^3)
        SELECT 1  
        SKIP
        ENDFOR

?s
=MESSAGEBOX("运行时间:"+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
运行结果:6*∏\((1+{8\over (P-2)^3})\)=8.471329459766746000,(P≥5,P是素数)
没有合适的位置,就这么着了。
#40
独木星空2023-07-08 21:35
SELECT 1
USE h:\素数式系数\素数表十亿.DBF ALIAS 素数表
SELECT 2
USE h:\素数式系数\孪中十亿.DBF ALIAS 孪中表
kssj=SECONDS()
SELECT  1
GO 1
For i=1 to 50847534
      @ 5,12 say i
    SELECT  1     &&打开盛放素数式的表
    jl=recno()
    PUBLIC A
    A=素数
    SKIP 1
    PUBLIC B
    B=素数
    cj1=B-A  &&求前后两个素数的差值
    IF cj1=2
         SELECT 2
         APPEND BLANK     &&增加一条空记录
         REPLACE 孪中 WITH A+1          &&将N值付给素数式
    ENDIF
      SELECT 1
      GO jl+1
   ENDFOR
    =MESSAGEBOX("运行时间:"+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
用素数表制作孪中程序。
发的位置不对,只是作下记录,给用心者学习参考。
#41
独木星空2023-07-11 22:23
SELECT 1
USE h:\二生素数串\孪串13表.DBF ALIAS 孪串13表
SELECT 3
USE h:\二生素数串\孪串17表.DBF ALIAS 孪串17表
kssj=SECONDS()                      &&取出开始时间
  FOR i=1 TO 17
     @22,20 SAY i
     SELECT 1
     GO 1
       FOR k=1 TO 1485
       lc13=孪串13
       bpz=lc13+(i-1)*30030 &&给被判断值置数,以循环值做变量。
        ys=MOD(bpz,17)
        IF ys=0 OR ys=2
        else                        
        SELECT 3              &&打开保存求解结果的信息表
        APPEND BLANK                        &&增加一条空记录
        REPLACE 孪串17 WITH bpz             &&将被判断值保存到被判断值字段中(被判断值数据栏)
        ENDIF
        SELECT 1
        SKIP
        ENDFOR
ENDFOR
=MESSAGEBOX("运行时间:"+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
这是对孪生素数串的一个研究,不知放在那儿,吹水佬版主莫怪。
#42
独木星空2023-07-11 22:36
SELECT 1
USE h:\二生素数串\孪串17表.DBF ALIAS 孪串17表
SELECT 3
USE h:\二生素数串\孪串19表.DBF ALIAS 孪串19表
kssj=SECONDS()                      &&取出开始时间
  FOR i=1 TO 19
     @22,20 SAY i
     SELECT 1
     GO 1
       FOR k=1 TO 22275
       lc17=孪串17
       bpz=lc17+(i-1)*510510 &&给被判断值置数,以循环值做变量。
        ys=MOD(bpz,19)
        IF ys=0 OR ys=2
        else                        
        SELECT 3              &&打开保存求解结果的信息表
        APPEND BLANK                        &&增加一条空记录
        REPLACE 孪串19 WITH bpz             &&将被判断值保存到被判断值字段中(被判断值数据栏)
        ENDIF
        SELECT 1
        SKIP
        ENDFOR
ENDFOR
=MESSAGEBOX("运行时间:"+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
#43
独木星空2023-07-11 22:36
SELECT 1
USE h:\二生素数串\孪串19表.DBF ALIAS 孪串19表
SELECT 3
USE h:\二生素数串\孪串23表.DBF ALIAS 孪串23表
kssj=SECONDS()                      &&取出开始时间
  FOR i=1 TO 23
     @22,20 SAY i
     SELECT 1
     GO 1
       FOR k=1 TO 378675
       lc19=孪串19
       bpz=lc19+(i-1)*9699690 &&给被判断值置数,以循环值做变量。
        ys=MOD(bpz,23)
        IF ys=0 OR ys=2
        else                        
        SELECT 3              &&打开保存求解结果的信息表
        APPEND BLANK                        &&增加一条空记录
        REPLACE 孪串23 WITH bpz             &&将被判断值保存到被判断值字段中(被判断值数据栏)
        ENDIF
        SELECT 1
        SKIP
        ENDFOR
ENDFOR
=MESSAGEBOX("运行时间:"+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
#44
独木星空2023-07-11 22:43
今天喝多了,发了不着边际的程序,不为别的,为出书做准备。
[附件]2[/附件
看来超过5M咋也发不了。
#45
独木星空2023-07-12 00:08
SELECT 1
USE h:\二生素数串\孪串23表.DBF ALIAS 孪串23表
SELECT 3
USE h:\二生素数串\孪串3表.DBF ALIAS 孪串3表
kssj=SECONDS()                      &&取出开始时间
     SELECT 1
     GO 1
       FOR k=1 TO 7952173
        @22,20 SAY k
        jlh=recno()&&jlh记录号
        PUBLIC A
        A=孪串23
        SKIP 2
        PUBLIC B
        B=孪串23
        C=B-A                    
        SELECT 3              &&打开保存求解结果的信息表
        APPEND BLANK                        &&增加一条空记录
        REPLACE 串3 WITH C          &&将被判断值保存到被判断值字段中(被判断值数据栏)
        SELECT 1
        GO jlh+1
        ENDFOR
=MESSAGEBOX("运行时间:"+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
#46
独木星空2023-07-12 00:09
sELECT 1
USE h:\二生素数串\孪串23表.DBF ALIAS 孪串23表
SELECT 3
USE h:\二生素数串\孪串4表.DBF ALIAS 孪串4表
kssj=SECONDS()                      &&取出开始时间
     SELECT 1
     GO 1
       FOR k=1 TO 7952172
        @22,20 SAY k
        jlh=recno()&&jlh记录号
        PUBLIC A
        A=孪串23
        SKIP 3
        PUBLIC B
        B=孪串23
        C=B-A                    
        SELECT 3              &&打开保存求解结果的信息表
        APPEND BLANK                        &&增加一条空记录
        REPLACE 串4 WITH C          &&将被判断值保存到被判断值字段中(被判断值数据栏)
        REPLACE 记录号4 WITH jlh  
        SELECT 1
        GO jlh+1
        ENDFOR
=MESSAGEBOX("运行时间:"+LTRIM(STR(INT((SECONDS()-kssj)/60)))+"分"+LTRIM(STR(MOD(SECONDS()-kssj,60),5,2))+"秒",64,"运行时间提示")
#47
独木星空2023-07-12 00:12
孪串4重磅出炉。
只有本站会员才能查看附件,请 登录
#48
laowan0012023-07-12 07:40
这是运行了将近24个小时?是46楼的程序吗?


[此贴子已经被作者于2023-7-12 07:51编辑过]

#49
独木星空2023-07-12 08:18
回复 48楼 laowan001
这是因为跨了夜,也就几分钟的事。那是负时间。
#50
laowan0012023-07-12 08:26
#51
独木星空2023-07-12 20:39
如果不显示循环号,运算速度还可以。
只有本站会员才能查看附件,请 登录
12