>>530
#include <stdio.h>
#include <math.h>
void GYAKU();
int P,R,II,BIT;
main(){
　　int N,IR,I,J2,K,K1,KK,L1,L;
　　float DEG,ARG,A,B,C,S,tmp;
　　printf("input N,R");
　　scanf("%d,%d",&N,&R);
　　float XR[N],XI[N];
　　XR[0]=1; XR[1]=1; XR[2]=1; XR[3]=1;
　　IR=1;
　　J2=N;
　　DEG=2*3.14159/N;
　　K=0; L1=R-1;
　　for(L=1; L<=R; L++){
　　　　J2=J2/2;
　　　　for(I=1; I<=J2; I++){
　　　　　　P=K/pow(2,L1);
　　　　　　GYAKU();
　　　　　　ARG=DEG*BIT;
　　　　　　C=cos(ARG); S=IR*sin(ARG);
　　　　　　K1=K+J2;
　　　　　　A=XR[K1]*C+XI[K1]*S;
　　　　　　B=XI[K1]*C-XR[K1]*S;
　　　　　　XR[K1]=XR[K]-A; XI[K1]=XI[K]-B;
　　　　　　XR[K]=XR[K]+A; XI[K]=XI[K]+B;
　　　　　　K=K+1;
　　　　}