/************************************
*
* Author:LOJA
* Date: 2009/08/25
*
*
*************************************/

#include<Shellapi.h>    //WinCE

//重建目录函数(删除整个目录及其子目录子文件,然后重新创建该目录)
BOOL RebuildDirectory(const CMzString str,bool BRebulid=true)
{
    SHFILEOPSTRUCT lpFileOp;
    lpFileOp.hwnd=NULL;
    lpFileOp.wFunc=FO_Delete;
    lpFileOp.fFlags=FOF_NOCONFIRMATION|FOF_NOERRORUI|FOF_SILENT;        //不提示覆盖，不显示错误，不显示进度
    lpFileOp.pFrom=str.C_Str();

    if(SHFileOperation(&lpFileOp)==0)
    {
        if (BRebulid)
        {
            return CreateDirectory(lpFileOp.pFrom,NULL);
        }
        else
        {
            return TRUE;
        }
    }
    return FALSE;
}


//M8清理垃圾
CMzString strFolder;    //使用魅族CMzString类前请包含相关头文件

strFolder=L"Temp";
RebuildDirectory(strFolder);

strFolder=L"Windows\\Profiles";
RebuildDirectory(strFolder);

strFolder=L"Windows\\Recent";
RebuildDirectory(strFolder);

strFolder=L"Windows\\Temp";
RebuildDirectory(strFolder);

strFolder=L"Windows\\www";
RebuildDirectory(strFolder);

strFolder=L"Disk\\Temp";
RebuildDirectory(strFolder);

strFolder=L"Disk\\Application Data\\WebBrowser\\cache";
RebuildDirectory(strFolder);

strFolder=L"Disk\\Application Data\\WebBrowser\\tmp";
RebuildDirectory(strFolder);



//M8清理日志
CMzString strFolder;
strFolder=L"Disk\\Application Data\\Log";
RebuildDirectory(strFolder);

    HDC hScrDC, hMemDC;
    int width, height;

    //the pointer will save all pixel point's color value
    BYTE *lpBitmapBits = NULL;

    //creates a device context for the screen device
    hScrDC = CreateDC(_T("DISPLAY"), NULL, NULL, NULL);

    //get the screen point size
    width = GetDeviceCaps(hScrDC, HORZRES);
    height = GetDeviceCaps(hScrDC, VERTRES);

    //creates a memory device context (DC) compatible with the screen device(hScrDC)
    hMemDC = CreateCompatibleDC(hScrDC);

    //initialise the struct BITMAPINFO for the bimap infomation，
    //in order to use the function CreateDIBSection
    //on wince os, each pixel stored by 24 bits(biBitCount=24)
    //and no compressing(biCompression=0)
    BITMAPINFO RGB24BitsBITMAPINFO;
    ZeroMemory(&RGB24BitsBITMAPINFO, sizeof(BITMAPINFO));
    RGB24BitsBITMAPINFO.bmiHeader.biSize = sizeof(BITMAPINFOHEADER);
    RGB24BitsBITMAPINFO.bmiHeader.biWidth = width;
    RGB24BitsBITMAPINFO.bmiHeader.biHeight = height;
    RGB24BitsBITMAPINFO.bmiHeader.biPlanes = 1;
    RGB24BitsBITMAPINFO.bmiHeader.biBitCount = 24;

    //use the function CreateDIBSection and SelectObject
    //in order to get the bimap pointer : lpBitmapBits
    HBITMAP directBmp = CreateDIBSection(hMemDC, (BITMAPINFO*) &RGB24BitsBITMAPINFO, DIB_RGB_COLORS, (void **) &lpBitmapBits, NULL, 0);
    HGDIOBJ previousObject = SelectObject(hMemDC, directBmp);

    // copy the screen dc to the memory dc
    BitBlt(hMemDC, 0, 0, width, height, hScrDC, 0, 0, SRCCOPY);

    //if you only want to get the every pixel color value,
    //you can begin here and the following part of this function will be unuseful;
    //the following part is in order to write file;

    //bimap file header in order to write bmp file
    BITMAPFILEHEADER bmBITMAPFILEHEADER;
    ZeroMemory(&bmBITMAPFILEHEADER, sizeof(BITMAPFILEHEADER));
    bmBITMAPFILEHEADER.bfType = 0x4d42; //bmp
    bmBITMAPFILEHEADER.bfOffBits = sizeof(BITMAPFILEHEADER) + sizeof(BITMAPINFOHEADER);
    bmBITMAPFILEHEADER.bfSize = bmBITMAPFILEHEADER.bfOffBits + ((width *height) *3); ///3=(24 / 8)

    //write into file
    FILE *mStream = NULL;
    if ((mStream = _wfopen(filename, TEXT("wb"))))
    {
        //write bitmap file header
        fwrite(&bmBITMAPFILEHEADER, sizeof(BITMAPFILEHEADER), 1, mStream);
        //write bitmap info
        fwrite(&(RGB24BitsBITMAPINFO.bmiHeader), sizeof(BITMAPINFOHEADER), 1, mStream);
        //write bitmap pixels data
        fwrite(lpBitmapBits, 3 *width * height, 1, mStream);
        //close file
        fclose(mStream);
    }

    //delete
    DeleteDC(hMemDC);
    DeleteDC(hScrDC);
    DeleteObject(directBmp);
    DeleteObject(previousObject);

/************************************
* 调用范例:
* SendSMS(_T("+8613xxxxxxxxx"),_T("骚扰~"));
*
* Author:LOJA
* Version: 1.0.1.0
* Date: 2009/07/30
*
* Smartphone Platforms: Smartphone 2002 and later
* OS Versions: Windows CE 3.0 and later
*
*************************************/
#include <sms.h>
#pragma comment(lib,"sms.lib")

bool SendSMS(IN LPCTSTR lpNumber,IN LPCTSTR lpszMessage)
{

    HRESULT hRes;
    SMS_HANDLE   smsHandle=NULL;
    SMS_ADDRESS   smsaDestination;
    SMS_MESSAGE_ID   smsmidMessageID=0;
    TEXT_PROVIDER_SPECIFIC_DATA   tpsd;

    hRes=SmsOpen(SMS_MSGTYPE_TEXT,SMS_MODE_SEND,&smsHandle,NULL);
    if   (FAILED(hRes))
    {
        return   false;
    }

    //地址方式注意国内号码前加"+86"
    smsaDestination.smsatAddressType = SMSAT_INTERNATIONAL;
    _tcsncpy(smsaDestination.ptsAddress, lpNumber,SMS_MAX_ADDRESS_LENGTH);

    tpsd.dwMessageOptions  = PS_MESSAGE_OPTION_NONE;
    //tpsd.dwMessageOptions =  PS_MESSAGE_OPTION_STATUSREPORT;//表示需要状态报告
    tpsd.psMessageClass  = PS_MESSAGE_CLASS1;
    //PS_MESSAGE_CLASS0表示短信在被接收后立即显示且不存储在收件箱(称为闪信)
    //PS_MESSAGE_CLASS1表示一般的情况，被接收后存储到收件箱并发送一个确认回短信中心，发送方收到一个已被接收的状态报告。

    ZeroMemory(tpsd.pbHeaderData, sizeof(tpsd.pbHeaderData));        
    tpsd.dwHeaderDataSize = 0;
    tpsd.fMessageContainsEMSHeaders = FALSE;
    tpsd.dwProtocolID = SMS_MSGPROTOCOL_UNKNOWN;
    tpsd.psReplaceOption = PSRO_NONE;


    hRes= SmsSendMessage(smsHandle,
        NULL,  
        &smsaDestination,  
        NULL,
        (PBYTE)lpszMessage,  
        _tcslen(lpszMessage) *  sizeof(TCHAR),  
        (PBYTE)&tpsd,
        sizeof(TEXT_PROVIDER_SPECIFIC_DATA),  
        SMSDE_OPTIMAL,  
        SMS_OPTION_DELIVERY_NONE,
        &smsmidMessageID);

    SmsClose(smsHandle);

    if   (SUCCEEDED(hRes))
    {
        return true;
    }
    else
    {
        return false;
    }
}

