コード例 #1
0
ファイル: coursera1.c プロジェクト: g1g19/codeforces
void countInversions(int* A,int left,int right)
{
	int mid;
	if(left < right)
	{
		mid = (left+right)/2;
		countInversions(A,left,mid);
		countInversions(A,mid+1,right);
		countSplit(A,left,mid,right);
	}
	return;
}
コード例 #2
0
//----------------------------------------------------------------------------------
void Simulation::measurePopulation(MeasurementsTable& table)
{
	if (mErrorCode!=EcNoError) return;

	checkPopSize("void Simulation::measurePopulation()");
	const unsigned int popSize = SoilMiteBase::getPopSize();

	std::vector<double> bodySizeTable(popSize);
	std::vector<double> reservesTable(popSize);
	std::vector<double> chanceToDieTable(popSize);
	std::vector<double> foodIntakeTable(popSize);
	std::vector<double> allocationTable(popSize);
	std::vector<double> fractionBiomassPerOffspringInitialTable(popSize);
	std::vector<double> fractionBiomassPerOffspringTable(popSize);

	std::vector<bool> isAdultTable(popSize);
	const double splitValue = SoilMiteBase::getBodySizeAdult();

	for (unsigned int i=0; i<popSize; ++i)
	{
		bodySizeTable[i]    = mPop[i]->getBodySize();
		reservesTable[i]    = mPop[i]->getReserves();
		foodIntakeTable[i]  = mPop[i]->getFeedRate();
		chanceToDieTable[i] = mPop[i]->getChanceToDie();
		allocationTable[i]  = mPop[i]->getAllocation();
		fractionBiomassPerOffspringInitialTable[i] = mPop[i]->getFractionBiomassPerOffspringInitial();
		fractionBiomassPerOffspringTable[i] = mPop[i]->getFractionBiomassPerOffspring();
		isAdultTable[i]     = ( bodySizeTable[i]>splitValue ? true : false);

		assert(allocationTable[i]>=0.0 && allocationTable[i]<=1.0);
		assert(fractionBiomassPerOffspringInitialTable[i]>=0.0 && fractionBiomassPerOffspringInitialTable[i]<=1.0);
		assert(fractionBiomassPerOffspringTable[i]>=0.0 && fractionBiomassPerOffspringTable[i]<=1.0);
		assert(mPop[i]->getChanceToDie()>=0.0 && mPop[i]->getChanceToDie()<=1.0);
	}

	double meanBodySize, stdDevBodySize;
	SafeMeanAndStdDev(bodySizeTable, meanBodySize, stdDevBodySize);
	table.mMeanBodySize.push_back(meanBodySize);
	table.mStdDevBodySize.push_back(stdDevBodySize);

	double meanReserves, stdDevReserves;
	SafeMeanAndStdDev(reservesTable, meanReserves, stdDevReserves);
	table.mMeanReserves.push_back(meanReserves);
	table.mStdDevReserves.push_back(stdDevReserves);

	double meanFoodIntake, stdDevFoodIntake;
	SafeMeanAndStdDev(foodIntakeTable, meanFoodIntake, stdDevFoodIntake);
	table.mMeanFoodIntake.push_back(meanFoodIntake);
	table.mStdDevFoodIntake.push_back(stdDevFoodIntake);
	
	double meanChanceToDie, stdDevChanceToDie;
	SafeMeanAndStdDev(chanceToDieTable, meanChanceToDie, stdDevChanceToDie);
	table.mMeanChanceToDie.push_back(meanChanceToDie);
	table.mStdDevChanceToDie.push_back(stdDevChanceToDie);

	//Plus
	int popSizeJuvenileInt, popSizeAdultInt;
	countSplit(splitValue, bodySizeTable, popSizeJuvenileInt, popSizeAdultInt);
	double popSizeJuvenile = static_cast<double>(popSizeJuvenileInt);
	double popSizeAdult    = static_cast<double>(popSizeAdultInt);
	table.mPopSizeJuvenile.push_back(popSizeJuvenile);
	table.mPopSizeAdult.push_back(popSizeAdult);
	
	double meanBodySizeJuvenile, meanBodySizeAdult;
	SafeMeanSplit(splitValue,bodySizeTable, meanBodySizeJuvenile, meanBodySizeAdult);
	assert(meanBodySizeJuvenile<=splitValue);
	assert(meanBodySizeAdult>=splitValue || meanBodySizeAdult==0.0);
	table.mMeanBodySizeJuvenile.push_back(meanBodySizeJuvenile);
	table.mMeanBodySizeAdult.push_back(meanBodySizeAdult);
	
	double meanReservesJuvenile, meanReservesAdult;
	SafeMeanSplit(isAdultTable,reservesTable, meanReservesJuvenile, meanReservesAdult);
	table.mMeanReservesJuvenile.push_back(meanReservesJuvenile);
	table.mMeanReservesAdult.push_back(meanReservesAdult);
	
	double meanChanceToDieJuvenile, meanChanceToDieAdult;
	SafeMeanSplit(isAdultTable,chanceToDieTable, meanChanceToDieJuvenile, meanChanceToDieAdult);
	table.mMeanChanceToDieJuvenile.push_back(meanChanceToDieJuvenile);
	table.mMeanChanceToDieAdult.push_back(meanChanceToDieAdult);
	
	double meanFoodIntakeJuvenile, meanFoodIntakeAdult;
	SafeMeanSplit(isAdultTable,foodIntakeTable, meanFoodIntakeJuvenile, meanFoodIntakeAdult);
	table.mMeanFoodIntakeJuvenile.push_back(meanFoodIntakeJuvenile);
	table.mMeanFoodIntakeAdult.push_back(meanFoodIntakeAdult);

	//Plus plus
	double meanAllocationJuvenile, meanAllocationAdult;
	SafeMeanSplit(isAdultTable,allocationTable, meanAllocationJuvenile, meanAllocationAdult);
	table.mMeanAllocation.push_back(meanAllocationAdult);

	double meanFractionBiomassPerOffspringInitialJuvenile, meanFractionBiomassPerOffspringInitialAdult;
	SafeMeanSplit(isAdultTable,fractionBiomassPerOffspringInitialTable, meanFractionBiomassPerOffspringInitialJuvenile, meanFractionBiomassPerOffspringInitialAdult);
	table.mMeanFractionBiomassPerOffspringInitial.push_back(meanFractionBiomassPerOffspringInitialAdult);

	double meanFractionBiomassPerOffspringJuvenile, meanFractionBiomassPerOffspringAdult;
	SafeMeanSplit(isAdultTable,fractionBiomassPerOffspringTable, meanFractionBiomassPerOffspringJuvenile, meanFractionBiomassPerOffspringAdult);
	table.mMeanFractionBiomassPerOffspring.push_back(meanFractionBiomassPerOffspringAdult);
}