C++博客-5D空间-随笔分类-我的开源库

PointerPool(指针池)

今晚打老虎 — Mon, 11 Apr 2011 08:13:00 GMT

摘要: 对象池的指针版本。池中本质上是一群对象，只构造一次，返回指针。主要支持频繁重建的内存空间不连续的自建结构，如树，链表。 /**//****************************************************************** *名称:PointerPool(指针池)(模板) *版本号:0.1 *作者:赵耀(中山大学2010级)&... 阅读全文

今晚打老虎 2011-04-11 16:13 发表评论

LHRODT(非递归求解度数为2的线性齐次方程的第n项）

今晚打老虎 — Mon, 04 Apr 2011 06:40:00 GMT

　　如题，利用非递归办法解决a(n) = r1*a(n-1) + r2*a(n-2)问题。其中斐波那契数列即为r1 = r2 = a1 = a2 = 1的特例。
　　函数支持两种版本的调用，一种是完全版，一种是简洁版。完全版需要５个参数r1 r2 a1 a2 n，简洁版只需要一个参数n。具体使用方法见代码注释。

/*********************************************************************************
*名称:LHRODT.h
*版本号:0.1
*作者:赵耀(中山大学2010级)
*时间:2011.4.4
*简介:
*        linear homogeneous relation of degree two的计算函数.形如:
*            a(n) = r1*a(n-1) + r2*a(n-2)
*    这样的递归数列，只需调用函数
*            LHRODT( r1, r2, a1, a2, n )
*    即可返回数列的第a(n)项(r1,r2,a1,a2均为double类型,n为int类型,返回类型为double).
*    这里r1,r2为公式中r1,r2,而a1,a2为数列的头两项,n为第几项.函数带缓存功能,即第一次
*    调用后,下次调用可以只输入
*            LHRODT( n )
*    .
*
*未完成特性:
*    1.不含有数据检测功能,如果输入的数据无解,则会返回0.
*    2.简洁版的调用不包含条件检测机制,如果不满足条件依然会调用,但是返回0.
*    3.未含有范围检测功能，如果数据函数结果增长很快，有可能出现数据溢出而没有任何提示.
*已知bug:
*    1.无法处理无解数据的输入.
*    2.简洁版在未调用完全版或者完全版调用失败的情况下只返回0.
*    3.可能在没有任何提示的情况下出现数据溢出.
*版权信息:
*        该代码为开源代码,原作者保留其所有权.你可以拷贝,修改,使用该代码,但是请保留必
*    要的版权信息.
*********************************************************************************/

#ifndef LHRODT_H
#define LHRODT_H
#include <cmath>
using std::pow;

class fsLHRODT //fs = function support
{
public:
    fsLHRODT();
    double operator()( double, double, double, double, int );
    double operator()( int );

private:
    double x1; //x1 x2 u1 u2 det 均为计算过程的中间变量
    double x2; //result 为最后结果的临时储存
    double u1;
    double u2;
    double det;
    double result;
    bool flag; //标记完全版的函数是否被调用过
};

fsLHRODT::fsLHRODT()
    : x1( 0 ), x2( 0 ), u1( 0 ), u2( 0 ), det( 0 ), result( 0 ), flag( false )
{
}

double fsLHRODT::operator()( double r1, double r2, double a1, double a2, int n )
{
    flag = true;
    det = r1*r1 + 4*r2;

    if ( det < 0 ) //det小于0说明输入的数据不合法,不能按照公式计算,并且下次不能直接调用简洁版函数
    {
        flag = false;
        return 0;
    }
    else if ( det > 0 )
    {
        det = sqrt( det );
        x1 = ( r1 + det ) / 2;
        x2 = ( r1 - det ) / 2;
        u1 = ( a1*x2 - a2 ) / ( x1*( x2 - x1 ) );
        u2 = ( a2*x1 - a1 ) / ( x2*( x1 - x2 ) );
        result = u1*pow( x1, n ) + u2*pow( x2, n );
        return result;
    }
    else
    {
        x1 = r1 / 2;
        u2 = ( a2 - x1*a1 ) / x1*x1;
        u1 = a1 / x1 - u2;
        result = ( u1 + u2*n ) * pow( x1, n );
        return result;
    }
}

