数学研发论坛

 找回密码
 欢迎注册
查看: 160|回复: 6

[讨论] 一类线性不定方程 \(x_1+2x_2+\cdots+kx_k=n\)非负整数解的数目

[复制链接]
发表于 2018-1-17 16:12:53 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有帐号?欢迎注册

x
一个猜想:
记方程$x_1+2x_2+\cdots+kx_k=n$的非负整数解个数为$f(n)$.
定义$g(q)=\frac{q^{-n-1}}{(1-q)(1-q^2)……(1-q^k)}$,

及$u(n)=Res[g(q),1],v(n)=Res[g(q),-1]$.

那么$f(n)=\floor(u(n)+v(n)-u(0)-(-1)^nv(0)+1)$.

点评

备注:本主题由 https://bbs.emath.ac.cn/thread-15213-1-1.html 分割而来  发表于 2018-9-12 14:55
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2018-1-17 18:04:21 | 显示全部楼层
用$f_k(n)$表示方程$x_1+2x_2+……+kx_k=n$的非负整数解的个数.
根据上面的猜想,有:
$f_1(n)=1.$
$f_2(n)=\floor(n/2+1].$
$f_3(n)=\floor(n^2/12+n/2+1].$
$f_4(n)=\floor(n^3/144+(5n^2)/48+((15+(-1)^n)n)/32+1].$
这几个公式已经验证过了,下面的公式有待验证:
$f_5(n)=\floor(n^4/2880+n^3/96+(31n^2)/288+((85+3(-1)^n)n)/192+1].$
$f_6(n)=\floor(n^5/86400 + (7 n^4)/11520 + (77 n^3)/6480 + (245 n^2)/2304 + (43981 n)/103680+(-1)^(-n) (n^2/768 + (7 n)/256) +1)$
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2018-1-17 19:27:59 | 显示全部楼层
由定义知,以上函数应满足递推关系式: $f_k(n)=f_{k-1}(n)+f_{k-1}(n-k)+f_{k-1}(n-2k)+\cdots$.
补充定义$f_k(n)=0(n<0).$
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2018-1-19 11:49:59 | 显示全部楼层
1#和2#猜想错误,k>5时会出现更多的重因子.
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2018-9-11 12:02:39 | 显示全部楼层
本帖最后由 葡萄糖 于 2018-9-11 17:43 编辑
lsr314 发表于 2018-1-16 13:09
$Res(q=1)=\frac{6n^2-6(a+b+c)n+a^2+b^2+c^2+3ab+3bc+3ca}{12abc}.$
令$\theta_a^k=e^{\frac{2\pi ik}{a}},(k=1,2,\cdots,a-1)$,如果a,b,c两两互素,则
$Res(q=\theta_a^k)=\frac{\theta_a^{b+c-n}}{a(1-\theta_a^b)(1-\theta_a^c)}.$
$f(n)=Req(q=1)+\sum_{k=1}^{a-1}Res(q=\theta_a^k)+\sum_{k=1}^{b-1}Res(q=\theta_b^k)+\sum_{k=1}^{c-1}Res(q=\theta_c^k)$
$f(n)=\frac{6n^2-6(a+b+c)n+a^2+b^2+c^2+3ab+3bc+3ca}{12abc}$
             $+\sum_{k=1}^{a-1}\frac{\theta_a^{b+c-n}}{a(1-\theta_a^b)(1-\theta_a^c)}+\sum_{k=1}^{b-1}\frac{\theta_b^{a+c-n}}{b(1-\theta_b^a)(1-\theta_b^c)}+\sum_{k=1}^{c-1}\frac{\theta_c^{b+a-n}}{c(1-\theta_c^b)(1-\theta_c^a)}$...


关于这个问题我暑假看到一篇小松尚夫(Takao Komatsu)写于2003年关于不定方程非负整数解数目的论文
On the number of solutions of the Diophantine equation of Frobenius – General case∗
https://core.ac.uk/download/pdf/14375587.pdf
好像就是用这个思路的
注:小松尚夫(Takao Komatsu)是武汉大学数学与统计学院聘请的外籍教授

另外希望把该帖子的标题修改为“一类线性不定方程$x_1+2x_2+……+kx_k=n$非负整数解的数目”


