J - Prime Game (Gym - 101981J)

原题链接:J - Prime Game

题目大意:

给你n个数,问 ∑ i = 1 n ∑ j = i n f a c ( i , j )   \sum_{i=1}^{n} \sum_{j=i}^{n}fac(i, j)\ i=1nj=infac(i,j) 
fac(i, j)是从下标 i i i j j j,对应的数累乘的结果的质因数的个数。

解题思路:

单个数的质因数一定是一堆数相乘的质因数(这个用纸分开写一下就能明白了)。那么这道题就可以转变为去找这个数的质因数对答案的贡献度了。这里需要一个算法:分解质因数(不懂得可以去了解一下,不是很难这里就不讲了)。这里给一个例子来帮助理解:
3     11 2     31 2     5 4     7 53     7 3 83 3 \ \ \ 11 \\ 2 \ \ \ 31 \\ 2 \ \ \ 5 \\ 4 \ \ \ 7 \\ 53 \ \ \ 7 \\ 3 \\ 83 3   112   312   54   753   7383
这组数里面,第一个3的贡献度就是n - 1 + 1,到第二个三就是上一个三的位置到这个位置的距离乘上(n - i + 1),原因就是 2 , 31 2,31 231那一列用了一次(n - i + 1), 2 , 5 2,5 25也用一次,以此类推,就是 ( i − p o s ( a [ i ] ) ∗ ( n − i + 1 ) (i - pos(a[i]) * (n - i + 1) (ipos(a[i])(ni+1)

代码:

#include <bits/stdc++.h>

using namespace std;

#define IOS                      \
    ios::sync_with_stdio(false); \
    cin.tie(0);                  \
    cout.tie(0);
//#define For(i, k, n) for (int i = k; i < n; i ++)
//#define FOR(i, k, n) for (int i = k; i <= n; i ++)
#define P acos(-1);
typedef long long LL;
#define PII pair<int, int>
#define PDD pair<double, double>
#define Pll pair<LL, LL>
typedef unsigned long long ULL;
const double eps = 1e-7;
const int N = 1e6 + 10, M = 1e6 + 10, mod = 1000000007;
const int inf = 0x3f3f3f3f;
int dx[4] = {0, 1, -1, 0}, dy[4] = {1, 0, 0, -1};
int lowbit(int x) { return (x & -x); };

int a[N], pos[N];
vector<int> v[N];

void solve()
{
    int n;
    cin >> n;
    for (int i = 1; i <= n; i ++ ) cin >> a[i];

    for (int i = 1; i <= n; i ++ )
    {
        for (int j = 2; j <= a[i] / j; j ++ )
        {
            if (a[i] % j == 0) v[i].push_back(j);
            while (a[i] % j == 0) a[i] /= j;
        }
        if (a[i] != 1) v[i].push_back(a[i]);
    }

    // for (int i = 1; i <= n; i ++ ) cout << v[i].size() << '\n';
    LL ans = 0;
    for (int i = 1; i <= n; i ++ )
    {
        for (int j = 0; j < v[i].size(); j ++ )
        {
            ans += (LL)(i - pos[v[i][j]]) * (n - i + 1);
            pos[v[i][j]] = i;
        }
    }

    cout << ans << '\n';
}   

int main()
{
    IOS 
    int T = 1;
    // cin >> T;
    while (T -- )
    {
        solve();
    }
    return 0;
}
  
  ```cpp
  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值