找回密码
 欢迎注册
楼主: 王守恩

[原创] 数字串的通项公式

 火... [复制链接]
 楼主| 发表于 2023-7-9 15:11:20 | 显示全部楼层
{1, 1, 2, 4, 7, 11, 16, 20, 21, 27, 34, 42, 51, 61, 72, 78, 85, 97, 110, 124, 139, 151, 160, 176, 189, 207, 226, 242, 263, 283, 296, 318, 337, 361, 386, 408, 435, 453, 476, 504}
  1. Table[Length@Flatten[Table[NSolve[{m/Sin[a]==n/Sin[b]==k/Sin[a+b]≥(mn/k)/Sin[Pi/3],Pi/3≥a>0,Pi/3≥b>0},{a,b}],{m,1,k},{n,1,k}]]/2,{k,1,30}]
复制代码

算式没问题,就是慢了一点。

点评

用NSolve实数范围内解方程,能快就怪了,速度估计会有至少几个数量级的差别  发表于 2023-7-9 16:35
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2023-7-9 17:29:21 | 显示全部楼层
更直观更高效的计算:
  1. Table[2Total@Table[Length@Select[Range[n-k+1,n],0<#^2-k^2<=n(n-k)&],{k,n}]+Ceiling[n/2],{n,176}]
复制代码


{1,1,2,4,7,11,16,20,21,27,34,42,51,61,72,78,85,97,110,124,139,151,160,176,189,207,226,242,263,283,296,318,337,361,386,408,435,453,476,504,525,555,578,610,639,663,692,726,753,789,822,860,887,921,956,992,1029,1071,1110,1144,1179,1219,1256,1298,1341,1385,1430,1466,1507,1553,1596,1648,1697,1751,1794,1836,1887,1935,1988,2046,2093,2153,2190,2242,2295,2353,2412,2476,2537,2585,2638,2696,2759,2823,2888,2958,3017,3067,3126,3190,3255,3325,3392,3464,3525,3581,3650,3724,3799,3875,3944,4024,4081,4155,4230,4310,4387,4465,4548,4618,4681,4761,4842,4924,5003,5095,5172,5252,5329,5419,5498,5586,5671,5765,5844,5930,6017,6105,6194,6288,6379,6473,6556,6644,6733,6835,6934,7030,7131,7223,7312,7414,7509,7617,7722,7828,7923,8021,8116,8216,8313,8427,8542,8650,8739,8843,8956,9062,9169,9289,9410,9510,9611,9725,9840,9956}
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2023-7-9 17:39:02 | 显示全部楼层
本帖最后由 northwolves 于 2023-7-9 18:00 编辑

或者
  1. Table[2Total@Table[Length@Select[Range[Max[k,n-k]+1,n],#^2+n*k<=n^2+k^2&],{k,n}]+Ceiling[n/2],{n,176}]
复制代码

  1. Table[2Total@Table[Floor@Sqrt[n^2+k^2-n*k]-Max[k,n-k],{k,n}]+Ceiling[n/2],{n,176}]
复制代码

  1. Table[2Sum[Floor@Sqrt[n^2+k^2-n*k]-Max[k,n-k],{k,n}]+Ceiling[n/2],{n,176}]
复制代码

  1. Table[2Sum[Floor@Sqrt[n^2+k^2-n*k],{k,n}]-n (3n-1)/2,{n,176}]
复制代码

点评

每道题(11道)中的#[[1]],#[[2]]都可以去掉吗?  发表于 2023-7-10 12:11

评分

参与人数 1威望 +24 金币 +24 贡献 +24 经验 +24 鲜花 +24 收起 理由
王守恩 + 24 + 24 + 24 + 24 + 24 吓坏了!

查看全部评分

毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2023-7-9 19:32:48 | 显示全部楼层
332:
  1. Table[2Total@Table[Length@Select[Range[n-k,n],(n-k)n≥#^2-k^2>0&],{k,n}]-Floor[(n-3)/2],{n,176}]
复制代码
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2023-7-10 16:33:00 | 显示全部楼层
每道题(11道)中的#[[1]],#[[2]]都可以去掉吗? 1—4题勉强可以。第5题不行, 第6题也不行, 7—11题更不行了。
  1. Table[2Total@Table[Length@Select[Range[n-k+1,n],(n-2Cos[a\[Pi]/12] k) n ≥ #^2-k^2>0 &],{k,n}]+Floor[(n Sec[a\[Pi]/12])/2]-Floor[n/2],{a,1,4},{n,30}]
复制代码

{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0},
{0, 0, 0, 0, 0, 0, 1, 0, 1, 2, 3, 4, 5, 5, 7, 7, 9, 11, 13, 15, 16, 19, 20, 23, 26, 30, 32, 32, 34, 38},
{0, 0, 1, 0, 1, 3, 5, 5, 8, 12, 12, 16, 21, 24, 25, 31, 38, 39, 44, 52, 52, 58, 67, 72, 79, 87, 98, 101, 110, 118},
{1, 1, 2, 4, 7, 11, 16, 20, 21, 27, 34, 42, 51, 61, 72, 78, 85, 97, 110, 124, 139, 151, 160, 176, 189, 207, 226, 242, 263, 283}
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2023-7-13 11:15:27 | 显示全部楼层
本帖最后由 王守恩 于 2023-7-13 16:52 编辑

正三角形(边长为整数)内(只要不超出边都可以)摆放三角形ABC(三边为整数), 这样的ABC有几个?

a(1)=01:{1,1,1},
a(2)=02:{1,1,1},{2,2,2},
a(3)=05:{1,1,1},{1,2,2},{2,2,2},{2,2,3},{3,3,3},
a(4)=10:{1,1,1},{1,2,2},{1,3,3},{2,2,2},{2,2,3},{2,3,3},{2,3,4},{3,3,3},{3,3,4},{4,4,4},
a(5)=18:{1,1,1},{1,2,2},{1,3,3},{1,4,4},{2,2,2},{2,2,3},{2,3,3},{2,3,4},{2,4,4},{2,4,5},{3,3,3},{3,3,4},{3,3,5},{3,4,4},{3,4,5},{4,4,4},{4,4,5},{5,5,5},

得到一串数。
{1, 2, 5, 10, 18, 29, 44, 60, 80, 106, 139,

点评

没有细致读过,进不了门。  发表于 2023-7-13 15:46
nyy
天灭你  发表于 2023-7-13 11:32
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2023-7-14 14:03:05 | 显示全部楼层
正三角形(边长为整数)内(只要不超出边都可以)摆放三角形ABC(三边为整数), 这样的ABC有几个?

a(1)=01:{1,1,1},
a(2)=02:{1,1,1},{2,2,2},
a(3)=05:{1,1,1},{1,2,2},{2,2,2},{2,2,3},{3,3,3},
a(4)=10:{1,1,1},{1,2,2},{1,3,3},{2,2,2},{2,2,3},{2,3,3},{2,3,4},{3,3,3},{3,3,4},{4,4,4},
a(5)=18:{1,1,1},{1,2,2},{1,3,3},{1,4,4},{2,2,2},{2,2,3},{2,3,3},{2,3,4},{2,4,4},{2,4,5},{3,3,3},{3,3,4},{3,3,5},{3,4,4},{3,4,5},{4,4,4},{4,4,5},{5,5,5},

得到一串数, (1)+(2)。
  1, 2, 5, 10, 18, 29, 44, 60, 80, 106, 139, 178, 224, 277, 335, 395, 467,549, 643, 746,

说明: 三角形ABC,   A,B,C的对边为a,b,c,   约定a≤b≤c,
A总可以与正三角形某顶点重合,AB总可以与正三角形某边重合,分2种情形统计
(1)   Pi/3≥B>0 :{1, 2, 4, 7, 12, 19, 29, 41, 54, 70, 90, 114, 143, 177, 217, 260, 307, 360, 420, 487}
  1. Table[Length@Flatten[Table[NSolve[{a/Sin[A]==b/Sin[B]==c/Sin[A+B],Pi/3≥A>0,Pi/3≥B>0},{A,B}],{a,1,n},{b,a,n},{c,b,n}]]/2,{n,1,20}]
复制代码

(2)Pi/2≥B>Pi/3:{0, 0, 1, 3,  6, 10, 15, 19, 26, 36, 49,   64,   81, 100, 118, 135, 160, 189, 223, 259}
  1. Table[Length@Flatten[Table[NSolve[{a/Sin[A]==b/Sin[B]==c/Sin[A+B],((n-c)Sqrt[3])/(2Sin[B-Pi/3])≥a,Pi/3≥A>0,Pi/2≥B>Pi/3},{A,B}],{a,1,n},{b,a+1,n},{c,b,n}]]/2,{n,1,20}]
复制代码

(1)=(3)前若干项的和
(3)1, 1, 2, 3, 5, 7, 10, 12, 13, 16, 20, 24, 29, 34, 40, 43, 47, 53, 60, 67, 75, 81, 86, 94, 101, 110, 120, 128, 139, 149}
Table[Sum[Length@Select[Range[n-k,n],(n-k)n≥#^2-k^2>0&],{k,0,n}],{n,1,30}]
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2023-7-15 19:04:13 | 显示全部楼层
一个正方形(边长=2,4,6,8,..., n+1个点)格点纸(1*1方格),
去掉一个最大的圆(半径=1,2,3,4,...),还剩几个格点(包括圆周上的格点)?
a(1)=4*2,
a(2)=4*4,
a(3)=4*6,
a(4)=4*9,
a(5)=4*11,

得到这样一串数: {2, 4, 6, 9, 11, 15, 20, 24, 28, 32, 39, 47, 51, 58, 64, 74, 82, 91, 99, 107}
  1. Table[n^2+n-Sum[SquaresR[2,k],{k,2,n^2}]/4,{n,1,20}]
复制代码
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2023-7-15 20:35:23 | 显示全部楼层
经久不衰的话题,不进门,难懂其道。
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2023-7-19 10:32:50 | 显示全部楼层
A135358       
       
6, 12, 19, 25, 32, 38, 45, 51, 58, 64, 71, 77, 83, 90, 96, 103, 109, 116, 122, 129, 135, 142, 148, 154,

\(a(n)=\frac{13n-GCD(n,n+2)}{2}\)   这样简单些。

点评

已达到,炉火纯青的地步!  发表于 2023-7-19 21:28
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

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

GMT+8, 2024-11-22 16:01 , Processed in 0.029134 second(s), 16 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

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