// UniformDesign.cpp : Defines the entry point for the console application.
//本程序計算的是由U18(18/4)表計算混料均勻設(shè)計5水平18次試驗的情況
#include <stdio.h>
#include <windows.h>
#include <math.h>
/////////////////////////
//#define N 18
//#define S 5
//定義U18(18/4)設(shè)計表
static unsigned int u_table[18][4] = 
{
	{1,5,9,7},
	{2,16,16,10},
	{3,13,5,17},
	{4,9,4,1},
	{5,2,14,15},
	{6,18,11,4},
	{7,7,1,13},
	{8,4,17,5},
	{9,11,6,9},
	{10,14,12,12},
	{11,1,7,3},
	{12,10,18,16},
	{13,15,2,6},
	{14,6,10,18},
	{15,12,15,2},
	{16,3,3,11},
	{17,17,8,14},
	{18,8,13,8}
};
///////////////////////////
int main(int argc, char* argv[])
{
	//中間變量
	double c_table[18][4];
	//臨時變量
	double cTmp;
	//用于遞推的變量
	double g_table[18][6];	
	//用于存儲結(jié)果
	double x_table[18][5];
	//用于循環(huán)計數(shù)
	int k,i,j;

	//輸出原表u_table
	printf("最初表u_table:\n");
	for (k=0;k<18;k++)
	{
		for (i=0;i<4;i++)
		{			
			printf("%d\t",u_table[k][i]);
		}
		printf("\n");
	}
	//計算并輸出c_tabale
	printf("中間表c_table:\n");
	for (k=0; k<18; k++ )
	{
		for (i=0; i<4; i++)
		{			
			c_table[k][i] = (double)((double)u_table[k][i] - 0.5)/18.0;
			printf("%6.5f\t",c_table[k][i]);
		}
		printf("\n");
	}	
	//遞推表
	for (k=0; k<18; k++ )
	{
		g_table[k][0] = 1.0;
		g_table[k][5] = 0.0;
		for (j=0; j<5; j++)
		{
			cTmp = pow(c_table[k][j],1.0/(double)(5-j-1));
			g_table[k][j+1] = g_table[k][j] * cTmp;
			x_table[k][j] = g_table[k][j] - g_table[k][j+1];
		}
	}
	//輸出遞推表
	printf("遞推表g_table:\n");
	for (k=0;k<18;k++)
	{
		for (i=0;i<6;i++)
		{			
			printf("%6.5f\t",g_table[k][i]);
		}
		printf("\n");
	}
	printf("最終結(jié)果x_table:\n");
	for (k=0;k<18;k++)
	{
		for (i=0;i<5;i++)
		{			
			printf("%6.5f\t",x_table[k][i]);
		}
		printf("\n");
	}
	return 0;
}