这是笔者根据小松尚夫的公式推出的一些三元情形下的结果(注:他给的公式只适用于线性不定方程的系数都互素的情形)
Let \(N\left(a_1,a_2,a_3\,;b\right)\) denote the number of  non-negative integer solutions of the equation:
\[a_1x_1+a_2x_2+a_3x_3=b\]
where  \(a_1,a_2,a_3,b\in\mathbb{Z_+}\) and \(\gcd(a_1,a_2)=\gcd(a_2,a_3)=\gcd(a_1,a_3)=1\)
\begin{align*}
N\left(1,2,3\,;b\right)&=\frac{1}{12}b^2+\frac{1}{2}b+\frac{47}{72}+\frac{1}{8}\cos(\pi b)+\frac{2}{9}\cos\left(\frac{2\pi}{3}b\right)\\
\\
N\left(1,2,5\,;b\right)&=\frac{1}{20}b^2+\frac{2}{5}b+\frac{27}{40}+\frac{1}{8}\cos(\pi b)+\frac{2\sqrt{5}}{25}\cos\left[\frac{2\pi}{5}(b-1)\right]+\frac{2\sqrt{5}}{25}\cos\left[\frac{\pi}{5}(4b+1)\right]\\
\\
N\left(1,3,4\,;b\right)&=\frac{1}{24}b^2+\frac{1}{3}b+\frac{83}{144}+\frac{2}{9}\cos\left[\frac{\pi}{3}(2b-1)\right]+\frac{1}{4}\cos\left(\frac{\pi}{2}b\right)+\frac{1}{16}\cos(\pi b)\\
\\
N\left(1,3,5\,;b\right)&=\frac{1}{30}b^2+\frac{3}{10}b+\frac{26}{45}+\frac{2}{9}\cos\left(\frac{2\pi}{3}b\right)+\frac{2\sqrt{5}}{25}\cos\left[\frac{\pi}{5}(2b-1)\right]+\frac{2\sqrt{5}}{25}\cos\left[\frac{2\pi}{5}(2b-1)\right]\\
\\
N\left(1,4,5\,;b\right)&=\frac{1}{40}b^2+\frac{1}{4}b+\frac{43}{80}+\frac{1}{4}\cos\left[\frac{\pi}{2}(b-1)\right]+\frac{1}{16}\cos(\pi b)+\frac{8}{25}\sin^2\left(\frac{2\pi}{5}\right)\cos\left(\frac{2\pi}{5}b\right)+\frac{8}{25}\sin^2\left(\frac{\pi}{5}\right)\cos\left(\frac{4\pi}{5}b\right)\\
\\
N\left(2,3,5\,;b\right)&=\frac{1}{60}b^2+\frac{1}{6}b+\frac{131}{360}+\frac{1}{8}\cos(\pi b)+\frac{2}{9}\cos\left[\frac{\pi}{3}(2b+1)\right]+\frac{8}{25}\sin^2\left(\frac{\pi}{5}\right)\cos\left(\frac{2\pi}{5}b\right)+\frac{8}{25}\sin^2\left(\frac{2\pi}{5}\right)\cos\left(\frac{4\pi}{5}b\right)\\
\\
N\left(3,4,5\,;b\right)&=\frac{1}{120}b^2+\frac{1}{10}b+\frac{191}{720}+\frac{2}{9}\cos\left(\frac{2\pi}{3}b\right)+\frac{1}{4}\cos\left(\frac{\pi}{2}b\right)+\frac{1}{16}\cos(\pi b)+\frac{2\sqrt{5}}{25}\cos\left[\frac{2\pi}{5}(b+1)\right]+\frac{2\sqrt{5}}{25}\cos\left[\frac{\pi}{5}(4b-1)\right]
\end{align*}
自己也整理了一些有公因子的情形:
where  \(a_1,a_2,a_3,b\in\mathbb{Z_+}\) and \(\exists\gcd(a_i,a_j)\ne1\,(i,j=1,2,3)\)
\begin{align*}
N\left(1,2,4\,;b\right)&=\frac{1}{16}b^2+\frac{7}{16}b+\frac{21}{32}+{\color{red}{\frac{1}{16}(b+1)\cos(\pi b)}}+\frac{5}{32}\cos(\pi b)+\frac{1}{8}\cos\left[\frac{\pi}{2}(b-1)\right]+\frac{1}{8}\cos\left(\frac{\pi}{2}b\right)\\
\\
N\left(2,3,4\,;b\right)&=\frac{1}{48}b^2+\frac{3}{16}b+\frac{107}{288}+{\color{red}{\frac{1}{16}(b+1)\cos(\pi b)}}+\frac{7}{32}\cos(\pi b)+\frac{2}{9}\cos\left(\frac{2\pi}{3}b\right)+\frac{1}{8}\cos\left(\frac{\pi}{2}b\right)+\frac{1}{8}\cos\left[\frac{\pi}{2}(b+1)\right]\\
\\
N\left(2,4,5\,;b\right)&=\frac{1}{80}b^2+\frac{11}{80}b+\frac{53}{160}+{\color{red}{\frac{1}{16}(b+1)\cos(\pi b)}}+\frac{9}{32}\cos(\pi b)+\frac{1}{8}\cos\left[\frac{\pi}{2}(b-1)\right]+\frac{1}{8}\cos\left(\frac{\pi}{2}b\right)+\frac{2\sqrt{5}}{25}\cos\left[\frac{\pi}{5}(2b+1)\right]+\frac{2\sqrt{5}}{25}\cos\left[\frac{2\pi}{5}(2b+1)\right]\\
\end{align*}
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2018-9-11 16:44:04 | 显示全部楼层
是否对于任意正整数$n$,总存在正整数$a,b,c,$使得$n^3+a^3=b^3+c^3$,且$b ne n,c ne n,gcd(a,b,c)=1$?
求解$a,b,c$的代码:
  1. f[n_] := (g = {};
  2.   Do[a = (3 x^2 + Sqrt[3] Sqrt[4 c^3 x - 4 n^3 x - x^4])/(6 x);
  3.    b = 1/6 (-3 x + (Sqrt[3] Sqrt[-x (-4 c^3 + 4 n^3 + x^3)])/x);
  4.    If[a != n && b != n && GCD[a, b, c] == 1, g = {n, a, b, c};
  5.     Break[]], {c, n + 1, 10^5}, {x,
  6.     Select[Divisors[
  7.       c^3 - n^3], (#^3 < c^3 - n^3 &&
  8.         IntegerQ[Sqrt[3] Sqrt[-# (-4 c^3 + 4 n^3 + #^3)]]) &]}]; g)
复制代码

找到公式了……

补充内容 (2018-9-12 16:26):
直接删除就行了。不过,为什么要限制发帖后可以编辑的时间呢,我觉得math.se这样随时可以编辑更好,甚至可以像math.se一样每个人都可以在登陆的情况下编辑别人的帖子,比如编辑公式之类的,发帖人可以撤回编辑就好了

点评

还是别在这个贴子发那个三次不定方程吧?不会显得太乱了吗?  发表于 2018-9-11 17:41
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

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

GMT+8, 2018-11-18 10:09 , Processed in 0.091771 second(s), 22 queries .

Powered by Discuz! X3.4

© 2001-2017 Comsenz Inc.

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