CodeForces - 615D Multipliers 费马小定理＋组合数学

最新推荐文章于 2021-11-22 10:25:10 发布

自ran而然

最新推荐文章于 2021-11-22 10:25:10 发布

阅读量579

点赞数

分类专栏：数论文章标签：费马小定理快速幂

本文链接：https://blog.csdn.net/S1766434333/article/details/53039536

版权

数论专栏收录该内容

15 篇文章 0 订阅

订阅专栏

题目传送门：点击打开链接

题意：求出整数n的所有因子的乘积对1e9+7的模。

从题目中的给出的例子1·2·3·6 = 36 可以看出前面的式子可以写成（1*6)*(2*3)=6*6=36,而n＝6，所以我们可以想到答案可以是 (n)^m这种类型。m＝n的因子个数除以2.

那么n的因子个数怎么求呢？根据题目给的数据，我们可以知道n的素因子种类可相应的个数，那我们可以假设 n的一个素因子为p，他的个数为np，那么这个素因子和其他的素因子组合的时候，我们可以取p的个数为 0-np，也就是(np+1)种可能，其他的也是这样，所以他们一共有(np1+1)(np2+1)....(npk+1).. 种可能，其中npi代表n中第i个素因子的个数，这个结果也就是n的因子个数。

但是有一种是例外的那就是n可以开方的时候，例如输入 2 5 5，由此可以知道n为25，1*5*25=125 和上面的情况不一样，我们可以直接25^(3/2)*sqrt(25), 但是n太大了，sqrt(n)不怎么好求。只好换个思路。。。。

所以我们现在要以n的素因子为思路，从(np1+1)(np2+1)....(npk+1).. 中我们还可以看到对于第一个素因子p1他的指数可以是 0-np1，也就是确定p1的指数为t时，n的因子数还有(np2+1)....(npk+1)..所以我们可以得到n的因子连乘中 p1 的指数＝(0+1+2+3...np1)*(np2+1)....(npk+1)..=(np1*(np1+1)/2)*(np2+1)....(npk+1)..那么其他的素因子的指数也能得出。现在的问题是怎么求出对于每个pi素因子(np1+1)..(np(i-1)+1)*(np(i+1)+1)..(npk+1)..

第一种思路就是可以用两个数组l[],r[] ,l[i]表示的是(np1+1)*(np2+1)...(npi+1)%mod 的值，r[i]表示的是(npi+1)*(np(i+1)+1)....%mod 的值，那么对于每个素因子pi 就可以得到(npi*(npi+1)/2)*r[i+1]*l[i-1]

第二种思路就是改变一下乘的顺序，这个跟代码一起看比较好理解，对于素因子pi，先求出pi^[(npi*(npi+1)/2)*sum],sum表示(np1+1)*(np2+1)....(np(i-1)+1)%mod,然后在乘ans在对mod取余就可以了

对与pi^[(npi*(npi+1)/2)*sum]我们用费马小定理和快速幂就可以求了 pi^[(npi*(npi+1)/2)*sum]＝pi^[(npi*(npi+1)/2)*sum％(mod-1)]

代码为第二种：

<span style="font-family:Times New Roman, Times, serif;font-size:14px;">#include <iostream>
#include <cstring>
using namespace std;
typedef long long  ll;
const int maxn=2e5+100;
const int mod=1e9+7;
int m,cnt;ll sum;ll tol;
ll data[maxn],num[maxn];
ll quickmod(ll a,ll b)
{
    ll ans=1;
    while(b)
    {
        if(b&1)
            ans=(ans*a)%mod;
        b/=2;
        a=a*a%mod;
    }
    return ans;
}

void init()
{
    cnt=0;
    memset(num,0,sizeof(num));
    for(int i=0;i<m;++i)
    {
        cin>>data[cnt];
        ++num[data[cnt]];
        if(num[data[cnt]]==1)
            ++cnt;
    }
    //cout<<cnt<<endl;
}
void solve()
{
    tol=1;
    sum=1;
    for(int i=0;i<cnt;++i)
    {
        ll tmp=quickmod(data[i],num[data[i]]*(num[data[i]]+1)/2);
        sum=(quickmod(sum,num[data[i]]+1)%mod)*quickmod(tmp,tol)%mod;
        tol=tol*(num[data[i]]+1)%(mod-1);
    }
    cout<<sum<<endl;
}
int main(int argc, const char * argv[])
{
    while(cin>>m)
    {
        init();
        solve();
    }
    return 0;
}
</span>

Multipliers

Time Limit: 2000MS

Memory Limit: 262144KB

64bit IO Format: %I64d & %I64u

Submit Status

Description

Ayrat has number n, represented as it's prime factorization p_i of size m, i.e. n = p₁·p₂·...·p_m. Ayrat got secret information that that the product of all divisors of n taken modulo 10⁹ + 7 is the password to the secret data base. Now he wants to calculate this value.