#include <sms.h>
#pragma comment(lib,"sms.lib")
void SendSMS1(LPCTSTR lpszMessage, LPCTSTR lpszRecipient,LPCTSTR lpszSMSC,BOOL bSendConfirmation, BOOL bUseDefaultSMSC)
{

SMS_HANDLE smshHandle; //短信句柄
SMS_ADDRESS smsaDestination;//接收号码
TEXT_PROVIDER_SPECIFIC_DATA tpsd;
SMS_MESSAGE_ID smsmidMessageID;

// 打开发送句柄
// 协议字符串,       SMS_MODE_SEND:发送模式 psmshHandle：返回的句柄 phMessageAvailableEvent:通知的事件
if(FAILED(SmsOpen(SMS_MSGTYPE_TEXT, SMS_MODE_SEND, &smshHandle, NULL)))
{
  return;
}


// Create the destination address
//SMSAT_NATIONAL：国内电话 SMSAT_INTERNATIONAL：国际电话
smsaDestination.smsatAddressType = SMSAT_UNKNOWN;
_tcsncpy(smsaDestination.ptsAddress, lpszRecipient, SMS_MAX_ADDRESS_LENGTH);

// Set up provider specific data
memset(&tpsd, 0, sizeof(tpsd));
//tpsd.dwMessageOptions =  PS_MESSAGE_OPTION_NONE;
tpsd.dwMessageOptions =  PS_MESSAGE_OPTION_STATUSREPORT;//表示需要状态报告
tpsd.psMessageClass = PS_MESSAGE_CLASSUNSPECIFIED;
//PS_MESSAGE_CLASS0表示短信在被接收后立即显示且不存储在收件箱(称为闪信)
//PS_MESSAGE_CLASS1表示一般的情况，被接收后存储到收件箱并发送一个确认回短信中心，发送方收到一个已被接收的状态报告。
tpsd.psReplaceOption = PSRO_NONE;
//手动清空头信息
ZeroMemory(tpsd.pbHeaderData, sizeof(tpsd.pbHeaderData));        
tpsd.dwHeaderDataSize = 0;
tpsd.fMessageContainsEMSHeaders = FALSE;
tpsd.dwProtocolID = SMS_MSGPROTOCOL_UNKNOWN;

// Send the message, indicating success or failure
/*
*smshHandle: SmsOpen打开的句柄
􀁺 psmsaSMSCAddress: SMSC中心
􀁺 *psmsaDestinationAddress:发送的目的地址
􀁺 pstValidityPeriod :发送时间的有效期
􀁺 *pbData: 信息的内容部分
􀁺 *dwDataSize: 信息的内容部分的长度
􀁺 *pbProviderSpecificData：针对运营商的附加数据
􀁺 *dwProviderSpecificDataSize ：附加数据的长度
􀁺 smsdeDataEncoding ：短信编码[SMSDE_OPTIMAL是其
推荐值]表示由短信中心决定字符的格式
􀁺 dwOptions :其他选项
􀁺 psmsmidMessageID: 用于得到系统回执的信息
SmsGetMessageStatus
*/


//SMSDE_GSM使用7-bit编码
//SMSDE_UCS2使用16-bit (Unicode)编码
if(SUCCEEDED(SmsSendMessage(smshHandle, NULL,
  &smsaDestination, NULL, (PBYTE) lpszMessage,
  _tcslen(lpszMessage) * sizeof(TCHAR), (PBYTE) &tpsd,
  sizeof(TEXT_PROVIDER_SPECIFIC_DATA), SMSDE_UCS2,
  SMS_OPTION_DELIVERY_NONE, &smsmidMessageID)))
{
  //Tstr=L"发送完成";


  //发送完成
  MessageBox(    NULL,L"SEND OK",L"MSG",MB_OK);
}
else
{

  //发送失败
  //Tstr=L"发送失败";
  MessageBox(NULL,L"发送失败",L"MSG",MB_OK);
}

// clean up
VERIFY(SUCCEEDED(SmsClose(smshHandle)));

}

注意：LPCTSTR 和 PCWSTR是宽字符，所以用16-bit (Unicode)编码。网上的HelloSMS的程序不够严谨，如果照用可能造成短信中心网关把短信给吞掉的情况，程序虽然调用成功，但是就是目标号码收不到短信。函数参数中的后三个参数可以不用或设默认值都可以。

使用Windows CE的SMS API函数SmsSendMessage可以发送7-bit编码和16-bit编码的短消息（未测试8-bit编码）。单条短消息的最大字符数分别是160和70。发送短消息全部为ASCII字符，SMS_DATA_ENCODING参数选择SMSDE_GSM或SMSDE_OPTIMAL都可以。发送短消息不全为ASCII字符，SMS_DATA_ENCODING参数选择SMSDE_UCS2或SMSDE_OPTIMAL都可以。建议选择SMSDE_OPTIMAL。

发送短消息全部为ASCII字符，且字符个数超过160个时，系统会自动将短消息进行分割，但并不是按每条160进行分割，而是按每条153。发送短消息不全为ASCII字符，且字符个数超过70个时，系统会自动将短消息进行分割，但并不是按每条70进行分割，而是按每条67。

发送短消息的字符数超过160或70时，应该自行进行分割，而不应交由系统进行分割，这样可以最大限度的利用单条短消息的最大字符数分别是160和70的条件.

上面的代码是通过MAPI来实现，任意一部有SIM卡的PDA/PPC上都已运行。


下面的代码是通过串口和AT指令来发送短信，需要有Gms Modem 才可以运行。
HANDLE hSerial;
int res;
unsigned long bytes;

//COM1 COM2 COM4
hSerial = CreateFile(_T("COM1:"), GENERIC_READ|GENERIC_WRITE,FILE_SHARE_READ,NULL,OPEN_EXISTING,0,NULL);



if(hSerial == NULL)
{
  ///L"串口打开失败";
  return;
}
///配置串口
DCB  PortDCB;  
PortDCB.DCBlength = sizeof(DCB);
// 默认串口参数
GetCommState(hSerial, &PortDCB);
PortDCB.BaudRate =115200 ; // baud
PortDCB.ByteSize = 8;     // Number of bits/byte, 4-8
PortDCB.Parity = NOPARITY;
PortDCB.StopBits = ONESTOPBIT;
PortDCB.fBinary=TRUE;
if (! SetCommState(hSerial, &PortDCB))
{
  ///L"配置串口失败";
  return;
}
////配置超时值
COMMTIMEOUTS  CommTimeouts;
GetCommTimeouts(hSerial, &CommTimeouts);
CommTimeouts.ReadIntervalTimeout = 100;
CommTimeouts.ReadTotalTimeoutMultiplier = 100;
CommTimeouts.ReadTotalTimeoutConstant = 100;  
CommTimeouts.WriteTotalTimeoutMultiplier =100;
CommTimeouts.WriteTotalTimeoutConstant = 100;  
if (!SetCommTimeouts(hSerial, &CommTimeouts))
{
  
  return;
}                

//设置串口“感兴趣”的事件
SetCommMask(hSerial, EV_RXCHAR);

//设置输入输出缓冲区
SetupComm(hSerial, 1024, 1024);

