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

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

 火... [复制链接]
发表于 2023-8-3 17:55:27 | 显示全部楼层
本帖最后由 northwolves 于 2023-8-3 18:06 编辑
王守恩 发表于 2023-8-3 17:51
我们这串数: 0, 4, 48, 1440, 65280, ..., 通项公式(A189849/2^n)与 A189849不大一样。
如果用这串数, ...


n对夫妻入住某酒店的n个标间,夫妻不得同屋
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2023-8-6 18:39:27 | 显示全部楼层
(1),
  1. Table[Minimize[{1/Sin[a] + n/Cos[a], \[Pi]/2 > a > 0}, {a}], {n, 1, 9}] // ToRadicals // FullSimplify
复制代码

   {2 Sqrt[2],   Sqrt[5 + 6 2^(1/3) + 3 2^(2/3)],    Sqrt[10 + 9 3^(1/3) + 3 3^(2/3)],    Sqrt[17 + 6 2^(1/3) + 12 2^(2/3)],
   Sqrt[26 + 15 5^(1/3) + 3 5^(2/3)],  Sqrt[37 + 18 6^(1/3) + 3 6^(2/3)],  Sqrt[50 + 21 7^(1/3) + 3 7^(2/3)],  5 Sqrt[5],   Sqrt[82 + 9 3^(1/3) + 27 3^(2/3)]}
(1)可以有通项公式:
  1. Table[Solve[{s == Sqrt[n^2 + 1 + 3 n Power[n, (3)^-1] + 3 Power[n^2, (3)^-1]]}, {s}], {n, 1, 9}] // ToRadicals // FullSimplify
复制代码

s -> 2 Sqrt[2], s -> Sqrt[5 + 6 2^(1/3) + 3 2^(2/3)], s -> Sqrt[10 + 9 3^(1/3) + 3 3^(2/3)], s -> Sqrt[17 + 6 2^(1/3) + 12 2^(2/3)],
s -> Sqrt[26 + 15 5^(1/3) + 3 5^(2/3)], s -> Sqrt[37 + 18 6^(1/3) + 3 6^(2/3)], s -> Sqrt[50 + 21 7^(1/3) + 3 7^(2/3)], s -> 5 Sqrt[5], s -> Sqrt[82 + 9 3^(1/3) + 27 3^(2/3)]}
(2),
  1. Table[Minimize[{2/Sin[a] + n/Cos[a], \[Pi]/2 > a > 0}, {a}], {n, 1, 9}] // ToRadicals // FullSimplify
