//---------------------------------------------------------------------------

#include <vcl.h>
#pragma hdrstop

#include "slnum_p.h"
//---------------------------------------------------------------------------


#pragma package(smart_init)
#pragma link "SHDocVw_OCX"
#pragma resource "*.dfm"
TForm1 *Form1;

long   a,b,c;
//---------------------------------------------------------------------------
__fastcall TForm1::TForm1(TComponent* Owner)
        : TForm(Owner)
{
}
//---------------------------------------------------------------------------
//---------------從string類中取出字符串轉換為char * -------------------------------
char * strtoa(String ss)
{
     char *sa;
     sa=new char[strlen(ss.c_str())+1];
     strcpy(sa,ss.c_str() );
     return sa;
}
/*-----------------------------------------------------------------------
版權所有：貴州省天柱縣二中  羅國文，郵編：556600，Email:xuefu998@yahoo.com.cn
本程序提供了一種關于素數判定的新算法，以這些基本函數為基礎，可幫助我們
處理關于素數的許多命題。
版權聲明：在您使用和推廣本程序時尊重作者版權，不準在不經許可的情況下發表
-------------------------------------------------------------------------*/

long Test_1(long x)            //求模2的倍數
{
     long n;
     n=x/2;
     return  n;
}

int Test_2(long x)             //判斷x模2的余數，返回值{0，1}
{
     int k;
     k=x%2;
     return k;
}

int Test_P(long n)        //P素性檢驗 當（2n+1）為素返回1，否則5返回0
{
    if(n==0)return 0;
    if((n>3)&& (n%3==1))  return 0;
     else   {
            long m,N;
            N=n/2;
            for(m=2;(m*m+m)<=N;m++)
                {
                if (Test_P(m))
                     if(n%(m*2+1)==m) return 0;
                }
            }
    return 1;
}

//數x的素性檢驗，x為素數返回1，否則為0
int Test_X(long x)
{
   if (x==1) return 0;
   if (x==2) return 1;
   if (Test_2(x)==0) return 0;
   long n;
   n=Test_1(x);
   return (Test_P(n));
}
//----------------------------------------------------------------------
void __fastcall TForm1::Button1Click(TObject *Sender)
{
  long x,A,B; int i=0;
  Form1->Label7->Caption ="";
  A=a;B=b;
  char *s,*ss;
  s=new char[65536];
  ss=new char[65536];
  Form1->RichEdit1->Lines->Clear();
//  if (A<=2) { Form1->RichEdit1->Lines->Append('2'),i=1;}
  for(x=A;x<=B;x=x+1)
          if(Test_X(x)==1)
                {++i;
                ltoa(x,s,10);
                Form1->RichEdit1->Lines->Append(s); }
  strcpy(ss,"計算完成，這個范圍內的素數有：");
  ltoa(i,s,10);
  strcat(ss,s);
  strcat(ss,"個");
  Form1->Label7->Caption =ss;
  delete s;
  delete ss;
}
//---------------------------------------------------------------------------



void __fastcall TForm1::Edit2Exit(TObject *Sender)
{
      b=atol(strtoa(Form1->Edit2->Text));
}
//---------------------------------------------------------------------------

void __fastcall TForm1::Edit1Exit(TObject *Sender)
{
      a=atol(strtoa(Form1->Edit1->Text));
}
//---------------------------------------------------------------------------
void __fastcall TForm1::Button6Click(TObject *Sender)
{
  long x,A,B; int i=0;
  A=a;B=b;
  Form1->Label7->Caption ="";
  char *s,*ss;
  s=new char[65536];
  ss=new char[65536];
  Form1->RichEdit1->Lines->Clear();
  for(x=A;x<=B;x=x+1)
          if(Test_P(x)==1)
                {++i;
                ltoa(x,s,10);
                Form1->RichEdit1->Lines->Append(s );}
  strcpy(ss,"計算完成，這個范圍內的Q數有：");
  ltoa(i,s,10);
  strcat(ss,s);
  strcat(ss,"個");
  Form1->Label7->Caption =ss;
  delete s;
  delete ss;
}
//---------------------------------------------------------------------------



void __fastcall TForm1::Edit3Exit(TObject *Sender)
{
      a=atol(strtoa(Form1->Edit3->Text));        
}
//---------------------------------------------------------------------------