//初始化缓冲区中的信息
PurgeComm(hSerial, PURGE_TXCLEAR|PURGE_RXCLEAR);

//-------------------------------------串口设置完毕----------------------------------------


UpdateData(TRUE);

char str_last[2000];
int num_Length=0;
int context_length=0;
TCHAR tmp[2];
int i=0;
char addr[100]="0891683108100005F0";//短信中心号码
char phone[100]="11000D91";
char msg[1000];
// unsigned char *msg_tmp;
unsigned char msg_tmp[1000];


char str_tmp[100]="";
char str_tmp2[50]="000800";

num_Length=m_num.GetLength();
context_length=m_context.GetLength();

m_num=L"86"+m_num;
if(num_Length%2){m_num=m_num+L"F";}


for (i=0; i <(m_num.GetLength()); )  //处理 目标电话号码 奇偶换位
{
  phone[i+9] = m_num.GetAt(i);
  phone[i+8] = m_num.GetAt(i+1);
  i+=2;
}

WideCharToMultiByte(CP_ACP, 0, m_context, m_context.GetLength(), str_tmp, 160, NULL, NULL);

// msg_tmp=(unsigned char *)malloc(sizeof(unsigned char ));
gsmEncodeUcs2(str_tmp,msg_tmp,context_length);//UCS2编码

for(i=0;i<context_length*2;i++)//将编后的码转换为16进制以进行传输 前两位预留
{
  wsprintf(tmp, L"%02X", msg_tmp[i]);
  msg[2*i+2]=tmp[0];
  msg[2*i+3]=tmp[1];

}


wsprintf(tmp, L"%02X", context_length*2);//前两位存储msg的长度
msg[0]=tmp[0];
msg[1]=tmp[1];

//-------------------------------addr phone msg 全部就绪----------------------------
// 最终信息为 str_last=addr+phone+000800+msg

for(i=0;i<18 ;i++)
{
  str_last[i]=addr[i];
}

for(i=0;i<22 ;i++)
{
  str_last[i+18]=phone[i];
}


for(i=0;i<6 ;i++)
{
  str_last[i+40]=str_tmp2[i];
}

for(i=0;i<context_length*4+2 ;i++)
{
  str_last[i+46]=msg[i];
}

str_last[i+46]='\32';


//--------------------------开始发送--------------------------------------------

char str1[100]="AT+CMGS=25\r";
wsprintf(tmp, L"%02d",(context_length*4+30)/2);//计算需要发送的字节

str1[8]=tmp[0];
str1[9]=tmp[1];

char str2[100]="";
res=WriteFile(hSerial,str1,100,&bytes,NULL);
res=ReadFile(hSerial,str2,100,&bytes,NULL);



char str3[100]="";
res=WriteFile(hSerial,str_last,1000,&bytes,NULL);
res=ReadFile(hSerial,str3,100,&bytes,NULL);
//----------------------------发送结束-------------------------------------------------------

// free(msg_tmp);
CloseHandle(hSerial);
}

