| 查看: 2077 | 回复: 11 | ||
theday0323木虫 (小有名气)
|
[求助]
急：求助求解不等式方程组 已有1人参与
|
|
|
N1，N2，N3，N4，N5，N6都是正整数 N1+N3+N5=N2+N4+N6; N1-N2>=0; N1-N2+N3-N4>=0; N1<=N; N1-N2+N3<=N; N1-N2+N3-N4+N5<=N; 求（N1+N2+N3+N4+N5+N6）/2的值分别在N到3N时解的个数，防止计算量太大，假定N为100. 求程序啊，急 |
回复此楼» 猜你喜欢
人气不行了
已经有9人回复
-
没有任何消息-是不是就凉了
已经有6人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有5人回复
-
2026国自然函评费到账
已经有19人回复
-
能否退出参与的面上项目解除限项
已经有21人回复
-
建议基金发布提前给出明确的时间点
已经有16人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有21人回复
-
估计是周四
已经有4人回复
-
我面上完蛋了
已经有9人回复
-
范进中举一文的中心思想
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助——求解三元二次多项式方程组
已经有9人回复
- matlab解方程遇到的
已经有7人回复
- 一个简单的运算，求助数学大神啊
已经有7人回复
- 求教一个关于不等式的问题!
已经有5人回复
- 求助，帮忙求解下面方程组，谢谢！
已经有13人回复
- mathematica求和为什么不出结果？
已经有7人回复
- 1stopt解方程组
已经有12人回复
- matlab解方程组求助攻
已经有4人回复
- 用Matlab的LMI工具计算线性不等式组
已经有6人回复
- 请教函数求解问题
已经有11人回复
- 半年了，还是解不了这个matlab..... 总是出错！！！
已经有18人回复
- 求教matlab线性优化求最小值的问题
已经有7人回复
- 请教MATLAB中的LMI求解问题，急！！！
已经有7人回复
- 求解一组四元一次方程
已经有4人回复
- 【求助】非线性方程组的求解问题
已经有6人回复
- 大神进，数学问题。。求解。
已经有7人回复
- matlab解不等式组
已经有6人回复
- mathematica 解方程组 有什么命令可以实现解集中没有复数解
已经有3人回复
- 求助解一道高中数学方程题
已经有8人回复
- 【求助】多元方程组求解及拟合
已经有9人回复
- 【求助】求助mathematica软件解方程，见图片
已经有5人回复
- 【求助】求解一元三次不等式
已经有4人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
walk1997
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4677.7
- 红花: 22
- 帖子: 1066
- 在线: 798.1小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
★ ★
theday0323(jjdg代发): 金币+2, 辛苦了 2014-06-11 12:25:31
theday0323(jjdg代发): 金币+2, 辛苦了 2014-06-11 12:25:31
|
比较死的循环代码，花了5000s的时间才算完....另外个代码可以用并行，不知道速度能提高多少,或者你编译成C算算 Clear["Global`*"]; n = 100; n6 = n1 + n2 + n3 - n4 - n5; j = 0; t1 = AbsoluteTime[]; Do[If[n <= (n1 + n2 + n2 + n3 + n4 + n5 + n6)/2 <= 3 n && n6 >= 1, j = j + 1], {n1, 1, n}, {n2, 1, n1}, {n3, 1, n - n1 + n2}, {n4, 1, n1 - n2 + n3}, {n5, 1, n - n1 + n2 - n3 + n4}]; t2 = AbsoluteTime[]; t2 - t1 Print[j] 5038.6471942 1022862022 |
4楼2014-06-09 11:08:06
aprillf
木虫 (正式写手)
- 应助: 12 (小学生)
- 金币: 3720.2
- 散金: 45
- 红花: 5
- 帖子: 562
- 在线: 248.1小时
- 虫号: 1688600
- 注册: 2012-03-13
- 性别: MM
- 专业: 地理信息系统
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
theday0323: 金币+8, ★有帮助, 摁 谢了，如果有能简化时间的算法就好了 2014-08-25 17:48:04
theday0323: 金币+20, ★★★★★最佳答案, 等了这么久，就给你了，只有你是应助贴就给你了，谢了 2014-08-28 10:26:56
theday0323: 金币+8, ★有帮助, 摁 谢了，如果有能简化时间的算法就好了 2014-08-25 17:48:04
theday0323: 金币+20, ★★★★★最佳答案, 等了这么久，就给你了，只有你是应助贴就给你了，谢了 2014-08-28 10:26:56
|
对于比较小的N可以高效计算; N=100时解的个数估计有上千万甚至更多, 所以下面的结果仅供参考了:  SetSystemOptions["ReduceOptions"->{"DiscreteSolutionBound"->10000000}]; f[n_]:=Module[{eqn},eqn=(n1+n3+n5==n2+n4+n6)&&(n1>=n2)&&(n1-n2+n3-n4>=0)&&(n1<=n)&&(n1-n2+n3<=n)&&(n1-n2+n3-n4+n5<=n)&&(2n<=n1+n2+n3+n4+n5+n6<=6n)&&(n1>0)&&(n2>0)&&(n3>0)&&(n4>0)&&(n5>0)&&(n6>0);Length@Reduce[eqn,{n1,n2,n3,n4,n5,n6},Integers]] Flatten[{#,Timing[f@#]}]&/@Range[5,30]//TableForm[#,TableHeadings->{None,{"n","timing","numbers"}}]& |
10楼2014-08-25 12:29:34
walk1997
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4677.7
- 红花: 22
- 帖子: 1066
- 在线: 798.1小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
★
jjdg: 金币+1, 感谢参与 2014-06-11 12:25:16
jjdg: 金币+1, 感谢参与 2014-06-11 12:25:16
|
这是考试题目还是实际碰到的问题呀？ 把约束条件重新排列一下 直接做个循环应该就可以 |
» 本帖已获得的红花（最新10朵）
theday03232楼2014-06-08 16:54:12
theday0323
木虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 2560
- 帖子: 149
- 在线: 214.3小时
- 虫号: 1075164
- 注册: 2010-08-15
- 性别: GG
- 专业: 光学和光电子材料
3楼2014-06-08 21:38:29
walk1997
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4677.7
- 红花: 22
- 帖子: 1066
- 在线: 798.1小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
5楼2014-06-09 11:10:12
theday0323
木虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 2560
- 帖子: 149
- 在线: 214.3小时
- 虫号: 1075164
- 注册: 2010-08-15
- 性别: GG
- 专业: 光学和光电子材料
6楼2014-06-10 21:37:38
walk1997
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4677.7
- 红花: 22
- 帖子: 1066
- 在线: 798.1小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
7楼2014-06-11 08:34:58
theday0323
木虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 2560
- 帖子: 149
- 在线: 214.3小时
- 虫号: 1075164
- 注册: 2010-08-15
- 性别: GG
- 专业: 光学和光电子材料
8楼2014-06-11 10:35:10
cooooldog
铁杆木虫 (著名写手)
ส็็็
- 应助: 237 (大学生)
- 金币: 6101.9
- 散金: 1114
- 红花: 39
- 帖子: 1380
- 在线: 553.8小时
- 虫号: 506699
- 注册: 2008-02-18
- 专业: 模式识别
9楼2014-08-24 15:02:16
