-
Notifications
You must be signed in to change notification settings - Fork 59
/
Copy pathevantai.cpp
61 lines (54 loc) · 1.31 KB
/
evantai.cpp
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
# include <fstream>
# include <algorithm>
# define LB(p) ((-p)&p)
# define NR 1005
# define MOD 30103
using namespace std;
ifstream f("evantai.in");
ofstream g("evantai.out");
struct elem {
int ci, cs, suma;
}v[NR*NR];
bool cmp (elem x, elem y) {
if (x.suma > y.suma) return 0;
else if (x.suma==y.suma && x.cs - x.ci < y.cs - y.ci) return 0;
else return 1;
}
int i,j,n,m,x,y,nr,S,sol;
int AIB[NR][NR], a[NR];
int query (int X, int Y) {
int sol=0;
for (int i=X; i>=1; i-=LB(i))
for (int j=Y; j>=1; j-=LB(j)) {
sol=sol + AIB[i][j];
if (sol >= MOD) sol-=MOD;
}
return sol;
}
void update (int X, int Y, int suma) {
for (int i=X; i<=1000; i+=LB(i))
for (int j=Y; j<=1000; j+=LB(j)) {
AIB[i][j]+=suma;
if (AIB[i][j] >= MOD) AIB[i][j]-=MOD;
}
}
int main ()
{
f>>n;
for (i=1; i<=n; ++i)
f>>a[i];
for (i=1; i<=n; ++i)
for (j=i+1; j<=n; ++j) {
++nr;
v[nr].ci=i; v[nr].cs=j; v[nr].suma=a[i]+a[j];
}
sort (v+1, v+nr+1, cmp);
for (i=1; i<=nr; ++i) {
S=(query (v[i].cs-1, v[i].cs-1) - query (v[i].ci, v[i].cs-1) + MOD)%MOD;
sol=sol + S + 1;
if (sol >= MOD) sol-=MOD;
update (v[i].ci, v[i].cs, S+1);
}
g<<sol<<"\n";
return 0;
}