d6r4f.java

针对金融方面的pso算法

import java.util.Random;
public class d6r4F
{
	static int inextp, inext;
	static double ma[] = new double[57];
	int ran3_idum;
	double  ran3(int idum)
	{
		double mk ,b, fac, mj;
		int i, ii, k, mbig, mz;
		long mseed;
		boolean iff; iff=true;
        mbig = 1000000000; mseed = 161803398; mz = 0; fac = 0.000000001;
        if ((idum < 0) || (iff == false))
		{
            iff = true;
            mj = mseed - Math.abs(idum);
            mj = mj - mbig * (int)(mj / mbig);
            ma[55] = mj;
            mk = 1;
            for (i = 1; i <= 54; i++)
			{
                ii = (21 * i) - 55 * (int)((21.0 * i)/ 55.0);
                ma[ii] = mk;
                mk = mj - mk;
                if (mk < mz) mk = mk + mbig;
                mj = ma[ii];
            }
            for (k = 1; k <= 4; k++) 
            {
                for (i = 1; i <= 55; i++)
                {
                    ma[i] = ma[i] - ma[1 + i + 30 - 55 * (int)((i + 30) / 55)];
                    if (ma[i] < mz) ma[i] = ma[i] + mbig;
				}
            }
            inext = 0;
            inextp = 31;
            idum = 1;
        }
        inext = inext + 1;
        if (inext == 56 ) inext = 1;
        inextp = inextp + 1;
        if (inextp == 56)  inextp = 1;
        mj = ma[inext] - ma[inextp];
        if (mj < mz ) mj = mj + mbig;
        ma[inext] = mj;
        b = mj * fac;
		ran3_idum = idum;
		return b;
	}

	double fnc(double x1, double x2, double x3, double x4)
	{
		double temp;
        temp = Math.sqrt(x1 * x1 + x2 * x2 + x3 * x3 + x4 * x4);
		return temp;
	}
}