【bzoj4146】 [AMPPZ2014]Divisors

最新推荐文章于 2018-08-16 17:03:00 发布

GEOTCBRL

最新推荐文章于 2018-08-16 17:03:00 发布

阅读量1.5k

点赞数

CC 4.0 BY-SA版权

分类专栏： AMPPZ 大暴力水题

本文链接：https://blog.youkuaiyun.com/GEOTCBRL/article/details/48935423

水题同时被 3 个专栏收录

25 篇文章

订阅专栏

大暴力

12 篇文章

订阅专栏

AMPPZ

9 篇文章

订阅专栏

本文深入探讨了一种在解决复杂算法问题时遇到困难的实例，从WA（Wrong Answer）到AC（Accepted）的转变过程。通过详细分析问题的根源，提出了解决方案并优化了算法实现，最终成功通过了评测。文章中还讨论了类似问题在CF（Codeforces）平台上的普遍性，以及如何有效利用空间优化技术，如筛法等，来提高程序效率。此外，文章强调了避免使用不切实际的方法，如错误的预处理步骤，并提供了实现过程中的一些技巧和经验教训。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

。。。。。。。
又WA又T一时爽

一开始以为 $O(n\log a)$ 过不了。。。。。。
后来发现。。。。。。是把小数据的时间也给用上了。。。。
bzoj的评测时间比较神奇。。。。。

直接像筛一样搞就行了。CF里面似乎很多这种题。
$H_{n}=\sum_{k=1}^{n}\frac{n}{k}=O(n\log n)$
直接统计

一开始T了之后YY了各种不靠谱的方法。。。比如先筛质数然后 $f[i]+=f[pri[j]*i]$ 等等。。。真是naiive…

咦那status里面他们是怎么做到空间这么小的……？

#include <bits/stdc++.h>
#define rep(i,a,b) for(int i=a;i<=b;i++)
#define per(i,a,b) for(int i=a;i>=b;i--)

inline int rd() {
    char c = getchar();
    while (!isdigit(c)) c = getchar() ; int x = c - '0';
    while (isdigit(c = getchar())) x = x * 10 + c - '0';
    return x;
}

const int maxm = 2000001;
const int maxn = 20000000;

typedef long long ll;

int n , cnt[maxn + 1] , mx , x;

void input() {
    n = rd();
    rep(i , 1 , n) cnt[x = rd()] ++ , mx = std::max(x , mx);
}

void solve() {
    ll ans = 0;
    rep(i , 1 , mx) {
        ans += (ll)(cnt[i] - 1) * (cnt[i]);
        for(int j = i + i;j <= mx;j += i) ans += (ll)cnt[j] * cnt[i];
    }
    printf("%lld\n" , ans);
}

int main() {
    #ifndef ONLINE_JUDGE
//      freopen("data.txt" , "r" , stdin);
    #endif
    input();
    solve();
    return 0;
}