复制代码

   {Sqrt[5 + 6 2^(1/3) + 3 2^(2/3)],  4 Sqrt[2],  Root[-2197 - 465 #^2 - 39 #^4 + #^6&,2,0],  2 Sqrt[5 + 6 2^(1/3) + 3 2^(2/3)],  {Root[-24389 - 177 #^2 - 87 #^4 + #^6&,2,0],
   2 Sqrt[10 + 9 3^(1/3) + 3 3^(2/3)],  Sqrt[53 + 21 2^(2/3) 7^(1/3) + 6 2^(1/3) 7^(2/3)],  2 Sqrt[17 + 6 2^(1/3) + 12 2^(2/3)],  Sqrt[85 + 18 6^(1/3) + 27 6^(2/3)]}
(3),
  1. Table[Minimize[{3/Sin[a] + n/Cos[a], \[Pi]/2 > a > 0}, {a}], {n, 1, 9}] // ToRadicals // FullSimplify
复制代码

   {Sqrt[10 + 9 3^(1/3) + 3 3^(2/3)],  Root[-2197 - 465 #^2 - 39 #^4 + #^6&,2,0],  6 Sqrt[2],   Sqrt[25 + 18 6^(1/3) + 12 6^(2/3)],  Sqrt[34 + 15 3^(2/3) 5^(1/3) + 9 3^(1/3) 5^(2/3)],
   3 Sqrt[5 + 6 2^(1/3) + 3 2^(2/3)],   Sqrt[58 + 21 3^(2/3) 7^(1/3) + 9 3^(1/3) 7^(2/3)],   Sqrt[73 + 36 3^(1/3) + 48 3^(2/3)],   3 Sqrt[10 + 9 3^(1/3) + 3 3^(2/3)]]}
......
(1)可以有通项公式;   (2),(3),(4)...可以有通项公式吗?
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2023-8-7 10:10:04 | 显示全部楼层
本帖最后由 王守恩 于 2023-8-7 10:16 编辑

(1),(2),(3),(4)...可以有统一的通项公式。

若  \(\pi/2 > a > 0,\frac{m}{\sin(a)}+\frac{n}{\cos(a)}\)  的最小值=\(\D(m^{2/3}+n^{2/3})^{3/2}\)

评分

参与人数 1威望 +8 金币 +8 贡献 +8 经验 +8 鲜花 +8 收起 理由
northwolves + 8 + 8 + 8 + 8 + 8 赞一个!

查看全部评分

毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2023-8-7 10:16:41 | 显示全部楼层
本帖最后由 northwolves 于 2023-8-7 13:44 编辑

$f(a)=k/Sin[a] + n/Cos[a]$
$f'(a)=\frac{kTan[a]}{Cos[a]}- \frac{n}{Tan[a]Sin[a]}$
$f'(a)=0->Tan[a]=(\frac{k}{n})^{\frac{1}{3}}$
此时 $k/Sin[a] + n/Cos[a]=(k+nTan[a])\sqrt{1+\frac{1}{Tan^2[a])}=(k^{\frac{2}{3}}+n^{\frac{2}{3}})^{\frac{3}{2}}$

评分

参与人数 1威望 +9 金币 +9 贡献 +9 经验 +9 鲜花 +9 收起 理由
王守恩 + 9 + 9 + 9 + 9 + 9 很给力!

查看全部评分

毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2023-8-9 10:32:54 | 显示全部楼层
有这样一串数:4, 8, 9, 16, 25, 27, 32, 36, 49, 64, 81, 100, 121, 125, 128, 144, 169, 196, 216, 225,...A001597
1,每个数-1,取倒数相加,和是多少?1/3+1/7+1/8+1/15+1/24+1/26+1/31+1/35+1/48+1/63+1/80+1/99+1/120+...
2,每个数-指数(较小者),取倒数相加,和是多少?1/2+1/5+1/7+1/14+1/23+24+1/27+1/34+1/47+1/62+1/79+1/98+...       
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2023-8-9 16:23:04 | 显示全部楼层
$\sum_{k=2}^{∞}\frac{1}{n^k}=\frac{1}{n-1}-\frac{1}{n}$
$1/3+1/7+1/8+1/15+1/24+1/26+1/31+1/35+1/48+1/63+1/80+1/99+1/120+...=1-\frac{1}{2}+\frac1 2-\frac1 3+...=1$

点评

https://www.maa.org/sites/default/files/pdf/upload_library/22/Ford/bibiloni206.pdf  发表于 2023-8-13 16:48
Goldbach showed that Sum 1/a(n) = 1, see A214390, A214391.  发表于 2023-8-13 12:56
Table[{m,N[Total[1/Union@Select[Flatten@Table[n^k-1,{n,2,10^m},{k,2,Floor[m*Log[100]/Log[n]]}],#<100^m&]],20]},{m,6}]  发表于 2023-8-10 15:39

评分

参与人数 1威望 +12 金币 +12 贡献 +12 经验 +12 鲜花 +12 收起 理由
王守恩 + 12 + 12 + 12 + 12 + 12 我就是缺少这些按钮!

查看全部评分

毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2023-8-9 16:31:48 | 显示全部楼层
  1. N[Total[1/Union@Select[Flatten@Table[n^k-k,{n,2,10000},{k,2,Floor[Log[10^8]/Log[n]]}],#<10^8&]],20]
复制代码


1.4085485959754331397
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2023-8-9 16:32:43 | 显示全部楼层
第2题猜测极限值$\sqrt2$
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2023-8-25 15:48:02 | 显示全部楼层
一个数列的通项式,据说难倒英雄汉!
每一行的和都是"兔子数列"。每一列都是"杨辉三角"中的列。
  1. Cn:Table[(n+k)!/((2k)!(n-k)!),{n,0,9},{k,0,n}]//TableForm
复制代码

{1},
{1, 1},
{1, 3, 1},
{1, 6, 5, 1},
{1, 10, 15, 7, 1},
{1, 15, 35, 28, 9, 1},
{1, 21, 70, 84, 45, 11, 1},
{1, 28, 126, 210, 165, 66, 13, 1},
{1, 36, 210, 462, 495, 286, 91, 15, 1},
{1, 45, 330, 924, 1287, 1001, 455, 120, 17, 1}
  1. An:Table[(n+k)!/((2k)!(n-k)!),{n,1,9},{k,0,n}]//TableForm
复制代码

{1, 1},
{1, 3, 1},
{1, 6, 5, 1},
{1, 10, 15, 7, 1},
{1, 15, 35, 28, 9, 1},
{1, 21, 70, 84, 45, 11, 1},
{1, 28, 126, 210, 165, 66, 13, 1},
{1, 36, 210, 462, 495, 286, 91, 15, 1},
{1, 45, 330, 924, 1287, 1001, 455, 120, 17, 1}}]
  1. Dn:Table[(n+k+1)!/((2k1)!(n-k)!),{n,0,9},{k,0,n}]//TableForm
复制代码

{1},
{2, 1},
{3, 4, 1},
{4, 10, 6, 1},
{5, 20, 21, 8, 1},
{6, 35, 56, 36, 10, 1},
{7, 56, 126, 120, 55, 12, 1},
{8, 84, 252, 330, 220, 78, 14, 1},
{9, 120, 462, 792, 715, 364, 105, 16, 1},
{10,165, 792, 1716, 2002, 1365, 560, 136, 18, 1}}]
  1. Bn:Table[(n+k+1)!/((2k+1)!(n-k)!),{n,1,9},{k,0,n}]//TableForm
复制代码

{2, 1},
{3, 4, 1},
{4, 10, 6, 1},
{5, 20, 21, 8, 1},
{6, 35, 56, 36, 10, 1},
{7, 56, 126, 120, 55, 12, 1},
{8, 84, 252, 330, 220, 78, 14, 1},
{9, 120, 462, 792, 715, 364, 105, 16, 1},
{10,165, 792, 1716, 2002, 1365, 560, 136, 18, 1}

点评

nyy
老同志  发表于 2023-8-25 15:54
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2023-8-28 17:57:19 | 显示全部楼层
用足够多的a,b,c凑出n,  一共有多少种方法?
  1. CoefficientList[Series[1/((1-x^a)(1-x^b)(1-x^c)),{x,0,n]],x]
复制代码
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

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

GMT+8, 2024-5-8 12:07 , Processed in 0.067595 second(s), 16 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

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