////////////////////////////////////////////////////////////////////////////////////////////////
//	Miu.cpp																					  //
//  作用：用于計算函數Miu表																	  //
//	編寫者：王小龍（00512087）																  //
//  完成時間：2006、5、5																	  //
////////////////////////////////////////////////////////////////////////////////////////////////

#include "stdafx.h"
#include "stdlib.h"
#include "math.h"

// 得到Miu表（應該在主函數中建立）
//int *Miu =(int*)malloc(sizeof(int)*x_13);

void GetMiu(int *PList ,int *Miu, int a, int x_13)
{
	register int i,j,k;

	for(i=0; i<x_13 ; i++)
	{
		Miu[i]=1;						// 先把miu 表全設為1
	}

	for(i=0; i<a ; i++)					// 這里用所有的P表中的素數去篩，篩到的變號
	{
		for (j=PList[i]-1 ; j<x_13 ; j+= PList[i])	
		{
			Miu[j] *=-1;
		}
	}

	Miu[0] =1;								// 1的miu (這里是Miu[0]) 這個數設為0

	int x_16 = floor( sqrt(x_13) );

	for(i=0 ; PList[i] <= x_16 ; i++)		// 用p^2 去篩，篩到的置0
	{
		for (j= (k=(PList[i]*PList[i]))-1 ; j<=x_13 ; j+=k)		
		{
			Miu[j] *=0;
		}
	}
}