BOOL CSendMessage2::OnInitDialog()
{
CDialog::OnInitDialog();

// TODO: Add extra initialization here
// m_num=L"15010657623";
m_num=L"13811916883";

// m_context=L"1192323 392486";
m_context=L"工作愉快！";
UpdateData(false);  
return TRUE;

网上的HelloSMS的发送短信代码，有不完善的地方。用它来发短信可能会有问题，下面这段代码是我对它的一个改进，基本可以运行在实际的项目中：
////////////////////////////////////////////////////////////////////////////
//                         发送短信
//lpszSMSC:短信中心号码  lpszRecipient:接收号码(目标号码)
//lpszMessage:发送短信内容
//bUseDefaultSMSC:目标地址  
//bSendConfirmation:消息发送选项
////////////////////////////////////////////////////////////////////////////
void CCellGuardApp::SendSMS( LPCTSTR lpszMessage, LPCTSTR lpszRecipient,LPCTSTR lpszSMSC,BOOL bSendConfirmation, BOOL bUseDefaultSMSC)
{
SMS_HANDLE smshHandle; //短信句柄
SMS_ADDRESS smsaDestination;//接收号码
TEXT_PROVIDER_SPECIFIC_DATA tpsd;
SMS_MESSAGE_ID smsmidMessageID;

// 打开发送句柄
// 协议字符串,       SMS_MODE_SEND:发送模式 psmshHandle：返回的句柄 phMessageAvailableEvent:通知的事件
if(FAILED(SmsOpen(SMS_MSGTYPE_TEXT, SMS_MODE_SEND, &smshHandle, NULL)))
{
  return;
}

// Create the destination address
smsaDestination.smsatAddressType = SMSAT_UNKNOWN;
_tcsncpy(smsaDestination.ptsAddress, lpszRecipient, SMS_MAX_ADDRESS_LENGTH);

// Set up provider specific data
memset(&tpsd, 0, sizeof(tpsd));
//tpsd.dwMessageOptions =  PS_MESSAGE_OPTION_NONE;
tpsd.dwMessageOptions =  PS_MESSAGE_OPTION_STATUSREPORT;
tpsd.psMessageClass = PS_MESSAGE_CLASSUNSPECIFIED;
tpsd.psReplaceOption = PSRO_NONE;
//tpsd.dwHeaderDataSize = 0;
ZeroMemory(tpsd.pbHeaderData, sizeof(tpsd.pbHeaderData));        
tpsd.dwHeaderDataSize = 0;
tpsd.fMessageContainsEMSHeaders = FALSE;
tpsd.dwProtocolID = SMS_MSGPROTOCOL_UNKNOWN;

// Send the message, indicating success or failure
if(SUCCEEDED(SmsSendMessage(smshHandle, NULL,
  &smsaDestination, NULL, (PBYTE) lpszMessage,
  _tcslen(lpszMessage) * sizeof(TCHAR), (PBYTE) &tpsd,
  sizeof(TEXT_PROVIDER_SPECIFIC_DATA), SMSDE_UCS2,
  SMS_OPTION_DELIVERY_NONE, &smsmidMessageID)))
{
  //Tstr=L"发送完成";
  //发送完成
  ;
}

// clean up
VERIFY(SUCCEEDED(SmsClose(smshHandle)));
}

#include<iostream>
#include<cmath>
using namespace std;
double hls(double * p,int n);

int main()
{
    int n,i,j;
    cout<<"----------------------N阶行列式计算器(LOJA Studio)----------------------"
        <<endl
        <<">>请输入要计算的行列式阶数：";
    cin>>n;
    if(n<1)
    {
        cout<<"输入阶数错误！"<<endl;cin.get();return 0;
    }
    cout<<endl
        <<">>请逐个输入行列式中的元素(实数)"<<endl
        <<"  a(行,列)"<<endl;
    double * p=new double[n*n];
    for (i=1;i<=n;++i)
    {
        for(j=1;j<=n;++j)
        {
            cout<<"  a("<<i<<","<<j<<")=";
            cin>>*(p+(i-1)*n+(j-1));

        }
    }
    cout<<">>正在计算……"<<endl;
    cout<<">>此"<<n<<"阶行列式计算结果为："<<hls(p,n)<<endl<<endl
        <<"                    欢迎访问http://www.loja.cn ^_^"<<endl;
    delete [] p;
    cin.get();
    cin.get();
    cin.get();
    return 0;
}


double hls(double * p,int n)                               //行列式计算函数，参数一为存储行列式的数组，参数二为阶数
{
    if(n<1)
    {
        return -1;
    }

    if(n==1)                                               //若行列式阶数为1，则返回其内的唯一元素值
    {
        return *p;
    }
    //分割符-----LOJA Studio-----请保留//
    double *curr=new double[(n-1)*(n-1)];                  //定义临时指针变量指向当前行列式余子式数组
    double ans=0;
    int i,r,c,cc;
    //分割符-----LOJA Studio-----请保留//
    for(i=1;i<=n;++i)                                      //求当前第一行第i列余子式(按第一行展开)
    {
        for(r=2;r<=n;++r)
        {
            cc=-1;
            for(c=1;c<=n;++c)
            {
                if(c!=i)
                {
                    *(curr+(r-2)*(n-1)+(++cc))=*(p+(r-1)*n+(c-1));
                }
            }
        }
        ans+=*(p+i-1)*pow(-1,1+i)*hls(curr,n-1);           //计算求代数余子式之和，使用递归调用计算其余子式的值
    }
    //分割符-----LOJA Studio-----请保留//
    delete [] curr;                                        //释放存储临时变量所占空间
    return ans;
}

#include<iostream>
using namespace std;
const int H=10 ;
double jc(int x);
int main()
{
    int m,n,i;
    cout<<"杨辉三角"<<endl;
    for (n=0;n<H;++n)
    {
        for (i=0;i<=H-n;++i)
        {
        cout<<" ";
        }
        for (m=0;m<=n;++m)
        {
            cout<<jc(n)/jc(n-m)/jc(m)<<" ";
        }
        cout<<endl;
    }
    cin.get();
    return 0;
}

double jc(int x)
{
    int i;
    double s=1;
    for (i=1;i<=x;++i)
    {
        s*=i;
    }
    return s;
}

#include <iostream>
using namespace std;
int main()
{
    int a[3][3],i,j;
    for(i=0;i<=2;i++)
    {
        for(j=0;j<=2;j++)
        {
            cout<<"Input  a["<<i<<"]["<<j<<"]=";
            cin>>a[i][j];
        }
    }
    system("cls");
    cout<<"原矩阵:"<<endl;
    for(i=0;i<=2;i++)
    {
        for(j=0;j<=2;j++)
        {
            cout<<a[i][j]<<" ";
        }
        cout<<endl;
    }
    cout<<endl<<"反矩阵:"<<endl;
    for(j=0;j<=2;j++)
    {
        for(i=0;i<=2;i++)
        {
            cout<<a[i][j]<<" ";
        }
        cout<<endl;
    }
    return 0;
}

10 CLEAR : CLS : GOSUB 2000: GOTO 200
50 CLEAR : CLS
60 PRINT "  是否欣赏歌曲?(Y/N)";
70 INPUT A$
80 IF A$ = "Y" or A$ = "y" THEN 100
90 IF A$ = "N" or A$ = "n" THEN 183
95 BEEP: CLS : GOTO 50
100 CLEAR : CLS : PRINT "请欣赏《长江之歌》 ..."
110 PLAY "MNT115": FOR I = 1 TO 2
120 PLAY "O4E8F8 GAE4C4O3G4 A2. O4D8E8F4O3G8. G16O4F4. E8 E2. E8F8G4E4C4O3G4"
130 PLAY "A2.O4D8E8 F4O3G8. G16O4E4. D8 C2. P4 O3B4. O4C8D4O3G8. G16G4B8O4C8D2C4.D8E 8O3G8O4C8D8 E1 D4.E8F4O3A8.A16"
140 PLAY "O4E4. D8D2 F4. E8D8. D16D8A8 G2. E8F8 G4E4C4O3G4A2. O4D8E8F4G8.G16F4.E8E2. E8F8"
150 PLAY "G4E4C4O3G4 A2.O4D8E8 O3G4O4F8.F16E4D4"
160 IF I = 1 THEN PLAY "O4C2."
170 NEXT I
180 PLAY "O4C2.E8F8G2.E8F8 F2.E4 O4E1"
183 CLS : P = 0: A$(1) = "感谢您的使用, 再见 ! ": A$(2) = "                      "
184 FOR J = 1 TO 4: P = P + 1: IF P > 2 THEN P = 1
185 FOR I = 1 TO LEN(A$(P)) STEP 2
186 LOCATE 12, 30: PRINT LEFT$(A$(P), I)
187 FOR K = 1 TO 25000: NEXT K: NEXT I
188 FOR L = 1 TO 7000: NEXT L: NEXT J
189 END
200 CLEAR : CLS
210 PRINT "", "", "================================<理科辅助工具>=================================="
220 PRINT "", "", "", "                               1.解一元二次方程  "
230 PRINT "                               2.求最简平方根    "
340 PRINT "                               3.求最简分数      "
350 PRINT "                               4.勾股定理        "
360 PRINT "                               5.公倍,公约求解   "
370 PRINT "                               6.质、素数判断    "
380 PRINT "                               7.求平面内两点距离"
382 PRINT "                               8.用素数表示偶数  "
383 PRINT "                               9.有关圆的计算    "
385 PRINT "                               0.退出本程序      "
386 PRINT "                                                              By LOJA"
387 PRINT "                                                       (C) LOJA Studio 2002 "
388 PRINT "                                                        http://www.loja.cn"
390 LOCATE 18, 5: PRINT "请选择:(0-9)"; : INPUT AJ
400 IF AJ < 0 or AJ > 9 THEN BEEP: LOCATE 18, 5: PRINT "选择错误! 按ENTER键后请重新选择~"; : INPUT A: GOTO 200
410 IF AJ = 1 THEN 500
420 IF AJ = 2 THEN 650
430 IF AJ = 3 THEN 800
440 IF AJ = 4 THEN 1400
450 IF AJ = 5 THEN 1200
460 IF AJ = 6 THEN 1100
465 IF AJ = 7 THEN 1600
468 IF AJ = 8 THEN 1700
469 IF AJ = 9 THEN 1900
470 IF AJ = 0 THEN 50
480 END
500 CLS : CLEAR
510 PRINT "================================<解一元二次方程>================================"
520 PRINT "                       aX*X+bX+C=0"
530 PRINT "                        a="; : INPUT BA
540 PRINT "                        b="; : INPUT BB
550 PRINT "                        c="; : INPUT BC
560 BD = BB * BB - 4 * BA * BC
570 IF BD < 0 or BA = 0 THEN PRINT "     方程无实根! 直接按ENTER键继续,输入88返回上级菜单"; : BEEP: INPUT A: IF A = 88 THEN 200 ELSE GOTO 500
580 X(1) = (-BB + SQR(BD)) / (2 * BA): X(2) = (-BB - SQR(BD)) / (2 * BA)
590 IF X(1) <> X(2) THEN PRINT "方程有两不同实数根", "        X1= "; X(1), "        X2= "; X(2) ELSE PRINT "方程有两相同实数根", "       X1=X2= "; X(1)
600 GOSUB 1000: GOTO 500: END
650 CLS : CLEAR
660 PRINT "================================<求最简平方根>=================================="
670 PRINT "                       X="; CHR$(253); " "; CHR$(251); " y"
680 PRINT "                     y="; : INPUT BY: IF BY < 0 THEN CLS : PRINT "被开方数不能小于为零! 直接按ENTER键重新输入,输入88返回上级菜单~"; : INPUT A: IF A = 88 THEN 200 ELSE BEEP: BEEP: GOTO 650
685 IF BY = 0 THEN 750
690 BKS = BY: BS = 1
700 IF BKS = 0 or BKS = 1 THEN BS = BKS: A$ = STR$(BKS): GOTO 750
710 FOR I = 2 TO INT(SQR(BKS)): IF BKS / (I * I) = INT(BKS / (I * I)) THEN BS = BS * I: BKS = BKS / (I * I)
720 NEXT I
730 A$ = "": IF BS <> 1 THEN A$ = A$ + STR$(BS)
740 IF BKS <> 1 THEN A$ = A$ + CHR$(251) + STR$(BKS)
750 PRINT "                 y="; A$
760 GOSUB 1000: GOTO 650
800 CLS : CLEAR
810 PRINT "=================================<求最简分数>==================================="
820 PRINT "                       X=a/b"
830 PRINT "                     a="; : INPUT BA: IF BA = 0 THEN CLS : PRINT "分子不能为零! 直接按ENTER键重新输入,输入88返回上级菜单~"; : INPUT A: IF A = 88 THEN 200 ELSE BEEP: BEEP: GOTO 800
840 PRINT "                     b="; : INPUT BB: IF BB = 0 THEN CLS : PRINT "分母不能为零! 直接按ENTER键重新输入,输入88返回上级菜单~"; : INPUT A: IF A = 88 THEN 200 ELSE BEEP: BEEP: GOTO 840
850 BC = INT(BA / BB): BA = INT(BA - BC * BB + .5): YS = 2: YS(1) = BA: YS(2) = BB
860 IF BA > 0 THEN 900 ELSE GOTO 950
900 GYS = 1: GBS = 1
905 FOR I = 2 TO YS(YS): KG = 0: FOR K = 1 TO YS
910 IF YS(K) / I = INT(YS(K) / I) THEN KG = KG + 1: YS(K) = YS(K) / I
920 NEXT K: IF KG = YS THEN GYS = GYS * I
930 IF KG > 0 THEN GBS = GBS * I: GOTO 905
940 NEXT I
950 PRINT "                   "; BA; "/"; BB; "=";
960 IF BC <> 0 THEN PRINT BC
970 IF BA <> 0 THEN PRINT "("; BA / GYS; "/"; BB / GYS; ")"
980 GOSUB 1000: GOTO 800
1000 PRINT "计算完毕,是否继续使用?(Y/N)";
1010 INPUT A$
1020 IF A$ = "Y" or A$ = "y" THEN RETURN
1030 IF A$ = "N" or A$ = "n" THEN 200
1040 BEEP: CLS : GOTO 1000: END
1100 CLS : CLEAR
1110 PRINT "=================================<质\素数判断>=================================="
1120 PRINT "                    请输入待测数值>"; : INPUT BA: LJ = BA: IF BA <> INT(BA) or BA < 2 THEN CLS : PRINT "输入的数值错误! 直接按ENTER键重新输入,输入88返回上级菜单~"; : INPUT A: IF A = 88 THEN 200 ELSE BEEP: BEEP: GOTO 1100
1130 A$ = "": IF BA = 2 THEN A$ = "2": GOTO 1160
1140 FOR I = 2 TO INT(SQR(BA)): IF BA / I = INT(BA / I) THEN A$ = A$ + STR$(I) + "*": BA = BA / I: GOTO 1140
1150 NEXT I: IF BA <> 1 THEN A$ = A$ + STR$(BA) ELSE A$ = LEFT$(A$, LEN(A$) - 1)
1160 IF A$ = STR$(BA) THEN PRINT "              "; LJ; "为质数": GOSUB 1000: GOTO 1100
1170 PRINT "          "; LJ; "="; A$: GOSUB 1000: GOTO 1100
1200 CLS : CLEAR : YS = 1
1210 PRINT "================================<公倍,公约求解>================================="
1230 PRINT "                 已知数"; RIGHT$("0" + STR$(YS), 2); "="; : INPUT YS(YS): IF YS(YS) <= 0 THEN CLS : PRINT "已知数必须不大于零!  直接按ENTER键重新输入,输入88返回上级菜单~"; : INPUT A: IF A = 88 THEN 200 ELSE CLS : BEEP: GOTO 1210
1240 IF YS(YS) < YS(YS - 1) THEN SWAP YS(YS), YS(YS - 1)
1250 PRINT "                 继续输入?(Y/N)"; : INPUT OO$: IF OO$ = "Y" or OO$ = "y" THEN YS = YS + 1: GOTO 1230
1260 IF OO$ = "n" or OO$ = "N" THEN 1270 ELSE CLS : PRINT "选择错误! 直接按ENTER键重新选择,输入88返回上级菜单~"; : INPUT A: IF A = 88 THEN 200 ELSE CLS : BEEP: GOTO 1210
1270 GYS = 1: GBS = 1
1280 FOR I = 2 TO YS(YS): KG = 0: FOR K = 1 TO YS
1290 IF YS(K) / I = INT(YS(K) / I) THEN KG = KG + 1: YS(K) = YS(K) / I
1300 NEXT K: IF GK = YS THEN GYS = GYS * I
1310 IF KG > 0 THEN GBS = GBS * I: GOTO 1280
1320 NEXT I
1330 FOR I = 1 TO YS: GBS = GBS * YS(I): NEXT I
1340 PRINT "                  公约数:"; GYS; "    公倍数:"; GBS
1350 GOSUB 1000: GOTO 1200
1400 CLS : CLEAR
1410 PRINT "==================================<勾股定理>===================================="
1420 PRINT "                        请输入已知边长"
1430 PRINT "                第一边:"; : INPUT BA: IF BA <= 0 THEN CLS : PRINT "输入的数不能小于等于零! 直接按ENTER键重新输入,输入88返回上级菜单~"; : INPUT A: IF A = 88 THEN 200 ELSE BEEP: CLS : GOTO 1400
1440 PRINT "                第二边:"; : INPUT BB: IF BB <= 0 THEN CLS : PRINT "输入的数不能小于等于零! 直接按ENTER键重新输入,输入88返回上级菜单~"; : INPUT A: IF A = 88 THEN 200 ELSE BEEP: GOTO 1440
1450 PRINT "                      直角边:"; : BKS = ABS(BA * BA - BB * BB): GOSUB 1456: IF A$ = "0" THEN PRINT " 不存在!" ELSE PRINT A$
1453 GOTO 1490
1456 BS = 1
1458 IF BKS = 0 or BKS = 1 THEN BS = BKS: A$ = STR$(BKS): GOTO 1475
1460 FOR I = 2 TO INT(SQR(BKS)): IF BKS / (I * I) = INT(BKS / (I * I)) THEN BS = BS * I: BKS = BKS / (I * I): GOTO 1460
1465 NEXT I
1467 A$ = "": IF BS <> 1 THEN A$ = A$ + STR$(BS)
1470 IF BKS <> 1 THEN A$ = A$ + CHR$(251) + STR$(BKS)
1475 RETURN
1490 PRINT "                      斜边:"; :  BKS = BA * BA + BB * BB: GOSUB 1456: IF A$ = "0" THEN PRINT " 不存在!" ELSE PRINT A$
1500 GOSUB 1000: GOTO 1400
1600 CLEAR : CLS
1610 PRINT "==================================<点距公式>===================================="
1620 PRINT "                      请输入一点的坐标(X1,Y1)  X1="; : INPUT XA
1625 PRINT "                                               Y1="; : INPUT YA
1630 PRINT "                      请输入另一点的坐标(X2,Y2) X2="; : INPUT XB
1635 PRINT "                                                Y2="; : INPUT YB
1640 S = SQR((XB - XA) ^ 2 + (YB - YA) ^ 2)
1650 PRINT "                          两点之间距离为 S="; S
1660 GOSUB 1000: GOTO 1600
1700 CLEAR : CLS
1710 PRINT "================================<用素数表示偶数>================================"
1720 PRINT "                      请输入偶数:"; : INPUT N: IF N <= 0 THEN CLS : PRINT "输入的偶数不能小于等于零! 直接按ENTER键重新输入,输入88返回上级菜单~"; : INPUT A: IF A = 88 THEN 200 ELSE BEEP: CLS : GOTO 1700
1730 CN$ = STR$(N / 2)
1740 FOR I = 1 TO LEN(CN$)
1750 CY$ = MID$(CN$, I, 1): IF CY$ = "." THEN PRINT "     输入错误,请输入偶数!  直接按ENTER键重新输入,输入88返回上级菜单~"; : INPUT A: IF A = 88 THEN 200 ELSE BEEP: CLS : GOTO 1700
1760 NEXT I
1770 FOR A = 6 TO N STEP 2
1780 FOR B = 3 TO A / 2 STEP 2
1790 FOR C = 2 TO SQR(B)
1800 IF B / C = INT(B / C) THEN 1880
1805 LOCATE 3, 8: PRINT INT(A / N * 100); "%"
1810 NEXT C
1820 D = A - B
1830 FOR E = 2 TO SQR(D)
1840 IF D / E = INT(D / E) THEN 1880
1850 NEXT E
1860 IF A = N THEN LOCATE 3, 5: PRINT "          "; A; "="; B; "+"; D
1870 GOTO 1890
1880 NEXT B
1890 NEXT A
1895 GOSUB 1000: GOTO 1700
1900 CLS : CLEAR
1910 PRINT "=================================<有关圆的计算>================================="
1920 PRINT "                请输入圆的半径:"; : INPUT R: IF R <= 0 THEN CLS : PRINT "输入的半径不能小于等于零! 直接按ENTER键重新输入,输入88返回上级菜单~"; : INPUT A: IF A = 88 THEN 200 ELSE BEEP: CLS : GOTO 1900
1930 L = 2 * 3.1415926# * R: S = 3.1415926# * R * R
1940 PRINT "            圆的面积为:"; S; "    圆的周长为:"; L
1950 GOSUB 1000: GOTO 1900
2000 CLS : P = 0: A$(1) = "欢迎您使用<理科辅助工具> ": A$(2) = "                          "
2010 FOR J = 1 TO 4: P = P + 1: IF P > 2 THEN P = 1
2005 FOR I = 1 TO LEN(A$(P)) STEP 2
2020 LOCATE 12, 26: PRINT LEFT$(A$(P), I)
2030 FOR K = 1 TO 15000: NEXT K: NEXT I
2040 FOR L = 1 TO 7000: NEXT L: NEXT J
2050 LOCATE 11, 30: PRINT "Now Loading ...": LOCATE 12, 30: PRINT "Please wait ...": FOR I = 1 TO 499999: NEXT I
2060 RETURN

0 CLEAR:DIM S$(100):DIM YS(100)
10 S$(1)=" 解一元二次方程":S$(2)=" 求最简平方根":S$(3)=" 求最简分数":S$(4)=" 勾股定理"
20 WR$="By LOJA":DIM BA(100)
30 S$(5)=" 公倍，公约求解":S$(6)=" 点到直线距离":S$(6)=" 质、素数判断"
1000 A=1
1010 B=A+4
1015 FOR I=A TO B:IF S$(I)="" THEN B=I:1020
1016 NEXT
1020 GRAPH:PRINT S$(A),S$(A+1),S$(A+2),S$(A+3),S$(A+4);
1030 LOCATE 1,1:PRINT " "," "," "," "," ";
1052 IF A>1 THEN LOCATE 1,19:PRINT CHR$(28);
1053 IF B/5=INT(B/5) AND S$(B+1)<>"" THEN LOCATE 5,19:PRINT CHR$(29);
1055 IF S$(B)="" THEN B=B-1:1055
1056 IF R>B-A THEN R=B-A
1058 LOCATE R+1,1:PRINT CHR$(16);
1060 M$=INKEY$:M=ASC(M$)
1065 IF M=27 THEN 1550
1070 IF M=20 AND R>0 THEN R=R-1:1030
1075 IF M=20 AND R=0 AND A>1 THEN A=A-1:1010
1080 IF M=21 AND R=4 AND S$(B+1)<>"" THEN A=A+1:1010
1085 IF M=21 AND R<4 AND R<I-A-1 THEN R=R+1:1030
1090 IF M=19 THEN IF A>1 THEN A=A-5:IF A<1 THEN 1000 ELSE 1010
1100 IF M=14 THEN IF S$(A+5)<>"" THEN A=A+5:1010
1110 IF M=13 THEN 1140
1120 BEEP: GOTO 1060
1130 LOCATE 5,19:PRINT "  ";RIGHT$("0"+STR$(I),2);".";S$(I);:A=A+1:I=I+1:B=B+1: GOTO 1050
1135 BOX 1,R*15+R,159,(R+1)*15+R,1,2:RETURN
1140 ON A+R GOTO 1400,1600,1700,1800,1900,2000
1199 GOTO 1056

1200 REM 开平方
1205 BS=1
1220 IF BKS=0 or BKS=1 THEN BS=BKS:A$=STR$(BKS):1270
1230 FOR I=2 TO INT(SQR(BKS)):IF BKS/(I*I)=INT(BKS/(I*I)) THEN BS=BS*I:BKS=BKS/(I*I)
1240 NEXT
1250 A$="":IF BS<>1 THEN A$=A$+STR$(BS)
1260 IF BKS<>1 THEN A$=A$+CHR$(7)+STR$(BKS)
1270 RETURN

1300 REM 公约
1310 GYS=1:GBS=1
1320 FOR I=2 TO YS(YS):KG=0: FOR K=1 TO YS
1330 IF YS(K)/I=INT(YS(K)/I) THEN KG=KG+1:YS(K)=YS(K)/I
1340 NEXT K:IF KG=YS THEN GYS=GYS*I
1350 IF KG>0 THEN GBS=GBS*I:1320
1370 NEXT I
1375 ON A+R GOTO 0,0,1750,0,1980
1380 REM 公倍
1390 FOR I=1 TO YS:GBS=GBS*YS(I):NEXT :RETURN
1400 T$="aX"+CHR$(11)+"+bX+c=0":TEXT :PRINT T$:LOCATE 5,1:PRINT WR$;
1410 LOCATE 2,1:INPUT "a=";BA
1420 LOCATE 2,1:INPUT "b=";BB
1430 LOCATE 2,1:INPUT "c=";BC
1440 BD=BB*BB-4*BA*BC:IF BD<0 or BA=0 THEN LOCATE 2,1:PRINT "方程无实根":BEEP: GOTO 1490
1450 X1=(-BB+SQR(BD))/(2*BA):X2=(-BB-SQR(BD))/(2*BA)
1460 LOCATE 2,1:IF X1<>X2 THEN PRINT "方程有两不同实数根" ELSE PRINT "方程有两相同实数根"
1470 IF X1<>X2 THEN PRINT "X";CHR$(129);"=";X1,"X";CHR$(130);"=";X2: GOTO 1490
1480 PRINT "X";CHR$(129);"=X";CHR$(130),"  =";X1
1490 S=1: GOTO 1500
1500 LOCATE 5,1:PRINT "继续[Y/N]";:LOCATE 5,6
1510 M$=INKEY$:IF M$=CHR$(13) or M$="y" THEN 1530
1520 IF M$="n" THEN 1010 ELSE BEEP:1510
1530 ON A+R GOTO 1400,1600,1700,1800,1900,2000
1540 GOTO 1510
1550 GRAPH:PRINT "MathTools By LOJA","(C)LOJA Studio","http://www.loja.cn/":END
1600 S=2:T$="X="+CHR$(11)+CHR$(7)+"y":TEXT :PRINT T$:LOCATE 5,1:PRINT WR$;
1610 LOCATE 2,1:INPUT "y=";BY:IF BY=0 or BY<0 THEN BEEP:1610
1620 BKS=BY:PRINT "计算中…":GOSUB 1200
1630 LOCATE 3,1:PRINT "X=";A$;: GOTO 1500
1700 S=3:T$="X=a/b":TEXT :PRINT T$:LOCATE 5,1:PRINT WR$;
1710 LOCATE 2,1:INPUT "a=";BA:IF BA=0 THEN BEEP:1710
1720 LOCATE 2,1:INPUT "b=";BB:IF BB=0 THEN BEEP:1720
1730 LOCATE 2,1:PRINT " ";BA;"/";BB,"计算中…"
1740 BC=INT(BA/BB):BA=INT(BA-BC*BB+0.5):YS=2:YS(1)=BA:YS(2)=BB:IF BA>0 THEN  GOTO 1300
1750 LOCATE 3,1:PRINT "=";
1760 IF BC<>0 THEN PRINT BC;
1770 IF BA<>0 THEN PRINT "(";BA/GYS;"/";BB/GYS;")";
1790 GOTO 1500
1800 S=4:T$="请输入已知边长":TEXT :PRINT T$:LOCATE 5,1:PRINT WR$;
1810 LOCATE 2,1:PRINT "第一边:";:INPUT "";BA:IF BA<=0 THEN BEEP:1810
1820 LOCATE 2,1:PRINT "第二边:";:INPUT "";BB:IF BB<=0 THEN BEEP:1820
1825 LOCATE 2,1
1830 PRINT "直角边: 计算中…":BKS=ABS(BA*BA-BB*BB):GOSUB 1200
1835 LOCATE 2,8:IF A$="0" THEN PRINT "不存在!" ELSE PRINT A$
1840 PRINT "斜  边: 计算中…":BKS=BA*BA+BB*BB:GOSUB 1200:LOCATE 3,8
1845 LOCATE 3,8:IF A$="0" THEN PRINT "不存在!" ELSE PRINT A$
1850 GOTO 1500
1900 T$="请输入已知数":TEXT :PRINT T$:LOCATE 5,1:PRINT WR$;:YS=1
1910 LOCATE 2,1:PRINT "已知数";RIGHT$("0"+STR$(YS),2);"="," "," ";
1920 LOCATE 2,10:INPUT "";YS(YS):IF YS(YS)<=0 THEN BEEP:1920
1930 IF YS(YS)<YS(YS-1) THEN SWAP YS(YS),YS(YS-1)
1940 LOCATE 4,1:PRINT "继续输入[Y/N]":LOCATE 4,10
1950 M$=INKEY$:IF M$=CHR$(13) or M$="y" THEN YS=YS+1:1910
1960 IF M$="n" THEN 1970 ELSE BEEP:1950
1970 LOCATE 2,1:PRINT "计算中…"," "," ";: GOTO 1300
1980 GOSUB 1380:LOCATE 2,1:PRINT "公约数:";GYS,"公倍数:";GBS
1990 GOTO 1500
2000 T$="请输入待测数值":TEXT :PRINT T$:LOCATE 5,1:PRINT WR$;
2010 LOCATE 2,1:INPUT ">";BA:IF BA<>INT(BA) or BA<2 THEN BEEP:2010
2020 LOCATE 3,1:PRINT "计算中…"
2030 A$="":IF BA=2 THEN A$="2":2060
2040 FOR I=2 TO INT(SQR(BA)):IF BA/I=INT(BA/I) THEN A$=A$+STR$(I)+CHR$(140):BA=BA/I:2040
2050 NEXT :IF BA<>1 THEN A$=A$+STR$(BA) ELSE A$=LEFT$(A$,LEN(A$)-1)
2060 IF A$=STR$(BA) THEN LOCATE 2,1:PRINT BA;"为质数"," ";: GOTO 1500
2070 LOCATE 3,1:PRINT "=";A$;: GOTO 1500

#include<iostream>
#include<cmath>
using namespace std;
void equation(void); //解一元二次方程函数定义
void square_root(void); //求最简平方根函数定义
void fraction(void); //求最简分数函数定义
void pythagorean_proposition(void); //勾股定理函数定义
void GCD_LCM(void); //公倍、公约求解函数定义
void prime_number(void);  //质数判断函数定义
void point_distance(void);  //平面点距函数定义
void circle(void);  //有关圆和球的计算函数定义
void finish(void);  //运算结束语句函数定义
const double PI=3.1415926; //定义常量PI

//主菜单
int main(void)
{
    int x=0;
    system("cls"); //清除屏幕
    cout<<"==============================<理科工具集>=============================="<<endl<<endl
        <<"                            1.解一元二次方程"<<endl
        <<"                            2.求最简平方根"<<endl
        <<"                            3.求最简分数"<<endl
        <<"                            4.用勾股定理求直角边"<<endl
        <<"                            5.公倍、公约求解"<<endl
        <<"                            6.质数判断"<<endl
        <<"                            7.求平面内两点距离"<<endl
        <<"                            8.有关圆和球的计算"<<endl
        <<"                            0.退出本程序"<<endl
        <<"                                                    By LOJA"<<endl
        <<"                                           (C)LOJA Studio 2002-2006"<<endl
        <<"                                               http://www.loja.cn "<<endl
        <<"请输入各项对应的数字(0-9)：";
    cin>>x;
    cout<<endl;

    if (x==0)
    {
        cout<<"感谢使用！再见～～"<<endl;
        cin.get(); //按键等待，获取回车
        cin.get();
        return 0;
    }
    else if (x==1)
    {
        equation();
        main();
        return 0;
        
    }
    else if (x==2)
    {
        square_root();
        main();
        return 0;
        
    }
    else if (x==3)
    {
        fraction();
        main();
        return 0;
        
    }
    else if (x==4)
    {
        pythagorean_proposition();
        main();
        return 0;
        
    }
    else if (x==5)
    {
        GCD_LCM();
        main();
        return 0;
        
    }
    else if (x==6)
    {
        prime_number();
        main();
        return 0;
        
    }
    else if (x==7)
    {
        point_distance();
        main();
        return 0;
        
    }
    else if (x==8)
    {
        circle();
        main();
        return 0;
        
    }
    cout<<"输入选择项错误！\a"<<endl;
    cin.get();
    cin.get();
    main();
    return 0;
}

//解一元二次方程函数
void equation(void)
{
    double a=0,b=0,c=0,drt=-1,x1,x2;
    cout.setf(ios_base::fixed,ios_base::floatfield); //强制使用定点表示法
    while (drt<0)
    {
    system("cls");
    cout<<"------------------------------[解一元二次方程]------------------------------"<<endl
        <<"                            一般式  a*X*X+b*X+C=0 (a>0)"<<endl<<endl;
    while (a==0)
        {
         cout<<"                   请输入a=";
         cin>>a;
         if (a==0){     cout<<"                   a不能为零,请重新输入a !\a"<<endl<<endl;}
        }
    cout<<"                   请输入b=";
    cin>>b;
    cout<<"                   请输入c=";
    cin>>c;
    drt=b*b-4*a*c;
    if (drt<0) {     cout<<"                   方程无实根,请重新输入!\a"<<endl<<endl; a=0; cin.get(); cin.get();}
    }
    x1=(-b+sqrt(drt))/(2*a);
    x2=(-b-sqrt(drt))/(2*a);
    if (x1==x2)
    {
        cout<<"        方程有两相等实数根 X="<<x1<<endl;
        finish();
    }
    else
    {
        cout<<"        方程有两不相等实数根 X1="<<x1<<endl<<"                     X2="<<x2<<endl;
        finish();
    }
    
    return;
}



//求最简平方根函数

void square_root(void)
{
    int i;
    double x=-1,ans=1;
    cout.unsetf(ios_base::fixed);//取消强制使用定点表示法
    while (x<0)
    {
    system("cls");
    cout<<"------------------------------[求最简平方根]------------------------------"<<endl
    <<"                            y=√ˉx"<<endl<<endl
    <<"                      请输入x=";
    cin>>x;
    if (x<0) {     cout<<"                   被开方数不能小于为零,请重新输入x !\a"<<endl<<endl; x=-1; cin.get(); cin.get();}
    }
    if (x<=1)
    {
        cout<<"                      结果  y="<<x<<endl;
        cin.get();
        cin.get();
        return;
    }
    for (i=2;i<=sqrt(x);++i)
    {
        while (int(x)%(i*i)==0)
        {
            ans*=i;
            x/=(i*i);
        }

    }
    cout<<"                      结果  y=";
    if (ans!=1)
    {
        cout<<ans;
    }
    if (x!=1)
    {
        cout<<"√ˉ"<<x<<endl;
    }
    finish();
    return;
}
//求最简分数函数
void fraction(void)
{
    unsigned int i,a=0,b=0;
    cout.unsetf(ios_base::fixed);//取消强制使用定点表示法
    while (a<1||b<1)
    {
    system("cls");
    cout<<"--------------------------------[求最简分数]--------------------------------"<<endl
        <<"                            y=a/b"<<endl<<endl
        <<"                      请输入a=";
    cin>>a;
    cout<<endl<<"                      请输入b=";
    cin>>b;
    cout<<endl;
    if (a<1) {     cout<<"           为了使“求最简分数”运算有意义，分子不能小于1。请重新输入a ！\a"<<endl<<endl; a=0; cin.get(); cin.get();}
    if (b<1) {     cout<<"           为了使“求最简分数”运算有意义，分母不能小于1。请重新输入b ！\a"<<endl<<endl; b=0; cin.get(); cin.get();}
    }
    for (i=2;i<=a;++i)
    {
        while (a%i==0&&b%i==0)
        {
            a/=i;
            b/=i;
        }

    }
    if (a==1&&b==1)
    {
        cout<<"                      结果  y=1"<<endl;
    }
    else if (b==1)
    {
        cout<<"                      结果  y="<<a<<endl;
    }
    else
    {
        cout<<"                      结果  y="<<a<<"/"<<b<<endl;
    }
    finish();
    return;
}

//勾股定理函数
void pythagorean_proposition(void)
{
    unsigned int c=0,a=0,b=0;
    cout.unsetf(ios_base::fixed);//取消强制使用定点表示法
    while (a<1||b<1)
    {
    system("cls");
    cout<<"----------------------------[用勾股定理求直角边]----------------------------"<<endl
        <<"                            c*c=a*a+b*b"<<endl<<endl
        <<"      请输入直角三角形的一斜边长a=";
    cin>>a;
    cout<<endl<<"      请输入直角三角形的一斜边长b=";
    cin>>b;
    cout<<endl;
    if (a<1) {     cout<<"           三角形边长不能小于1。请重新输入a ！\a"<<endl<<endl; a=0; cin.get(); cin.get();}
    if (b<1) {     cout<<"           三角形边长不能小于1。请重新输入b ！\a"<<endl<<endl; b=0; cin.get(); cin.get();}
    }
    cout<<"                        结果  直角边长c="<<sqrt(float(a*a+b*b))<<endl;
    finish();
    return;
}
//公倍、公约求解函数
void GCD_LCM(void)
{
    unsigned int tmp,a=0,b=0,ba,bb;
    cout.unsetf(ios_base::fixed);//取消强制使用定点表示法
    while (a<1||b<1)
    {
    system("cls");
    cout<<"------------------------------[公倍、公约求解]------------------------------"<<endl<<endl
        <<"                请输入一个正整数a=";
    cin>>a;
    cout<<endl<<"                请输入一个正整数b=";
    cin>>b;
    cout<<endl;
    if (a<1) {     cout<<"             请重新输入一个正整数a ！\a"<<endl<<endl; a=0; cin.get(); cin.get();}
    if (b<1) {     cout<<"             请重新输入一个正整数b ！\a"<<endl<<endl; b=0; cin.get(); cin.get();}
    }
    ba=a;
    bb=b;
    if (a>b)
    {
        tmp=a;
        a=b;
        b=tmp;
    }


    //展转相除法
    while(b)
    {
        tmp=a%b;
        a=b;
        b=tmp;
    }
    cout<<"                a和b的最大公约数为："<<a<<endl<<endl
        <<"                a和b的最小公倍数数为："<<ba*bb/a<<endl;
    finish();
    return;
}

//质数判断函数
void prime_number(void)
{
    unsigned int a=0,i;
    bool b=1;
    cout.unsetf(ios_base::fixed);//取消强制使用定点表示法
    while (a<1)
    {
    system("cls");
    cout<<"--------------------------------[质数判断]--------------------------------"<<endl<<endl
        <<"                请输入一个正整数a=";
    cin>>a;
    cout<<endl;
    if (a<1) {     cout<<"             请重新输入一个正整数a ！\a"<<endl<<endl; a=0; cin.get(); cin.get();}
    }
    if (a==1) { cout<<"                1既不是质数，也不是合数。"<<endl; finish(); return;}
    for (i=2;i<a;++i)
    {
        if (a%i==0) {b=0; break;}
    }
    if (b)
    {
    cout<<"                a 是一个质数"<<endl;
    }
    else
    {
    cout<<"                a 不是一个质数"<<endl;
    }
    finish();
    return;
}

//平面点距函数
void point_distance(void)
{
    double x1,y1,x2,y2;
    cout.setf(ios_base::fixed,ios_base::floatfield); //强制使用定点表示法
    system("cls");
    cout<<"----------------------------[求平面内两点距离]----------------------------"<<endl<<endl
        <<"                请输入一个点的横坐标x1=";
    cin>>x1;
    cout<<"                请输入一个点的纵坐标y1=";
    cin>>y1;
    cout<<"                请输入另一个点的横坐标x2=";
    cin>>x2;
    cout<<"                请输入另一个点的纵坐标y2=";
    cin>>y2;
    cout<<"                这两点之间的距离为："<<sqrt(pow((x1-x2),2)+pow((y1-y2),2))<<endl;

    finish();
    return;
}



//有关圆和球的计算函数
void circle(void)
{
    long int r=-1;
    cout.setf(ios_base::fixed,ios_base::floatfield); //强制使用定点表示法
    while (r<0)
    {
    system("cls");
    cout<<"------------------------------[有关圆和球的计算]------------------------------"<<endl<<endl
        <<"                请输入圆或者球的半径 r=";
    cin>>r;
    cout<<endl;
    if (r<0) {     cout<<"             请重新输入一个非负整数 r ！\a"<<endl<<endl; r=0; cin.get(); cin.get();}
    }
    cout<<endl<<"                圆的周长："<<2*PI*r<<endl
        <<"                圆的面积："<<PI*r*r<<endl<<endl
        <<"                球的面积："<<4*PI*r*r<<endl
        <<"                球的体积："<<4/3.0*PI*r*r*r<<endl;


    finish();
    return;
}


//运算结束语句
void finish(void)
{
    cout<<endl<<endl<<endl<<endl<<"                                          运算完成！按<Enter>返回主菜单。"<<endl;
    cin.get();
    cin.get();
    return;
}