#pragma GCC optimize ("Ofast")
#define _USE_MATH_DEFINES
#include <bits/stdc++.h>
#define FOR(i, a, b) for (auto i=(a); i<(b); i++)
#define FORD(i, a, b) for (int i=(a); i>(b); i--)
#define SZ(x) ((int)(x).size())
#define ALL(x) (x).begin(), (x).end()
#define PPC(x) __builtin_popcountll(x)
#define MSB(x) (63 - __builtin_clzll(x))
#define LSB(x) __builtin_ctzll(x)
#define ARG(x, i) (get<i>(x))
#define ithBit(m, i) ((m) >> (i) & 1)
#define pb push_back
#define ft first
#define sd second
#define kw(a) ((a) * (a))
#define CLR(x) x.clear(), x.shrink_to_fit()
#ifdef DEBUG
#include "debug.h"
#else
#define dbg(...) 0
#endif
using namespace std;
template <typename T1, typename T2> inline void remin(T1& a, T2 b) { a = min(a, (T1)b); }
template <typename T1, typename T2> inline void remax(T1& a, T2 b) { a = max(a, (T1)b); }
const int maxN = 1 << 9;
map <int, int> pos, neg;
int T[maxN];
long long npo2(long long n)
{
return n * (n-1) / 2;
}
long long npo3(long long n)
{
return n * (n-1) * (n-2) / 6;
}
long long findCustom(const auto& M, int k)
{
auto it = M.find(k);
return it == M.end()
? 0
: it->sd;
}
long long f(const auto& A, const auto& B)
{
long long ret = 0;
for (auto it = A.begin(); it != A.end(); it++)
{
for (auto jt = A.begin(); jt != it; jt++)
ret += findCustom(B, it->ft + jt->ft) * it->sd * jt->sd;
ret += findCustom(B, it->ft * 2) * npo2(it->sd);
}
return ret;
}
void solve()
{
int n;
scanf ("%d", &n);
FOR(i, 1, n+1)
scanf ("%d", T+i);
int zeros = 0;
FOR(i, 1, n+1)
{
int s = 0;
FOR(j, i, n+1)
{
s += T[j];
if (s < 0)
neg[-s]++;
if (s == 0)
zeros++;
if (s > 0)
pos[s]++;
}
}
long long res = npo3(zeros);
for (auto& [s, cnt] : pos)
res += findCustom(neg, s) * zeros * cnt;
res += f(pos, neg);
res += f(neg, pos);
printf("%lld\n", res);
// check(n, res);
}
int main()
{
int t = 1;
// scanf ("%d", &t);
FOR(tid, 1, t+1)
{
//printf("Case #%d: ", tid);
solve();
}
return 0;
}
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 | #pragma GCC optimize ("Ofast") #define _USE_MATH_DEFINES #include <bits/stdc++.h> #define FOR(i, a, b) for (auto i=(a); i<(b); i++) #define FORD(i, a, b) for (int i=(a); i>(b); i--) #define SZ(x) ((int)(x).size()) #define ALL(x) (x).begin(), (x).end() #define PPC(x) __builtin_popcountll(x) #define MSB(x) (63 - __builtin_clzll(x)) #define LSB(x) __builtin_ctzll(x) #define ARG(x, i) (get<i>(x)) #define ithBit(m, i) ((m) >> (i) & 1) #define pb push_back #define ft first #define sd second #define kw(a) ((a) * (a)) #define CLR(x) x.clear(), x.shrink_to_fit() #ifdef DEBUG #include "debug.h" #else #define dbg(...) 0 #endif using namespace std; template <typename T1, typename T2> inline void remin(T1& a, T2 b) { a = min(a, (T1)b); } template <typename T1, typename T2> inline void remax(T1& a, T2 b) { a = max(a, (T1)b); } const int maxN = 1 << 9; map <int, int> pos, neg; int T[maxN]; long long npo2(long long n) { return n * (n-1) / 2; } long long npo3(long long n) { return n * (n-1) * (n-2) / 6; } long long findCustom(const auto& M, int k) { auto it = M.find(k); return it == M.end() ? 0 : it->sd; } long long f(const auto& A, const auto& B) { long long ret = 0; for (auto it = A.begin(); it != A.end(); it++) { for (auto jt = A.begin(); jt != it; jt++) ret += findCustom(B, it->ft + jt->ft) * it->sd * jt->sd; ret += findCustom(B, it->ft * 2) * npo2(it->sd); } return ret; } void solve() { int n; scanf ("%d", &n); FOR(i, 1, n+1) scanf ("%d", T+i); int zeros = 0; FOR(i, 1, n+1) { int s = 0; FOR(j, i, n+1) { s += T[j]; if (s < 0) neg[-s]++; if (s == 0) zeros++; if (s > 0) pos[s]++; } } long long res = npo3(zeros); for (auto& [s, cnt] : pos) res += findCustom(neg, s) * zeros * cnt; res += f(pos, neg); res += f(neg, pos); printf("%lld\n", res); // check(n, res); } int main() { int t = 1; // scanf ("%d", &t); FOR(tid, 1, t+1) { //printf("Case #%d: ", tid); solve(); } return 0; } |
English