轨道力学中二体问题下求解兰伯特方程需要解决一个迭代问题。这是一个老外写的,有很多注释,相信大家应该能看懂,经实际检测,切实可用function
[v1,v2]=solve_lambert(r1,r2,t,GM,lw,N,branch)%This routine implements a new
al...
分类:
其他好文 时间:
2014-05-16 20:36:40
阅读次数:
1531
问题:
有N件物品和一个容量为V的背包。第i件物品的费用是c[i],价值是w[i]。求解将哪些物品装入背包可使价值总和最大。
思路:
这是最基础的背包问题,特点是:每种物品仅有一件,可以选择放或不放。用子问题定义状态:即f[i][v]表示前i件物品恰放入一个容量为v的背包可以获得的最大价值。则其状态转移方程便是:f[i][v]=max{f[i-1][v],f[i-1][v-c...
分类:
其他好文 时间:
2014-05-15 23:50:57
阅读次数:
421
连分数法解佩尔方程特解一、佩尔方程的形式:二、关于佩尔方程的特解:特解是指佩尔方程的最小整数解,容易发现当x最小的时候y也同样达到最小。在一般情况下,佩尔方程的特解是通过暴利枚举方法得到的,本文将介绍如何用应用连分数法求特解。本文将不涉及证明,只介绍方法。三、连分数法:一个实数的简单连分数表示,是指...
分类:
其他好文 时间:
2014-05-15 22:47:37
阅读次数:
537
【程序3】题目:打印出所有的"水仙花数",所谓"水仙花数"是指一个三位数,其各位数字立方和等于该数本身。例如:153是一个"水仙花数",因为153=1的三次方+5的三次方+3的三次方。1.程序分析:利用for循环控制100-999个数,每个数分解出个位,十位,百位。packagetest;
publicc..
分类:
编程语言 时间:
2014-05-15 19:05:20
阅读次数:
256
线性常微分方程解法一阶线性微分方程\[\frac{{dy}}{{dx}} + P(x)y =
Q(x)\]对应的齐次线性方程\[\frac{{dy}}{{dx}} + P(x)y = 0\]此齐次方程可以用分离变量法求得通解:\(y =
C{e^{-\int {P(x)dx} }}\)常数变易法求非...
分类:
其他好文 时间:
2014-05-14 08:52:13
阅读次数:
314
轨道力学中二体问题下求解兰伯特方程。老外写的Matlab程序,我把它转成了Fortran程序。!*****************************************************************
subroutine solve_lambert(r1,r2,tt...
分类:
其他好文 时间:
2014-05-14 08:37:56
阅读次数:
399
1 Accepted 8508K 579MS C++ 2237B/** 2
hash的强大,,还是高次方程,不过要求n不一定是素数 3 **/ 4 #include 5 #include 6 #include 7 #include 8
#inclu...
分类:
其他好文 时间:
2014-05-13 20:51:41
阅读次数:
330
1 Accepted 406MS 8576K 2379 B C++/** 2
这里加了一点限制,,大体还是一样的,, 3 4 **/ 5 #include 6 #include 7 #include 8 #include 9
#include 10...
分类:
其他好文 时间:
2014-05-13 20:15:52
阅读次数:
301
一,反演理论的初识以一个方程组的例子来引出主题:a*x+b*y=p,c*x+d*y=q.在已知a,b,c,d,p,q的前提下我们如何来求出x和y呢?在计算机解决这类问题的一个最基本方法就是高斯消元法。像这样通过已知的结果来求未知数的方法就叫做反演。用数学语言来说:为得到某个组合计数问题的解,我们首先...
分类:
其他好文 时间:
2014-05-13 10:04:15
阅读次数:
318
可以用递归简洁的写出,但是会超时。
dp嘛。这个问题需要从后往前算,最右下角的小规模是已知的,边界也很明显,是最后一行和最后一列,行走方向的限制决定了这些位置的走法是唯一的,可以先算出来。然后不断的往前推算。
用distance[i][j]保存从当前位置走到最右下角所需的最短距离,状态转移方程是从distance[i+1][j]和distance[i][j+1]中选一个小的,然后再加上自身的。...
分类:
其他好文 时间:
2014-05-13 00:13:56
阅读次数:
339