牛客周赛 46 F 祥子拆团

原题链接:F-祥子拆团

题目大意:多测,每次给a,b,要将a分解成b个数相乘,问有多少种分的方法。

思路:对a进行质因数分解,对每一个质数计数,然后分到b个篮子里面,允许篮子里面没有数,对全部的方案数相乘就是答案。d代表每个质数的数量,使用隔板法,如果允许篮子里面没有数的话,那么公式就是C\binom{b-1}{d+b-1},如果不允许没有数可以理解成先给每个篮子放一个数,那么公式就是C\binom{b-1}{d-1}。因为允许没有数所以使用第一个公式。

//冷静,冷静,冷静
//调不出来就重构
#include<bits/stdc++.h>
#define endl '\n'
using namespace std;
typedef long long ll;
typedef long double ld;
typedef pair<ll,ll> pii;
const int N=1e6+10,mod=1000000007;
ll ksm(ll a,ll b)
{
	ll ans=1;
	do
	{
		if(b&1)ans*=a;
		a*=a;b>>=1;
		a%=mod;ans%=mod;
	}while(b);
	return ans;
}
int main()
{
	ios::sync_with_stdio(0),cin.tie(0),cout.tie(0);
	ll t;cin>>t;
	while(t--)
	{
		ll a,b;cin>>a>>b;
		ll ans=1;
		for(int i=2;i<=a/i;i++)
		{
			ll d=0;
			while(a%i==0)
			{
				d++;
				a/=i;
			}
			ll sum=1;
			for(int j=b,k=1;j<=d+b-1;j++,k++)
			{
				sum=sum*j;sum%=mod;
				sum=sum*ksm(k,mod-2)%mod;
			}
			ans*=sum;ans%=mod;
		}
		if(a>1)
		{
			ll d=1,sum=1;
			for(int j=b,k=1;j<=d+b-1;j++,k++)
			{
				sum=sum*j;sum%=mod;
				sum=sum*ksm(k,mod-2)%mod;
			}
			ans*=sum;ans%=mod;
		}
		cout<<ans<<endl;
	}
    return 0;
}

  • 13
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值