[Luogu OJ P1619]解一元二次方程的烦恼

题目背景

JosephZheng在写数学作业的预习。他往往使用Casio来帮忙解一元二次方程。但是Casio有一个问题，就是当Δ=b²-4ac为一个大素数或大合数时，其开平方的结果会以小数显示，而不是老师要求的二次根式形式。JosephZheng很是苦恼，一遇到这种情况就要手动解方程。一天他再也忍不住了，于是打开了电脑，编了一个prime程序……于是悲剧的OI们就要跟着疯狂的JosephZheng一起编这个程序，呵呵……

 

题目描述

废话少说，给你一个大数N（不一定在int64范围内），让你进行素性判断，然后分解质因数。当然，初中数学题不可能有大于int64的数让你判断素性，因此超过范围的数可以忽略不计。为了让程序更加贴心，JosephZheng多了一些要求，会在输入输出中给出具体情况。

 

输入格式

一个大数N（N为非负整数），其中这个数的各个位数之间可以插入各种符号，例如1234可以为1 - 2alsdkjf3！@￥%！@@@##￥……！￥#-4等。你需要在这一长串乱码中找出这个要判断的数。输入数据可能有多组，如果读到一行没有数字的字串即结束。字符串长度可能大于255。

 

输出格式

在读入数据之前先输出“Enter the number=”，不要换行。
读入数据完毕后writeln。
然后输出“Prime? ”，问号后有一个空格，但不要换行。
如果是质数则输出“Yes!”否则输出“No!”。此时换行。
若果是质数就halt，若是小于2的数则在输出“No!”后也halt。若是合数则分解质因数。如果该数大于四千万则输出“The number is too large!”，然后halt。否则分解质因数。输出结果的方式在输出样例中会详细给出。每组数据之间空一行。

 

样例输入 

【输入样例1】
4
eed
【输入样例2】
2
end
【输入样例3】
-1
adfs
【输入样例4】
1234###24#@13#@￥！1
hehe
【输入样例5】
1.5
1
1234324123512343123
@~@~@~@

样例输出 

【输出样例1】
Enter the number=
Prime? No!
4=2^2

Enter the number=
【输出样例2】
Enter the number=
Prime? Yes!

Enter the number=
【输出样例3】
Enter the number=
Prime? No!

Enter the number=
【输出样例4】
Enter the number=
Prime? No!
The number is too large!

Enter the number=
【输出样例5】
Enter the number=
Prime? No!
15=3^1*5^1

Enter the number=
Prime? No!

Enter the number=
Prime? No!
The number is too large!

Enter the number=

此题本水题，只是很烦人、、、

以下来自本人题解：

个人认为题里的一堆halt除了判断输出末尾那个应该是halt，
别的（“若果是质数就halt，若是小于2的数则在输出“No!”后也halt。若是合数则分解质因数。如果该数大于四千万则输出“The number is too large!”，然后halt。”）
都不该那么写。
话说读入你本也可以写成while not eof do的、、

Var
  st:AnsiSTring;
  ch:char;
  Tt:Int64;
Function Prime(Tt:int64):Boolean;
  //注意细节，一开始习惯性写成Tt:longint，然后202半天不知道为什么
  Var
    i:longint;
  Begin
    if Tt<=1 Then Exit(False);
    For i:=2 to Trunc(Sqrt(Tt)) do
      if Tt mod i=0 then Exit(False);
    Exit(True);
  End;
Procedure Task();
  Var
    i,j:longint;
    xx:longint;
    Flag:boolean;
  Begin
    Write(Tt,'=');
    For i:=2 to Tt do
      Begin
        Flag:=False;
        if Tt<=1 Then Break;
        xx:=0;
        While Tt Mod i=0 Do
          Begin
            Inc(xx);
            Tt:=Tt Div i;
            Flag:=True;
          End;
        if Flag Then
          Begin
            Write(i,'^',xx);
            Break;//输出一个后退出for循环，进入带*输出阶段
          End;
      End;
   {一开始想怎么判断前面是否有*，后来发现分成两端是个不错的方法，方便快捷一次见效不反复}
    For j:=i to Tt do
      Begin
        Flag:=False;
        if Tt<=1 Then Break;
        xx:=0;
        While Tt Mod j=0 Do
          Begin
            Inc(xx);
            Tt:=Tt Div j;
            Flag:=True;
          End;
        if Flag Then
          Begin
            Write('*',j,'^',xx);
    ｛下面的标程是拿数组存的，我感觉直接输出比较好，省空间，尤其是解决了*号问题之后、、、｝
          End;
      End;
    Writeln;
  End;
Begin
  While True do //话说原来while true do 都写成一些好玩的东西、、谁有兴趣自己改着玩啊、、
    Begin
      sT:='';
      Write('Enter the number=');
      While Not Eoln do
        Begin
          Read(ch);
          if ch in ['0'..'9'] then
            st:=st+ch;
        End;
      Readln;
      Writeln;
      if st='' Then Halt;
      Write('Prime? ');
      Val(st,Tt);//字符变数字
      if Prime(Tt) Then Writeln('Yes!')
      Else
        Begin
          Writeln('No!');
          if Tt>40000000 Then Writeln('The number is too large!')
          //数清楚0，我一开始写成4000000了、、、、、
          else if Tt>=2 Then Task;
        End;
      Writeln;
    End;
End.

这个题交了4次才AC，第一次是Task()里面循环写了个如果Tt<=1 Then Exit 导致没有回车换行，0
第二次是又发现Prime(Tt)里面范围int64的问题，20
第三次是又发现40000000写成了4000000，80
第四次AC、、
总体说，确实是水题，但也确实很考察细心、、、、、
各位认真了···
（*）Done By Catch-22.S.Iris in Lensent CST Team.（*）