洛谷【P1619】 解一元二次方程的烦恼

我对模拟的理解：https://www.cnblogs.com/AKMer/p/9064018.html

题目传送门：https://www.luogu.org/problemnew/show/P1619

写完后我觉得我该告诉大家的第一句话就是：

不要深夜在洛谷写这种蓝色难度的模拟题（当然您如果够神写深蓝色难度的也是体会不到我这种蒟蒻今晚尝到的痛苦的。）

这题写的我真是有一种“日了出题人亲娘”爽快的感觉

换行就是这题最难的模拟，一下要换，一下不要换，特殊情况要特判。

然后判质数和分解质因数应该是入门选手都会的了，因为数字不大，所以可以根号\(n\)求法,不需要线性筛质数。

反正这题挺\(fake\)的，我一开始看完题目觉得今晚\(10\)点就可以睡觉了，然后就从\(21:50\)一直\(Wa\)到了\(22:41\)才\(A\)掉\(emmm\)

还有一句劝：不管写什么题都应该先想清楚再动手，不要觉得自己好像可以\(A\)就莽莽撞撞的开了，在码代码的过程中如果漏想了什么思路后面可能会很难想到。如果依靠数据发现了自己漏想的情况改起来也不如一开始就想到那种情况那么简单。而且万一在大型比赛遇到这种情况然后官方给的大数据比较水，那么你可能会以为自己\(A\)了然后开开心心的玩了几天然后就发现自己爆零了…………

总之，希望自己和大家都不要再有这种情况了吧。\(OI\)有风险，开码需谨慎。

时间复杂度：\(O(n+\sqrt{n})\)

空间复杂度：\(O(n)\)

代码如下：

#include <cmath>
#include <cstdio>
#include <cstring>
using namespace std;
#define ll long long
int top;
char s[25154];//25154是某人学号emmm
int pri[100],sum[100];
bool ispri(ll x) {
	if(x==1)return 0;//1特判
	for(int i=2;i<=trunc(sqrt(x));i++)
		if(x%i==0)return 0;//sqrt(x)复杂度判质数
	return 1;
}
void check(ll x) {
	top=0;ll num=x;
	if(x<2)return;
	for(int i=2;i<=trunc(sqrt(x));i++) {
		if(num%i==0&&ispri(i)) {
			pri[++top]=i;sum[top]=0;
			while(num%i==0)
				num/=i,sum[top]++;
		}
		if(num%(x/i)==0&&ispri(x/i)) {
			pri[++top]=x/i;sum[top]=0;
			while(num%(x/i)==0)
				num/=(x/i),sum[top]++;
		}//记得要判即使质数又是因数才能算
	}
}
int main() {
	while(1) {
		gets(s);//数字之间可以插入各种符号，包括空格，scanf直接挂。
		int len=strlen(s);ll num=0;//就算num>4e7也要判负数，所以必须用long long存下来
		printf("Enter the number=\n");
		for(int i=0;i<len;i++)
			if(s[i]>='0'&&s[i]<='9')num=(num<<1)+(num<<3)+s[i]-'0';
		if(!num)break;
		printf("Prime? ");
		if(ispri(num))puts("Yes!");else puts("No!");
		if(num>4e7) {
			printf("The number is too large!\n\n");//记得每一组数据之间都有一个空格，所以要换两次行
			continue;
		}check(num);//分解质因数
		for(int i=1;i<=top;i++) {
			if(i==1)printf("%lld=",num);
			printf("%d^%d",pri[i],sum[i]);
			if(i!=top)printf("*");
			if(i==top)puts("");//记得换行，这题恶心就恶心在换行
		}
		printf("\n");//还要换行，至于哪里要换行哪里不要自己可以好好想想
	}
	return 0;
}