摘要: 了解C++的童鞋都知道algorithm里面有個(gè)next_permutation可以求下一個(gè)排列數(shù)，通過《STL 源碼剖析》（或者自己讀代碼）可以知道其實(shí)現(xiàn)，比如：

abcd next_permutation -> abdc

那么，為什么abcd的下一個(gè)是abdc而不是acbd呢？  閱讀全文

     摘要: 怎么取得析構(gòu)函數(shù)的地址？？？  閱讀全文

     摘要: C++下的垃圾回收機(jī)制可能會(huì)在下個(gè)版本加入，我只是想通過實(shí)例，分析垃圾回收器的內(nèi)部機(jī)制，深入了解以后，在以后
的項(xiàng)目中，就可以對(duì)是否需要垃圾回收功能做出準(zhǔn)確的判斷。  閱讀全文

條件：1K內(nèi)存，1MHzCPU，每秒可以改變2^20次狀態(tài)。問：一個(gè)程序最長的運(yùn)行時(shí)間是多少？  答： 首先程序是確定性的，就說明內(nèi)存的狀態(tài)不會(huì)重復(fù)，否則就永遠(yuǎn)結(jié)束不了。從這一點(diǎn)出發(fā)，可以知道內(nèi)存的狀態(tài)共有 2^8k , 然后CPU每秒改變 2^20 個(gè)狀態(tài)，所以這臺(tái)計(jì)算機(jī)最長出現(xiàn)不重復(fù)的狀態(tài) 2^(8k-20）秒。

題目描述：設(shè)有n個(gè)正整數(shù)，將它們聯(lián)接成一排，組成一個(gè)最小的多位整數(shù)。
程序輸入：n個(gè)數(shù)程序輸出：聯(lián)接成的多位數(shù)
例如：n=2時(shí)，2個(gè)整數(shù)32，321連接成的最小整數(shù)為：32132，n=4時(shí)，4個(gè)整數(shù)55，31，312, 33 聯(lián)接成的最小整數(shù)為：312313355
[題目要求]1. 給出偽代碼即可，請(qǐng)給出對(duì)應(yīng)的文字說明，并使用上面給出的例子試驗(yàn)?zāi)愕乃惴ā?. 給出算法的時(shí)間空間復(fù)雜度。3. 證明你的算法。（非常重要）

證明：

假設(shè)： 排序后的 a0a1...an不是最小的，那么存在a0a1...ajai....an<a0a1...an,且ajai > aiaj.
那么交換ajai會(huì)使可以使a0a1...an更小，與假設(shè)a0a1...ajai....an<a0a1...an矛盾。

證明完畢。

     摘要: 這一塊主要是講move語義的，我認(rèn)為這是在C++0x中，最好的特性之一，因?yàn)樗鼛缀蹩梢酝耆该鞯奶岣咝省Ｔ赟tephan T. Lavavej這篇帖子之后，有很多評(píng)論，大體上認(rèn)為C++因?yàn)檫@些特性而變得更復(fù)雜了，而難以掌握，另初學(xué)者望而生畏。但是我認(rèn)為這是值得的，因?yàn)镃++的宗旨是：“don't pay for what you don't use 不要為你不使用的東西而付出代價(jià)&#...  閱讀全文

     摘要:   最近最大的新聞莫過于微軟發(fā)布Visual Studio2010了，對(duì)c++的支持更進(jìn)一步，其intellsence的解析也使用了和以前完全不同的方法（以前是靠編譯器，現(xiàn)在是獨(dú)立inellsence單元），番茄可能要被打入冷宮了。Stephan T. Lavavej在Visual c++ Team Blog上發(fā)布了VC10對(duì)C++0x標(biāo)準(zhǔn)的支持情況，包括：lambdas, auto，...  閱讀全文

     摘要: 類實(shí)例能做做模板參數(shù)嗎？

MyClass t;  閱讀全文

posted @ 2009-05-24 20:21 尹東斐 閱讀(2376) | 評(píng)論 (10) | 編輯 收藏

最近項(xiàng)目里總是要對(duì)很龐大的公式求導(dǎo)，很煩人，手工求導(dǎo)容易出錯(cuò)。
當(dāng)然MATLAB是個(gè)好選擇，不過當(dāng)它要錢的時(shí)候，您可能就不這么認(rèn)為了。
于是，實(shí)現(xiàn)了一個(gè)可以編譯期求導(dǎo)(不用擔(dān)心運(yùn)行時(shí)負(fù)擔(dān))的小型庫，還不完全，僅支持多項(xiàng)式，sin，cos，pow，exp，log等函數(shù)求導(dǎo)。
后期的表達(dá)式優(yōu)化做的不是很好。
下面是一些測(cè)試代碼，完整的源碼在http://www.boostpro.com/vault/index.php?action=downloadfile&filename=[math]AD.zip
實(shí)現(xiàn)部分很復(fù)雜，請(qǐng)多多指教。

只有1個(gè)函數(shù)， d(...)
支持高階，多元求導(dǎo)。

d(exp, var)(value1, value2, ...)
exp內(nèi)可以有多個(gè)變量，var表示要對(duì)其求導(dǎo)的變量，value表示求導(dǎo)以后用于計(jì)算表達(dá)式的變量的值。

比如:
d(d(x*x*x, x),x)(3.0)                           表示對(duì)x*x*x求二階導(dǎo)數(shù)在x=3.0時(shí)候的值。
d(d(x*x*y, x), y)(3.0, 4.0)                   表示d(x*x*y)/(dxdy)在x=3.0，y=4.0的值。
d(d(x*x*x, x) +d(y*x, y), y) (2.0)        則表示 (d(x*x*x)/dx + d(y*x)/dy)/dy == 0。

可以直接用cout把求導(dǎo)后的表達(dá)式輸出，不用給變量給值。
cout<<d(x*x, x)     //  結(jié)果是：2*x

這里沒有用任何迭代，是直接對(duì)表達(dá)式求導(dǎo)的。返回值是求導(dǎo)后的表達(dá)式，本質(zhì)是一個(gè)仿函數(shù)。可以用boost::function保存起來使用。
例如：

boost::function<double (double)> df = d(pow(x, const_<10>::type()), x); //df 參數(shù)為1個(gè)double，返回double

然后就可以在任何地方使用 df 了：

double res = df(3.0)   // res == pow(3, 9)


輸出結(jié)果如下：