void __fastcall TForm1::Edit4Exit(TObject *Sender)
{
      b=atol(strtoa(Form1->Edit4->Text));        
}
//---------------------------------------------------------------------------

void __fastcall TForm1::Edit5Exit(TObject *Sender)
{
      a=atol(strtoa(Form1->Edit5->Text));        
}
//---------------------------------------------------------------------------

void __fastcall TForm1::Edit6Exit(TObject *Sender)
{
      b=atol(strtoa(Form1->Edit6->Text));
}
//---------------------------------------------------------------------------



void __fastcall TForm1::Button3Click(TObject *Sender)
{
  long x,A,B,i=0;
  char *s,*ss;
  Form1->Label8->Caption ="";
  s=new char[65536];
  ss=new char[65536];
  A=Test_1(a);B=Test_1(b)-1;
  Form1->RichEdit2->Lines->Clear();
  for(x=A;x<=B;x=x+1)
          if(Test_P(x)*Test_P(x+1))
                {++i;
                ltoa(2*x+1,s,10);
                strcpy(ss,s);
                strcat(ss,"------>");
                ltoa(2*x+3,s,10);
                strcat(ss,s);
                Form1->RichEdit2->Lines->Append(ss );}
  strcpy(ss,"這個范圍內的孿生素數對有數有：");
  ltoa(i,s,10);
  strcat(ss,s);
  strcat(ss,"對");
  Form1->Label8->Caption =ss;
  delete s;
  delete ss;
}
//---------------------------------------------------------------------------

void __fastcall TForm1::Button5Click(TObject *Sender)
{
  Form1->RichEdit3->Lines->Clear();
  int k,flag;
  long A,i,x,r,B;
  A=a;B=b;
  Form1->Label9->Caption ="";
  char *s,*ss;
  s=new char[65536];
  ss=new char[65536];
  for (i=A-1;i<B;i++)
  {
    k=Test_2(i);
    x=Test_1(i) ;
    r=0; flag=1;
     while  (flag)
          {
          if (r>x){ Form1->RichEdit3->Lines->Append("找不到這個偶數的1+1解");return;}
          if(Test_P(x-r)*Test_P(x+r+k))
            {   ltoa(2*(x-r)+1,s,10);
                strcpy(ss,s);
                strcat(ss,"+");
                ltoa(2*(x+r+k)+1,s,10);
                strcat(ss,s);
                strcat(ss,"=");
                ltoa(2*i+2,s,10);
                strcat(ss,s);
                Form1->RichEdit3->Lines->Append(ss );
                flag=0;
            }
            r++;
          }
   }
  Form1->Label9->Caption ="計算已經完成，請查看數據……";
}
//---------------------------------------------------------------------------

void __fastcall TForm1::Button2Click(TObject *Sender)
{
  Form1->RichEdit3->Lines->Clear();
  Form1->Label9->Caption ="";
  int k;
  long i,j=0,x,r;
  char *s,*ss;
  s=new char[65536];
  ss=new char[65536];
  i=a-1;
    k=Test_2(i);
    x=Test_1(i) ;
    r=0;
     while  (r<x)
          {
          if (r>x){ Form1->RichEdit3->Lines->Append("找不到這個偶數的1+1解");return;}
          if(Test_P(x-r)*Test_P(x+r+k))
            {   j++;
                ltoa(2*(x-r)+1,s,10);
                strcpy(ss,s);
                strcat(ss,"+");
                ltoa(2*(x+r+k)+1,s,10);
                strcat(ss,s);
                strcat(ss,"=");
                ltoa(2*i+2,s,10);
                strcat(ss,s);
                Form1->RichEdit3->Lines->Append(ss );
            }
            r++;
          }
  strcpy(ss,"這個偶數2a的1+1解有：");
  ltoa(j,s,10);
  strcat(ss,s);
  strcat(ss,"個");
  Form1->Label9->Caption =ss;
  delete s;
  delete ss;
}
//---------------------------------------------------------------------------






void __fastcall TForm1::FormResize(TObject *Sender)
{

    Form1->RichEdit1->Height =Form1->PageControl1->Height-100;
    Form1->RichEdit2->Height =Form1->PageControl1->Height-100;
    Form1->RichEdit3->Height =Form1->PageControl1->Height-100;
}
//---------------------------------------------------------------------------