2019.12月第四周练习总结

最新推荐文章于 2024-04-30 23:53:05 发布

weixin_44115550

最新推荐文章于 2024-04-30 23:53:05 发布

阅读量93

点赞数

文章标签：算法

本文链接：https://blog.csdn.net/weixin_44115550/article/details/103835321

版权

12月第四周练习总结

第一题洛谷P4931
题意:
　　　有n排座位,每排有两个座位,现在有 n对情侣入座,恰好有K对情侣坐一起的情况数
题解:
设 f(x) 为 x 对情侣错开的方案数;
所以 ans(n,k) = C(n,k) * A(n,k) * 2^k * f(n-k);
组合数和阶乘数和K次幂都好求,关键点为求f(n-k);
f(x) 为错位排序,我们可以考虑第一排的情况
第一排有3种情况,两男,两女,一男一女;
假设第一排坐两个男生,那么的人可以分为2两种情况,即
　　第一种,两男生的女朋友坐在一起, 方案数为 f(x-2)*2(x-1)
　　第二种,不坐在一起,那我们可以把他们当做一对情侣,所以现在的　　
　　　　　　　方案数为　f ( x - 1 );
可以想到剩下的两个女生和一男一女和两个男生的方法是一样的
所以 f(x) = 4 x (x-1) * [ f ( x - 1) + 2 ( x - 1) * f ( x -2 ) ] ;
可以将数据预处理，所以最后的时间复杂度为O(nT)

using namespace std;
typedef long long LL;
const int N = 5e6+10;
const int lim = 5e6;
const int p = 998244353;
LL fac[N],inv[N],f[N],ksm[N];
void init() //预处理 阶乘,组合数,2的幂次和f(x);
{
	fac[0] = inv[0] = inv[1] = 1;
	f[0] = ksm[0] = 1;	f[1] = 0;
	for(int i=1;i<=lim;i++)
		fac[i] = 1ll*fac[i-1]*i%p;
	for(int i=2;i<=lim;i++)
		inv[i] = (p-p/i)*inv[p%i]%p;
	for(int i=1;i<=lim;i++)
		inv[i] = inv[i-1]*inv[i]%p; 
	for(LL i=2;i<=lim;i++)
		f[i] = (4*i%p*(i-1)%p*((2*(i-1)*f[i-2]+f[i-1])%p)%p)%p;
	for(int i=1;i<=lim;i++)
		ksm[i] = 1ll*2*ksm[i-1]%p;
}


LL A(int n,int m)
{
	return fac[n]*inv[n-m]%p;
}

LL C(int n,int m)
{
	return A(n,m)*inv[m]%p;
}

int main()
{
	init();
	int T;
	scanf("%d",&T);
	while(T--)
	{
		int n,k;
		scanf("%d%d",&n,&k);
		LL ans = C(n,k)*A(n,k)%p*ksm[k]%p*f[n-k]%p;
		printf("%lld\n",ans);
	}
	return 0;
}

第二题洛谷P1891

题意：
　　LCM(1,N)+LCM(2,N)+…+LCM(N,N);
题解:
　在这里插入图片描述

#include<bits/stdc++.h>
using namespace std;
#define re register
typedef long long LL;
const int N = 1e6+10;
const int lim = 1e6;
LL phi[N],sum[N];
int prime[N];
bool vis[N];

void init()　//打表
{
	phi[1] = 1;
	int cnt = 0;
	for(re int i=2;i<=lim;i++)
	{
		if( ! vis[i] )
		{
			prime[++cnt] = i;
			phi[i] = i-1;
		}
		for(re int j=1;j<=cnt&&prime[j]*i<=lim;j++)
		{
			vis[prime[j]*i] = 1;
			if( i%prime[j] == 0 )
			{
				phi[i*prime[j]] = phi[i] * prime[j];
				break;
			}
			phi[i*prime[j]] = phi[prime[j]] * phi[i];
		}
	}
	for(re int i=2;i<=lim;i++)
		phi[i] = 1ll*(phi[i]*i)/2;
	for(re int i=1;i<=lim;i++)
		sum[i] = 1;
	for(re int i=2;i<=lim;i++) // 类似于素数的O(nlogn)筛
		for(re int j=i;j<=lim;j+=i)
			sum[j] += phi[i];
}

