【数学-高精度压位】中山纪念中学暑期游Day15——【NOIP2013模拟联考3】恭介的法则(rule)

顾玥_浅笑

于 2019-08-15 20:02:16 发布

阅读量293

点赞数

分类专栏：数论高精度文章标签：数学高精度压位

本文链接：https://blog.csdn.net/qq_36294918/article/details/99641351

版权

数论同时被 2 个专栏收录

6 篇文章 0 订阅

订阅专栏

高精度

1 篇文章 0 订阅

订阅专栏

前言

做完这题体会到：数学推导能力超重要啊！

题目

Time Limits: 1000 ms Memory Limits: 131072 KB

终于，在众亲们的奋斗下，最终boss 恭介被关进了库特设计的密室。正当她们松了一口气时，这个世界却发生了天翻覆地的变化：地面开始下沉，天空开始变成血红色，海水沸腾……一幅世界末日的图景。美鱼从她手中的古籍《若山牧水诗歌集》中发现了原因：白鸟はかなしからずや空の青海のあをにも染まずただよふ。大（xia）意（shuo）就是狡猾的恭介在创造这个世界的时候就篡改了法则。而这个法则的起源，就是一只生死之间的猫。这个猫被关在一个黑盒子里，盒子里有两个毒气罐，如果有任意一个毒气罐被打开那么猫将会被杀死，法则也能得到纠正。然而外界能控制的仅仅是这两个毒气罐被打开的概率。假设第一个毒气罐被打开的概率为1/x，第二个毒气罐为1/y（x,y 为正整数），那么当两个概率和为1/(n!)时，猫将会被莫名其妙地杀死。现在美鱼想知道，有多少对(x,y)可以让猫被莫名其妙杀死。

Input

一行，一个正整数n

Output

一行，满足题意的(x,y)对数。

Sample Input
6
Sample Output
135

【数据范围】

对于30%的数据 n<=6

对于60%的数据 n<=50

对于100%的数据 n<=700000

分析

由题目可得出等式： $\frac{1}{x}+\frac{1}{y}=\frac{1}{n!}$ ，要你求出（x,y）的方案数

变形方式有很多种，说说我自己推导的过程：

左右同乘 x 得： $1+\frac{x}{y}=\frac{x}{n!}$

移项得： $1=\frac{x}{n!}-\frac{x}{y}$

通分得： $\frac{x*(y-n!)}{y*n!}=1$

去分母得： x*(y-n!)=y*n!

左边只留x： $x=\frac{y*n!}{y-n!}$

因为x，y，n为正整数，所以 y*n! 为 (y-n!) 的倍数，且y>n!

令y=n!+k，则原始变为： $x=\frac{(n!+k)*n!}{(n!+k)-n!}$ ，即 $x=\frac{n!^{2}+k*n!}{k}$

因为k*n!一定是k的倍数，所以只需保证n!^2为k的倍数

问题转化为：求n!^2的因数个数

一个数可表示为 $p_{1}^{k1}*p_{2}^{k2}*...*p_{x}^{kx}$ ，

该数平方即为 $p_{1}^{2k1}*p_{2}^{2k2}*...*p_{x}^{2kx}$ ，

对于每个质因数pi，有选0个、选1个、...选2ki个共（2ki+1）种方案

总方案即为 $(2k_{1}+1)*(2k_{2}+1)*(2k_{3}+1)*...*(2k_{x}+1)$

就用线性筛，筛出质数，算出n!的各质因数的质数，把它们*2+1，最后全部乘起来即可

做乘法时用【高精度+压位】，不然直接【高精度】会TLE

（考试时直接打的高精度，不会压位QAQ）

丢些【高精度压位】的博客：

高精度压位

算法笔记-压位高精

考试无压位TLE代码

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
const int MAXN=7e4;
bool vis[MAXN+5];
ll prime[MAXN+5],cnt[MAXN+5];
ll a[MAXN+5],b[MAXN+5],c[MAXN+5],d[MAXN+5];
ll n,tot,lena,lenb,lend;
void read(ll &x)
{
    int f=1;x=0;char s=getchar();
    while(s<'0'||s>'9'){if(s=='-')f=-1;s=getchar();}
    while(s>='0'&&s<='9'){x=x*10+s-'0';s=getchar();}
    x*=f;
}
void Sieve()//线性欧拉筛 
{
	for(register int i=2;i<=n;i++)
	{
		if(!vis[i])
			prime[++tot]=i;
		for(register int j=1;j<=tot;j++)
		{
			if(i*prime[j]>n)
				break;
			vis[i*prime[j]]=1;
			if(i%prime[j]==0)
				break;
		}
	}
}
void Cal(ll x)
{
	for(register int i=1;i<=tot;i++)
	{
		if(prime[i]>x)
			break;
		while(x%prime[i]==0)
		{
			cnt[i]++;
			x/=prime[i];
		}
	}
}
void Change(ll tmp)
{
	memset(a,0,sizeof(a));
	lena=0;
	while(tmp>0)
	{
		lena++;
		a[lena]=tmp%10;
		tmp/=10;
	}
}
void Plus(ll tmp)//高精度乘法 
{
	Change(tmp);//A
	lenb=lend;
	memcpy(b,d,sizeof(d));
	memset(c,0,sizeof(c));//c暂存答案 
	ll res,t,result;
	for(register int i=1;i<=lenb;i++)
	{
		res=0;
		for(register int j=1;j<=lena;j++)
		{
			t=j+i;
			result=a[j]*b[i]+res+c[t-1];
			res=result/10;
			c[t-1]=result%10;
		}
		while(res)
		{
			c[t++]=res%10;
			res=res/10;
		}
	}
	while(c[t]==0)
		t--;
	memcpy(d,c,sizeof(c));
	lend=t;
}
int main()
{
	read(n);
	Sieve();
	for(register int i=2;i<=n;i++)
		Cal(i);
	lend=1;
	d[1]=1;
	for(register int i=1;i<=tot;i++)
		Plus(1LL*cnt[i]*2+1);
	for(int i=lend;i>=1;i--)
		printf("%lld",d[i]);
	return 0;
}

AC代码

高精度压八位

#include<cstdio>
#include<cmath>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
const int MAXN=7e5,mod=1e8;
ll a[MAXN+5];
bool vis[MAXN+5];
int n,len;
void work(int x)
{
	int k=0;
	for(int i=1;i<=len;i++)
	{
		a[i]=a[i]*x+k;
		k=a[i]/mod;
		a[i]%=mod;
	}
	if(k)
		a[++len]=k;
	while(a[len]>=mod)
	{
		a[len+1]=a[len]/mod;
		a[len]=a[len]%mod;
		len++;
	}
}
int main()
{
	scanf("%d",&n);
	len=1;
	a[1]=1;
	ll x=1;
	for(int i=2;i<=n;i++)
	{
		if(!vis[i])//筛质数 
		{
			ll y=n,z=0;
			while(y>=i)//算出质数指数 
			{
				y/=i;
				z+=y;
			}
			z=2*z+1;
			if(x*z>mod)//如果需要压位就压位 
			{
				work(x);
				x=1;
			}
			x*=z;
			for(int j=1;j<=n/i;j++)
				vis[i*j]=1;
		}
	}
	if(x>1)//还有没乘完的，继续乘 
		work(x);
	printf("%d",a[len]);
	for(int i=len-1;i>=1;i--)
		printf("%08lld",a[i]);
	return 0;
}