Friday, March 27, 2020

[Codeforces] E. Subset Sums

Author            : Dipu Kumar Mohanto 
                    CSE, Batch - 6
                    BRUR.
Problem Statement : E. Subset Sums
Source            : Codeforces
Category          : Data Structure
Algorithm         : Heavy Light Technique
Verdict           : Accepted

  1. #include "bits/stdc++.h"  
  2.   
  3. using namespace std;  
  4.   
  5. #define ll           long long int   
  6. #define endl         '\n'  
  7.   
  8. static const int maxn = 1e5 + 5;  
  9. static const int maxb = 320;  
  10.   
  11. vector <int> light_set[maxn];  
  12. vector <int> heavy_set[maxb];  
  13. int light_heavy_intersection[maxn][maxb];  
  14. int heavy_heavy_intersection[maxb][maxb];  
  15. int num_ele[maxn];  
  16. int mapper[maxn];  
  17. long long lazy[maxn];  
  18. long long ans_heavy[maxn];  
  19. long long arr[maxn];  
  20.   
  21. signed main()  
  22. {  
  23.     ios_base::sync_with_stdio(false);  
  24.     cin.tie(nullptr);  
  25.       
  26.     int n, m, q;  
  27.     cin >> n >> m >> q;  
  28.     for (int i = 1; i <= n; i++) cin >> arr[i];  
  29.     int light = 0;  
  30.     int heavy = 0;  
  31.     for (int i = 1; i <= m; i++)  
  32.     {  
  33.         cin >> num_ele[i];  
  34.         if (num_ele[i] > maxb)  
  35.         {  
  36.             mapper[i] = ++heavy;  
  37.             long long sum = 0;  
  38.             heavy_set[heavy].resize(num_ele[i]);  
  39.             for (int &x : heavy_set[heavy])   
  40.             {  
  41.                 cin >> x;  
  42.                 sum += arr[x];        
  43.             }  
  44.             sort(heavy_set[heavy].begin(), heavy_set[heavy].end());  
  45.             ans_heavy[heavy] = sum;  
  46.         }  
  47.         else   
  48.         {  
  49.             mapper[i] = ++light;  
  50.             light_set[light].resize(num_ele[i]);  
  51.             for (int &x : light_set[light]) cin >> x;  
  52.             sort(light_set[light].begin(), light_set[light].end());  
  53.         }  
  54.     }  
  55.     auto get = [&](vector <int> &vec, int key)  
  56.     {  
  57.         auto fnd = lower_bound(vec.begin(), vec.end(), key);  
  58.         return fnd != vec.end() && *fnd == key;   
  59.     };  
  60.     for (int lgt = 1; lgt <= light; lgt++)  
  61.     {  
  62.         for (int x : light_set[lgt])  
  63.         {  
  64.             for (int hvy = 1; hvy <= heavy; hvy++)   
  65.                 light_heavy_intersection[lgt][hvy] += get(heavy_set[hvy], x);  
  66.         }  
  67.     }  
  68.     for (int hvy = 1; hvy <= heavy; hvy++)  
  69.     {  
  70.         for (int x : heavy_set[hvy])  
  71.         {  
  72.             for (int i = 1; i <= heavy; i++) if (i != hvy)  
  73.                 heavy_heavy_intersection[hvy][i] += get(heavy_set[i], x);    
  74.         }  
  75.     }  
  76.     while (q--)  
  77.     {  
  78.         char type;  
  79.         cin >> type;  
  80.         if (type == '?')  
  81.         {  
  82.             int x;  
  83.             cin >> x;  
  84.             if (num_ele[x] > maxb)  
  85.             {  
  86.                 int hvy = mapper[x];  
  87.                 long long res = ans_heavy[hvy] + 1LL * num_ele[x] * lazy[hvy];  
  88.                 for (int i = 1; i <= heavy; i++) if (hvy != i)  
  89.                     res += heavy_heavy_intersection[hvy][i] * lazy[i];  
  90.                 cout << res << endl;  
  91.             }  
  92.             else  
  93.             {  
  94.                 int lgt = mapper[x];  
  95.                 long long res = 0;  
  96.                 for (int p : light_set[lgt]) res += arr[p];  
  97.                 for (int z = 1; z <= heavy; z++) res += (light_heavy_intersection[lgt][z] * lazy[z]);  
  98.                 cout << res << endl;  
  99.             }  
  100.         }  
  101.         else   
  102.         {  
  103.             int x;  
  104.             long long val;  
  105.             cin >> x >> val;  
  106.             if (num_ele[x] > maxb)   
  107.             {  
  108.                 int hvy = mapper[x];  
  109.                 lazy[hvy] += val;  
  110.             }  
  111.             else   
  112.             {  
  113.                 int lgt = mapper[x];  
  114.                 for (int p : light_set[lgt]) arr[p] += val;  
  115.                 for (int hvy = 1; hvy <= heavy; hvy++) ans_heavy[hvy] += (light_heavy_intersection[lgt][hvy] * val);  
  116.             }  
  117.         }  
  118.     }  
  119. }  

No comments:

Post a Comment

Note: Only a member of this blog may post a comment.