int main()
{
	init();
	int T;
	scanf("%d",&T);
	while(T--)
	{
		int n;
		scanf("%d",&n);
		printf("%lld\n",1ll*n*sum[n]);
	}
	return 0;
}

题目三洛谷P2557

题意:
　　输入ｎ和 k 两个数,求出 n 的 k 次的所有因子和;
题解：
　　根据唯一分解定理将 n 分解成 m 个质因数的乘积,即
　　n=p1^a1 * p2^a2 . . .*pn^an;
　　所以因子数的和为
　　q1^0+q11+q1^2…q1a1）＊（q2^0+q21+q2^2…q2a2）＊
　　…＊（qn^0+qnn+qn^2…qnan）
　　由于本题 n <= 2e16,k<=20;
　　所以这里需要用到高精度的乘法,加法和除法

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N = 230;
const int M = 65536;
LL n,k;

string Gmul(string aa,string bb)
{
	int a[N],b[N],c[N]={};
	int lena = aa.length();
	int lenb = bb.length();
	for(int i=0;i<lena;i++)
		a[lena-i] = aa[i] - '0';
	for(int i=0;i<lenb;i++)
		b[lenb-i-1] = bb[i] - '0';
	for(int i=1;i<=lena;i++)
	{
		for(int j=0;j<lenb;j++)
			c[i+j] += a[i]*b[j];
	}
	for(int i=1;i<lena+lenb;i++)
	{
		c[i+1] += c[i]/10;
		c[i] %= 10;
	}
	int k=lena+lenb;
	while( c[k] >= 10 )
	{
		c[k+1] += c[k]/10;
		c[k] %= 10;
		k++;
	}
	bool f = 0;
	string ans;
	for(int i=k;i>=1;i--)
	{
		if( c[i] != 0 )
			f = 1;
		if( f == 1 )
			ans += c[i]+'0';
	}
	return ans;
}

string Gdiv(string a,int b)
{
	string ans;
	int d=0,f=0;
	int len = a.length();
	for(int i=0;i<len;i++)
		a[i] -= '0';
	for(int i=0;i<len;i++)
	{
		d = d*10+a[i];
		if( d >= b )
			f = 1;
		if( f == 1 )
		{
			ans += d/b + '0';
			d %= b;
		}
	}
	return ans;
}

string Gksm(string a,int b)
{
//	cout<<a<<" "<<b<<"..."<<endl;
	string ans = "1";
	while( b )
	{
		if( b&1 )
			ans = Gmul(a,ans);
		b >>= 1;
		a = Gmul(a,a);
	}
	return ans;
}

string Gsub( string aa )
{
	int len = aa.length();
	aa[len-1] -= 1;
	return aa;
}

int main()
{
	scanf("%lld%lld",&n,&k);
	string ans = "1";
	for(LL i=2; i<=M && n>1;i++)
	{
		if( n % i == 0 )
		{
			LL b = 0, t = i;
			while( n%i == 0 )
			{
				b++;
				n /= i;
			}
//			printf("%lld___%lld\n",i,b);
			string aa;
			stack<int> ss;
			while( t )
			{
				ss.push(t%10);
				t /= 10;
			} 
			while( !ss.empty() )
			{
				aa += ss.top()+'0';
				ss.pop();
			}
			b *= k;
//			printf("%lld\n",b);
			aa = Gksm(aa,b+1);
//			cout<<aa<<".._"<<endl;
			aa = Gsub(aa);
//			cout<<aa<<"__"<<endl;
			aa = Gdiv(aa,i-1);
//			cout<<aa<<"_//"<<endl;
			ans = Gmul(ans,aa);
		}
	}
	cout<<ans<<endl;
	return 0;
}

题目四洛谷P2675

