CF893E Counting Arrays

最新推荐文章于 2024-09-15 10:20:23 发布

minato_yukina

最新推荐文章于 2024-09-15 10:20:23 发布

阅读量167

点赞数

文章标签：图论算法 c++

本文链接：https://blog.csdn.net/qq_36018057/article/details/126651487

版权

这篇博客探讨了如何利用数学中的分解因式和二项式系数解决实际问题。通过分析数字的质因数分解，将问题转化为多重集组合，巧妙地处理负数情况，得到正数结果。博主展示了如何使用C++实现相关算法，并通过实例解释了计算过程。文章以一段寓意深长的诗句开头，吸引读者进入数学的世界。

摘要由CSDN通过智能技术生成

在这里插入图片描述
不难想到对 $x$ 分解后,得到 $x=p_1^{c1}*..p_i^{c_i}$
每一项 $c_i$ 的分配都是独立的,问题转化为对每一项 $c_i$ 分配位置,类似多重集的组合数。
那么负数怎么理解呢,因为最后得到的数字是正数,就是挑0,2,4,6,8…个位置取符号,该系数是二项式系数的偶数项之和 $2^{n-1}$
在这里插入图片描述

/*
Stairs upon the temple
I climb and I crawl in
People travel millions of miles just to see it
But they never conquer this way
*/
#include<bits/stdc++.h>
using namespace std;
const int maxn = 1e6+5;
const int INF = 1e9+7;
typedef long long ll;
typedef pair<int,int> pii;
#define all(a) (a).begin(), (a).end()
#define pb(a) push_back(a)
vector<int> G[maxn];
const int mod = INF;
int vis[maxn];ll inv[2*maxn];
ll pw(ll a,ll b){
	ll ans = 1;a%=mod;
	while(b){
		if(b&1) ans = ans*a%mod;
		a =a*a%mod;b>>=1;
	}
	return ans%mod;
}
ll fac[maxn*2];
void init(int n){
	vis[1] = 1;
	for(int i=2;i<=n;i++){
		if(!vis[i]){
			for(int j=i;j<=n;j+=i) vis[j] = i;
		} 
	}
	fac[1] = 1;fac[0] = 1;
	for(int i=2;i<=2*n;i++) fac[i] = fac[i-1]*i%mod;
	inv[2*n] = pw(fac[2*n],mod-2);inv[0] = 1;
	for(int i=2*n-1;i>=1;i--) inv[i] = inv[i+1] * (i+1) %mod;
} 
ll C(int n,int m){
	return fac[n] * inv[m] % mod * inv[n-m] %mod;
}
ll get(int &x){
	int cur = vis[x];ll cnt = 0;
	while(x!=0&&x%cur==0) cnt+=1,x/=cur;
	return cnt;
}
int main(){
    ios::sync_with_stdio(false);
	cin.tie(0);
	cout.tie(0);
	int MX = 1e6;
	init(MX);
	int T;cin>>T;
	while(T--){
		int x,n;cin>>x>>n;
		ll ans = 1;
		while(x!=0&&x!=1){
			ll cnt = get(x);
			if(cnt==0) break;
			ans = ans * C(n+cnt-1,n-1) %mod;	
	    }
	    ans = ans * pw(2,n-1)%mod;
	    cout<<ans<<"\n";
	}
}