P1009阶乘和 P1217回文质数

P1009阶乘和

#include<stdio.h>

int a[1001],b[1001],c[1001];

int main(){
	a[0]=1;
	int n;
	scanf("%d",&n);
	for(int i=1;i<=n;i++){
		for(int j=0;j<=1000;j++){
			b[j]=a[j]*i;
		}
		for(int j=0;j<=1000;j++){
			if(b[j]>=10){
				b[j+1]+=b[j]/10;
				b[j]%=10;
			}
		}
		for(int j=0;j<=1000;j++){
			a[j]=b[j];
		}
		for(int j=0;j<=1000;j++){
			c[j]+=b[j];
			if(c[j]>=10){
				c[j+1]+=c[j]/10;
				c[j]%=10;
			}
		}
	}
	int here;
	for(int j=999;j>=0;j--){
		if(a[j]){
			here=j;
			break;
		}
	}
	for(;here>=0;here--){
		printf("%d",c[here]);
	}
	return 0;
}

高精度乘法的核心公式：

c[i+j]=a[i]*b[j]

高精度乘低精度时不必另开数组，直接乘并处理进位即可

P1217[USACO1.5]回文质数https://www.luogu.com.cn/problem/P1217

#include<bits/stdc++.h>
using namespace std;

bool vis[100000005];
int prime[5800000];

bool find(int k){
	char c[10];
	int h=k,count=0;
	while(h){
		h/=10;
		count++;
	}
	h=k;
	for(int i=0;i<count;i++){
		c[count-i-1]=h%10;
		h/=10;
	}
	int left=0,right=count-1;
	while(left<=right){
		if(c[left++]!=c[right--]){
			return false;
		}
	}
	return true;
} 

int main(){
	int left,right;
	cin>>left>>right;
	int count=0;
	for(int i=2;i<=right;i++){
		if(!vis[i]){
			prime[count++]=i;
		}
		for(int j=0;j<count&&i*prime[j]<=right;j++){
			vis[i*prime[j]]=true;
			if(i%prime[j]==0){
				break;
			}
		}
	}
	int here=2;
	for(int i=0;i<count;i++){
		if(prime[i]>=left){
			here=i;
			break;
		}
	}
	for(int j=here;j<count;j++){
		if(find(prime[j])){
			cout<<prime[j]<<'\n';
		}
	}
	return 0;
}

 但是这里完全不需要这么麻烦

我原本是嫌先构造回文数再判断的话情况太多

就试着用欧拉筛然后逐个判断了

虽然AC了，但是涉及到一些数学知识需要记录一下：

1. 偶数位回文数，除11外，必不为质数，因为其可以被11整除。

证明方法，我不会（理直气壮），自行查询

这样的话就可以把数据从三种状况转换为一种情况，这样的话就不必列举那么多种情况了。

（因为特判5，7，11）（从5开始是题目限制）

这个方法理论上不开O2优化的话是会TLE最后一个点的

当然最后一个点我可以提前算出来然后特判，这样也能过

也就是说本题并非一个完美的解法，最好的方法还是构造回文数然后判断。

2. 回文数构造方法，这个东西题目下面给提示了

//产生长度为5的回文数
for (d1 = 1; d1 <= 9; d1+=2) {    // 只有奇数才会是素数
     for (d2 = 0; d2 <= 9; d2++) {
         for (d3 = 0; d3 <= 9; d3++) {
           palindrome = 10000*d1 + 1000*d2 +100*d3 + 10*d2 + d1;//(处理回文数...)
         }
     }
 }