「一本通 6.2 练习 5」樱花

最新推荐文章于 2024-09-12 19:02:20 发布

Dh_virgil

最新推荐文章于 2024-09-12 19:02:20 发布

阅读量55

点赞数

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

本文链接：https://blog.csdn.net/Dh_virgil/article/details/133830742

版权

文章介绍了两段C++代码，用于解决一个关于给定数对(x,y)满足特定条件的问题，即它们的因子个数乘积对模M的结果。代码涉及求质数、因子分解和模运算。

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

传送门：「一本通 6.2 练习 5」樱花 - 题目 - Liuser's OJ

样例

样例输入

复制2

样例输出

复制3

样例说明

共有三个数对（x, y）满足条件，分别是 (3, 6)(4, 4) 和(6, 3) 。

代码1：

#include <iostream>
#include <cstring>
#include <vector>
using namespace std;
typedef long long LL;
typedef pair<int,int> PII;
const int N = 1000010,M = 1e9+7;
int prime[N],cnt;
bool st[N];
vector<int> v;
void GetPrime(int n)
{
    for(int i=2;i<=n;i++)
    {
        if(!st[i]) prime[++cnt] = i;
        for(int j=1;prime[j]<=n/i;j++)
        {
            st[i*prime[j]] = 1;
            if(i%prime[j]==0) break;
        }
    }
}
void get(int n)
{
    for(int i=1;i<=cnt;i++)
    {
        int t = n,sum = 0;
        while(t) sum+=t/prime[i],t/=prime[i];
        if(!sum) break;
        v.push_back(sum*2);
    }
}
int main() {
    int n;
    scanf("%d",&n);
    GetPrime(n);
    get(n);
    int res = 1;
    for(auto it:v)
        res = (LL)res*(it+1)%M;
    printf("%d\n",res);
    return 0;
}

代码2：

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#define int long long

using namespace std;
int vis[1000001];
int prime[1000001];
int cnt;
int n;
void ol(){
  for(int i=2;i<=n;++i){
    if(!vis[i]){
      prime[++cnt]=i;
    }
    for(int j=1;j<=cnt&&i*prime[j]<=n;++j){
      vis[i*prime[j]]=1;
      if(i%prime[j]) break;
    }
  }
}
int cntt[10000001];
void divide(int x){
  for(int i=1;prime[i]*prime[i]<=x;++i){
    while(x%prime[i]==0){
      cntt[prime[i]]++;
      x/=prime[i];
    }
  }
  if(x!=1){
  cntt[x]++;
  }
  return ;
}
int ans=1;
int mod=1000000007;
signed main(){
  scanf("%d",&n);
  ol();
  for(int i=1;i<=n;++i){
    divide(i);
  }
  for(int i=1;i<=cnt;++i){
    ans=ans*((cntt[prime[i]]*2+1));
    ans%=mod;
  }
  cout<<(ans);
  return 0;
}