题意:
　　输入 n ,将1-n个数按任意顺序放入a[1][ i ], a[i][j] = a[i-1][j-1]+a[i-1][j];
　使得 a[n][1]最大

题解:
　　当我们手模一下会发现每个数对最终答案的贡献系数为杨辉三角系
　　数
　　例如:
　　a　　　　　　ｂ　　　　　　　ｃ　　　　　　　ｄ
　　　　ａ＋ｂ　　　　ｂ＋ｃ　　　　　ｃ＋ｄ
　　　　　　　ａ＋２ｂ＋ｃ　　ｂ＋２ｃ＋ｄ
　　　　　　　　　　ａ + 3b + 3c + d
　　所以本题只要将大的数与杨辉三角大系数相乘就可得出答案

#include<bits/stdc++.h>
using namespace std;
const int p = 10007;
const int N = 10010;
int fac[N],inv[N];
void init()
{
	inv[1] = fac[1] = fac[0] = inv[0] = 1;
	for(int i=2;i<p;i++)
		fac[i] = 1ll*fac[i-1]*i%p;
	for(int i=2;i<p;i++)
		inv[i] = (p-p/i)*inv[p%i]%p;
	for(int i=2;i<p;i++)
		inv[i] = inv[i-1]*inv[i]%p;
}
int C(int n,int m)
{
	if( m > n )
		return 0;
	if( m < p && n < p )
		return fac[n]*inv[m]%p*inv[n-m]%p;
	return  C(n/p,m/p)%p*C(n%p,m%p)%p;
}

int main()
{
	init();
	int n,ans=0;
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
	{
		if( i&1 )
			ans = (ans+1ll*i*C(n-1,(i+1)/2-1)%p+p)%p;
		else
			ans = (ans+1ll*i*C(n-1,i/2-1)%p+p)%p;
//		printf("%d__%d\n",i,ans);
	}
	printf("%d\n",ans%p);
	return 0;
}

第五题洛谷P1390

题意:
　　算出1~N这N个数中每任意两个不同的数的最大公约数的和

题解：
　　在这里插入图片描述
我们这里经过反演将函数简化，在将欧拉函数打表就可以求出公式的答案，由于每组数出现一次且不出现 i == j 的情况,所求我们先减去(1+n)*n/2在整体除2就是答案

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

const int N = 2e6+10;
const int lim = 2e6;
LL phi[N];
int prime[N];
bool vis[N];
int w[N],v[N];
void init()
{
	int cnt = 0;
	phi[1] = 1;
	for(int i=2;i<=lim;i++)
	{
		if( vis[i] == 0 )
		{
			phi[i] = i-1;
			prime[++cnt] = i;
		}
		for(int j=1;j<=cnt&&prime[j]*i<=lim;j++)
		{
			vis[i*prime[j]] = 1;
			if( i%prime[j] == 0 )
			{
				phi[i*prime[j]] = phi[i]*prime[j];
				break;
			}
			phi[i*prime[j]] = phi[prime[j]]*phi[i];
		}
	}
	for(int i=1;i<=lim;i++)
		phi[i] += phi[i-1];
}
int main()
{
	init();
	int n;
	scanf("%d",&n);
	LL ans=0;
	for(LL l=1,r;l<=n;l=r+1)
	{
		r = n/(n/l);
		ans += 1ll*(n/l)*(n/l)*(phi[r]-phi[l-1]);
	}
	printf("%lld\n",(ans-1ll*(n+1)*n/2)/2);
	return 0;
}

weixin_44115550

关注

0
点赞
踩
0

收藏

觉得还不错? 一键收藏
0
评论
2019.12月第四周练习总结

12月第四周练习总结第一题洛谷P4931题意:　　　有n排座位,每排有两个座位,现在有 n对情侣入座,恰好有K对情侣坐一起的情况数题解:设 f(x) 为 x 对情侣错开的方案数;所以 ans(n,k) = C(n,k) * A(n,k) * 2^k * f(n-k);组合数和阶乘数和K次幂都好求,关键点为求f(n-k);f(x) 为错位排序,我们可以考虑第一...
复制链接

扫一扫