double fsLHRODT::operator()( int n )
{
    if ( flag )
    {
        if ( det < 0 )
            return 0;
        else if ( det > 0 )
        {
            result = u1*pow( x1, n ) + u2*pow( x2, n );
            return result;
        }
        else
        {
            result = ( u1 + u2*n ) * pow( x1, n );
            return result;
        }
    }
    else
        return 0;
}

fsLHRODT LHRODT;

#endif

今晚打老虎 2011-04-04 14:40 发表评论

probability(概率发生器)

今晚打老虎 — Sun, 03 Apr 2011 11:56:00 GMT

　　只有一个double参数，以百分比的形式给出发生概率。则函数每次调用按照这个概率返回true。
　　提示：函数内部会调用srand( time( 0 ) )。特别情况下需谨慎。

#ifndef PROBABILITY_H
#define PROBABILITY_H

#include <cstdlib>
#include <iostream>
#include <ctime>
using namespace std;

bool probability( double p ) //大概能精确到小数点后一位
{
    double happeness;
    static bool flag = false;

    if ( flag == false )
    {
        srand( time( NULL ) );
        flag = true;
    }

    /*
    happeness = rand() % 10001;
    happeness /= 100;
    这样的运算得到的值概率是不满足题意的。
    所以rand() % 这样的运算其实并没有实现真正的公平随机。
    */

    happeness = rand();
    happeness /= RAND_MAX;
    happeness *= 100;

    return ( happeness <= p );
}

#endif

今晚打老虎 2011-04-03 19:56 发表评论

findAnWithDegreeOfTwo(计算度数为2的齐次递归数列的第n项)

今晚打老虎 — Sun, 03 Apr 2011 11:53:00 GMT

　　递归数列a(n) = r1*a(n-1) + r2*a(n-2)
　　其中r1 r2是常数，a1 a2已知，按照顺序a1 a2 r1 r2 n输入参数，则返回这样一个递归数列的第n项值。类型均为double

#ifndef FINDAN_H
#define FINDAN_H

#include <cmath>

// a(n) = r1*a(n-1) + r2*a(n-2), give the a1, a2, r1, r2 and n
double findAnWithDegreeOfTwo( double a1, double a2, double r1, double r2, int n )
{
    double x1;
    double x2;
    double u1;
    double u2;
    double det = r1*r1 + 4*r2;
    double result;

    if ( det < 0 )
        return 0;
    else if ( det > 0 )
    {
        det = sqrt( det );
        x1 = ( r1 + det ) / 2;
        x2 = ( r1 - det ) / 2;
        u1 = ( a1*x2 - a2 ) / ( x1*( x2 - x1 ) );
        u2 = ( a2*x1 - a1 ) / ( x2*( x1 - x2 ) );
        result = u1*pow( x1, n ) + u2*pow( x2, n );
        return result;
    }
    else
    {
        x1 = r1 / 2;
        u2 = ( a2 - x1*a1 ) / x1*x1;
        u1 = a1 / x1 - u2;
        result = ( u1 + u2*n ) * pow( x1, n );
        return result;
    }
}

#endif

今晚打老虎 2011-04-03 19:53 发表评论

Zeller(计算某一天是星期几）

今晚打老虎 — Sun, 03 Apr 2011 11:47:00 GMT

　　输入三个参数分别表示年、月、日，返回一个int数字，代表这一天是星期几。计算过程采用蔡勒公式。

#ifndef ZELLER_H
#define ZELLER_H

int Zeller( int y, int m, int d )
{
    int c;

    if ( m <= 2 )
    {
        m += 12;
        y--;
    }

    c = y / 100;
    y %= 100;

    int w = ( y+(y/4)+(c/4)-2*c+(26*(m+1)/10)+d-1 ) % 7;
    if ( w <= 0 )
        w += 7;

    return w;
}

#endif

今晚打老虎 2011-04-03 19:47 发表评论