Advertisement
El_GEMMY

dream team 3

Aug 20th, 2022
127
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C++ 4.49 KB | None | 0 0
  1. // Those who cannot remember the past are
  2. // condemned to repeat it (use DP -_-)
  3. // - George Santayana
  4.  
  5. #include <bits/stdc++.h>
  6. #include <ext/pb_ds/assoc_container.hpp>
  7. #include <ext/pb_ds/tree_policy.hpp>
  8.  
  9. using namespace std;
  10. using namespace __gnu_pbds;
  11.  
  12. #define all(v) v.begin(), v.end()
  13. #define rall(v) v.rbegin(), v.rend()
  14. #define ll long long
  15. #define ull unsigned long long
  16. #define MOD 1000000007
  17. #define PI acos(-1)
  18. #define ceil(a, b) (((a) / (b)) + ((a) % (b) ? 1 : 0))
  19. #define imin INT_MIN
  20. #define imax INT_MAX
  21. #define llmax LLONG_MAX
  22. #define llmin LLONG_MIN
  23. #define inf 2000000000
  24. #define nl '\n'
  25. #define ppcnt __builtin_popcount
  26. #define ppcntll __builtin_popcountll
  27. #define clz __builtin_clz
  28. #define clzll __builtin_clzll
  29. #define ctz __builtin_ctz
  30. #define ctzll __builtin_ctzll
  31. #define modulo(a, b, mod) ((((a) % (mod)) + ((b) % (mod))) % (mod))
  32. #define cnte(v, x) count(all(v), (x))
  33. #define mine(v) min_element(all(v))
  34. #define maxe(v) max_element(all(v))
  35. #define updmin(a, b) a = min(a, b)
  36. #define updmax(a, b) a = max(a, b)
  37. #define findmod(x, m) x = ((x) % (m) + (m)) % m
  38. #define getmod(x, m) ((x) % (m) + (m)) % (m)
  39. #define debug(x) cout << "x: " << (x) << nl;
  40. #define debug2(x, y) cout << "x: " << (x) << " y: " << y << nl;
  41. #define ordered_set tree<pair<int, int>, null_type, less<>, rb_tree_tag, tree_order_statistics_node_update>
  42. #define ordered_map tree<int, int, less<>, rb_tree_tag, tree_order_statistics_node_update>
  43.  
  44. //vector<int> dx = {0, 0, 1, -1, 1, 1, -1, -1}, dy = {1, -1, 0, 0, 1, -1, 1, -1};
  45. vector<int> dx = {0, 0, 1, -1}, dy = {1, -1, 0, 0};
  46.  
  47. template<typename T = int>
  48. istream &operator>>(istream &in, vector<pair<int, int>> &v) {
  49.     for (auto &[x, y]: v) in >> x >> y;
  50.     return in;
  51. }
  52.  
  53. template<typename T = int>
  54. istream &operator>>(istream &in, vector<T> &v) {
  55.     for (T &i: v) in >> i;
  56.     return in;
  57. }
  58.  
  59. template<typename T = int>
  60. ostream &operator<<(ostream &out, const vector<T> &v) {
  61.     for (const T &x: v)
  62.         out << x << ' ';
  63.     return out;
  64. }
  65.  
  66. template<typename T = pair<int, int>>
  67. ostream &operator<<(ostream &out, const vector<pair<int, int>> &v) {
  68.     for (auto &[x, y]: v) {
  69.         out << x << ' ' << y << nl;
  70.     }
  71.     return out;
  72. }
  73.  
  74. void Start_Crushing() {
  75.     ios_base::sync_with_stdio(false);
  76.     cin.tie(nullptr);
  77.     cout.tie(nullptr);
  78. #ifndef ONLINE_JUDGE
  79.     freopen("input.txt", "r", stdin);
  80.     freopen("output.txt", "w", stdout);
  81. #endif
  82. }
  83.  
  84. struct dsu {
  85.     vector<int> p, rank;
  86.  
  87.     explicit dsu(int size) {
  88.         p.resize(size + 1);
  89.         rank.resize(size + 1);
  90.  
  91.         iota(all(p), 0);
  92.     }
  93.  
  94.     int get(int x) {
  95.         return p[x] = x == p[x] ? x : get(p[x]);
  96.     }
  97.  
  98.     void join(int u, int v) {
  99.         u = get(u), v = get(v);
  100.  
  101.         if (u == v)
  102.             return;
  103.  
  104.         if (rank[u] == rank[v])
  105.             rank[u]++;
  106.  
  107.         if (rank[u] > rank[v])
  108.             p[v] = u;
  109.         else
  110.             p[u] = v;
  111.     }
  112. };
  113.  
  114. vector<vector<int>> factors;
  115.  
  116. void build(int limit){
  117.     factors.assign(limit + 5, vector<int>());
  118.  
  119.     for(int i = 1; i <= limit; i++){
  120.         for(int j = i; j <= limit; j += i)
  121.             factors[j].emplace_back(i);
  122.     }
  123. }
  124.  
  125. void solve(){
  126.     int n; cin >> n;
  127.  
  128.     dsu d(n);
  129.  
  130.     vector<int> vec;
  131.     int mx = imin;
  132.     for(int i = 0; i < n; i++){
  133.         int x; cin >> x;
  134.  
  135.         vec.emplace_back(x);
  136.         updmax(mx, x);
  137.     }
  138.     vector<vector<int>> adj(mx + 5, vector<int>());
  139.  
  140.     for(int i = 0; i < n; i++){
  141.         for(auto& factor : factors[vec[i]])
  142.             adj[factor].emplace_back(i);
  143.     }
  144.  
  145.     ll ans = 0;
  146.     for(int i = (int)adj.size() - 1; i >= 1; i--){
  147.         for(int j = 1; j < adj[i].size(); j++){
  148.             int u = adj[i][j], v = adj[i][j - 1];
  149.             if(d.get(u) != d.get(v)){
  150.                 d.join(u, v);
  151.                 ans += i;
  152.             }
  153.         }
  154.     }
  155.     cout << ans;
  156. }
  157.  
  158. void set_file(string &file_name) {
  159.     freopen((file_name + ".in").c_str(), "r", stdin);
  160. //    freopen((file_name + ".out").c_str(), "w", stdout);
  161. }
  162.  
  163. int main() {
  164.     Start_Crushing();
  165.  
  166.     string file_name = "dream";
  167.     set_file(file_name);
  168.  
  169.     build(100000);
  170.  
  171.     int t = 1;
  172.     /*Multiple test cases?*/ cin >> t;
  173. //    while (t--) {
  174. //        solve();
  175. //        if (!t)
  176. //            break;
  177. //        cout << nl;
  178. //    }
  179.  
  180.     for(int tc = 1; tc <= t; tc++){
  181.         cout << "Case " << tc << ": ";
  182.         solve();
  183. //        if(tc != t)
  184.             cout << nl;
  185.     }
  186.  
  187.     return 0;
  188. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement