Editorial of Codeforces Round #672 (Div. 2)

#	User	Rating
1	tourist	3985
2	jiangly	3741
3	jqdai0815	3682
4	Benq	3529
5	orzdevinwang	3526
6	ksun48	3489
7	Radewoosh	3483
8	Kevin114514	3442
9	ecnerwala	3392
9	Um_nik	3392

#	User	Contrib.
1	cry	169
2	maomao90	162
2	Um_nik	162
2	atcoder_official	162
5	djm03178	158
6	-is-this-fft-	157
7	adamant	155
8	awoo	154
8	Dominater069	154
10	nor	150

Hint

Solution

Code in C++ (Wind_Eagle)

#include<iostream>

using namespace std;

int a[1000000+5];

int main()
{
    ios_base::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    int t;
    cin>>t;
    while (t--)
    {
        int n;
        cin>>n;
        for (int i=0; i<n; i++)
        {
            cin>>a[i];
        }
        bool can=false;
        for (int i=1; i<n; i++)
        {
            if (a[i]>=a[i-1])
            {
                can=true;
                break;
            }
        }
        if (can) cout<<"YES"<<'\n';
        else cout<<"NO"<<'\n';
    }
}

1420B - Rock and Lever

Hint

Solution

Code in C++ (Wind_Eagle)

#include<iostream>
#include<vector>
#include<algorithm>
#include<ctime>
#include<random>

using namespace std;

mt19937 rnd(time(NULL));

int a[1000000+5];

int main()
{
    ios_base::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    int t;
    cin>>t;
    while (t--)
    {
        int n;
        cin>>n;
        for (int i=0; i<n; i++)
        {
            cin>>a[i];
        }
        int64_t ans=0;
        for (int j=29; j>=0; j--)
        {
            int64_t cnt=0;
            for (int i=0; i<n; i++)
            {
                if (a[i]>=(1<<j)&&a[i]<(1<<(j+1)))
                {
                    cnt++;
                }
            }
            ans+=cnt*(cnt-1)/2;
        }
        cout<<ans<<'\n';
    }
}

1420C1 - Pokémon Army (easy version)

Hint

Solution

Code in C++ (gepardo)

#include <iostream>
#include <vector>

using namespace std;

inline int64_t calc(const vector<int> &a) {
	int n = a.size();
	vector<int64_t> d1(n+1), d2(n+1);
	d1[0] = -static_cast<int64_t>(1e18);
	d2[0] = 0;
	for (int i = 0; i < n; ++i) {
		d1[i+1] = max(d1[i], d2[i] + a[i]);
		d2[i+1] = max(d2[i], d1[i] - a[i]);
	}
	return max(d1.back(), d2.back());
}

int main() {
	ios_base::sync_with_stdio(false);
	int t; cin >> t;
	while (t--) {
		int n, q; cin >> n >> q;
		vector<int> a(n);
		for (int i = 0; i < n; ++i) {
			cin >> a[i];
		}
		cout << calc(a) << "\n";
		for (int i = 0; i < q; ++i) {
			int l, r; cin >> l >> r; --l; --r;
			swap(a[l], a[r]);
			cout << calc(a) << "\n";
		}
	}
	return 0;
}

1420C2 - Pokémon Army (hard version)

Hint

Solution

Code in Java (gepardo)

import java.io.*;
import java.util.*;

public class C2_Java {
    int[] a;
    int n;
    long ans;
    
    void add(int i, int sign) {
        if (i == 0 || i == n+1) return;
        if (a[i-1] < a[i] && a[i] > a[i+1]) ans += sign * a[i];
        if (a[i-1] > a[i] && a[i] < a[i+1]) ans -= sign * a[i];
    }
    
    void addTwo(int l, int r, int sign) {
        if (l == r) return;
        if (l+1 == r) {
            add(l-1, sign);
            add(l, sign);
            add(r, sign);
            add(r+1, sign);
            return;
        }
        add(l-1, sign);
        add(l, sign);
        add(l+1, sign);
        if (l+2 != r) add(r-1, sign);
        add(r, sign);
        add(r+1, sign);
    }
    
    public void run() throws IOException {
        BufferedReader in = new BufferedReader(new InputStreamReader(System.in));
        BufferedWriter out = new BufferedWriter(new OutputStreamWriter(System.out));
        
        StringTokenizer tok = new StringTokenizer(in.readLine());
        int t = Integer.parseInt(tok.nextToken());
        while (t --> 0) {
            tok = new StringTokenizer(in.readLine());
            n = Integer.parseInt(tok.nextToken());
            int q = Integer.parseInt(tok.nextToken());
            a = new int[n + 2];
            tok = new StringTokenizer(in.readLine());
            for (int i = 1; i <= n; ++i) {
                a[i] = Integer.parseInt(tok.nextToken());
            }
            a[0] = -1;
            a[n+1] = -1;
            ans = 0;
            StringBuilder builder = new StringBuilder();
            for (int i = 1; i <= n; ++i) {
                if (a[i-1] < a[i] && a[i] > a[i+1]) ans += a[i];
                if (a[i-1] > a[i] && a[i] < a[i+1]) ans -= a[i];
            }
            builder.append(ans).append("\n");
            while (q --> 0) {
                tok = new StringTokenizer(in.readLine());
                int l = Integer.parseInt(tok.nextToken());
                int r = Integer.parseInt(tok.nextToken());
                if (l == r) {
                    builder.append(ans).append("\n");
                    continue;
                }
                addTwo(l, r, -1);
                int tmp = a[l];
                a[l] = a[r];
                a[r] = tmp;
                addTwo(l, r, +1);
                builder.append(ans).append("\n");
            }
            out.write(builder.toString());   
        }
        
        in.close();
        out.close();        
    }

    public static void main(String[] args) throws IOException {
        (new C2_Java()).run();
    }
}

Code in C++ (Wind_Eagle)

#include<iostream>

using namespace std;

#define int long long

int a[1000000+5];

int n;

int ans=0;

inline void insert(int i)
{
    if (i==0||i==n+1) return;
    if (a[i-1]<a[i]&&a[i]>a[i+1]) ans+=a[i];
    if (a[i-1]>a[i]&&a[i]<a[i+1]) ans-=a[i];
}

inline void erase(int i)
{
    if (i==0||i==n+1) return;
    if (a[i-1]<a[i]&&a[i]>a[i+1]) ans-=a[i];
    if (a[i-1]>a[i]&&a[i]<a[i+1]) ans+=a[i];
}

int32_t main()
{
    ios_base::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    int t;
    cin>>t;
    while (t--)
    {
        int q;
        cin>>n>>q;
        for (int i=1; i<=n; i++)
        {
            cin>>a[i];
        }
        a[0]=-1;
        a[n+1]=-1;
        ans=0;
        for (int i=1; i<=n; i++)
        {
            if (a[i-1]<a[i]&&a[i]>a[i+1]) ans+=a[i];
            if (a[i-1]>a[i]&&a[i]<a[i+1]) ans-=a[i];
        }
        cout<<ans<<'\n';
        while (q--)
        {
            int l,r;
            cin>>l>>r;
            erase(l-1);
            erase(l);
            erase(l+1);
            if (l!=r)
            {
                if (r-1!=l+1&&r-1!=l) erase(r-1);
                if (r!=l+1) erase(r);
                erase(r+1);
            }
            swap(a[l],a[r]);
            insert(l-1);
            insert(l);
            insert(l+1);
            if (l!=r)
            {
                if (r-1!=l+1&&r-1!=l) insert(r-1);
                if (r!=l+1) insert(r);
                insert(r+1);
            }
            cout<<ans<<'\n';
        }
    }
}

1420D - Rescue Nibel!

Hint

Solution

Code in Java (gepardo)

import java.io.*;
import java.util.*;

public class D_Java {
	public static final int MOD = 998244353;
	
	public static int mul(int a, int b) {
		return (int)((long)a * (long)b % MOD);
	}
	
	int[] f;
	int[] rf;
	
	public int C(int n, int k) {
		return (k < 0 || k > n) ? 0 : mul(f[n], mul(rf[n-k], rf[k]));
	}
	
	public static int pow(int a, int n) {
		int res = 1;
		while (n != 0) {
			if ((n & 1) == 1) {
				res = mul(res, a);
			}
			a = mul(a, a);
			n >>= 1;
		}
		return res;
	}
	
	static void shuffleArray(int[] a) {
		Random rnd = new Random();
		for (int i = a.length-1; i > 0; i--) {
			int index = rnd.nextInt(i + 1);
			int tmp = a[index];
			a[index] = a[i];
			a[i] = tmp;
		}
	}
	
	public static int inv(int a) {
		return pow(a, MOD-2);
	}
	
	public void doIt() throws IOException {
		BufferedReader in = new BufferedReader(new InputStreamReader(System.in));
		
		StringTokenizer tok = new StringTokenizer(in.readLine());
		int n = Integer.parseInt(tok.nextToken());
		int k = Integer.parseInt(tok.nextToken());
		
		f = new int[n+42];
		rf = new int[n+42];
		f[0] = rf[0] = 1;
		for (int i = 1; i < f.length; ++i) {
			f[i] = mul(f[i-1], i);
			rf[i] = mul(rf[i-1], inv(i));
		}
		
		int[] events = new int[2*n];
		for (int i = 0; i < n; ++i) {
			tok = new StringTokenizer(in.readLine());
			int le = Integer.parseInt(tok.nextToken());
			int ri = Integer.parseInt(tok.nextToken());
			events[i] = le*2;
			events[i + n] = ri*2 + 1;
		}
		shuffleArray(events);
		Arrays.sort(events);
		
		int ans = 0;
		int balance = 0;
		for (int r = 0; r < 2*n;) {
			int l = r;
			while (r < 2*n && events[l] == events[r]) {
				++r;
			}
			int added = r - l;
			if (events[l] % 2 == 0) {
				// Open event
				ans += C(balance + added, k);
				if (ans >= MOD) ans -= MOD;
				ans += MOD - C(balance, k);
				if (ans >= MOD) ans -= MOD;
				balance += added;
			} else {
				// Close event
				balance -= added;
			}
		}
		
		in.close();
		System.out.println(ans);
	}
	
	public static void main(String[] args) throws IOException {
		(new D_Java()).doIt();
	}
}

1420E - Battle Lemmings

Hint 1

Hint 2

Hint 3

Solution

Code in C++ (gepardo)

#include <cassert>
#include <climits>
#include <iostream>
#include <vector>
#include <numeric>

using namespace std;

template<typename T>
inline void umin(T &a, const T &b) {
	a = min(a, b);
}

template<typename T>
inline void umax(T &a, const T &b) {
	a = max(a, b);
}

int main() {
	ios_base::sync_with_stdio(false);
	
	int n; cin >> n;
	vector<int> a(n);
	for (int &x : a) {
		cin >> x;
	}
	a.push_back(1);
	
	vector<int> gg;
	for (int i = 0; i <= n; ++i) {
		int s = i;
		while (a[i] == 0) ++i;
		gg.push_back(i - s);
	}
	
	vector<int> p = gg;
	partial_sum(begin(p), end(p), begin(p));
	
	constexpr int mx = 103;
	constexpr int dx = mx + 1, dy = mx + 1, dz = mx * (mx + 1) / 2 + 1;
	static int dp[dx][dy][dz] = {};
	
	for (int i = 0; i < dx; ++i) {
		for (int j = 0; j < dy; ++j) {
			for (int k = 0; k < dz; ++k) {
				dp[i][j][k] = INT_MAX;
			}
		}
	}
	dp[0][0][0] = 0;
	
	int k = gg.size(), l = p.back();
	for (int i = 0; i < k; ++i) {
		for (int j = 0; j <= l; ++j) {
			for (int s = 0; s <= n * (n-1) / 2; ++s) {
				if (dp[i][j][s] == INT_MAX) continue;
				for (int q = j; q <= l; ++q) {
					umin(dp[i+1][q][s + abs(q - p[i])], dp[i][j][s] + (q-j)*(q-j));
				}
			}
		}
	}
	
	int mn = INT_MAX;
	for (int s = 0; s <= n * (n-1) / 2; ++s) {
		umin(mn, dp[k][l][s]);
		int val = l*l - mn;
		assert(val % 2 == 0);
		cout << val/2 << " ";
	}
	cout << endl;
	
	return 0;
}

#include<iostream> #include<vector> #include<algorithm> #include<ctime> #include<random> using namespace std; mt19937 rnd(time(NULL)); int a[1000000+5]; int main() { ios_base::sync_with_stdio(false); cin.tie(0); cout.tie(0); int t; cin>>t; while (t--) { int n; cin>>n; for (int i=0; i<n; i++) { cin>>a[i]; } int64_t ans=0; for (int j=29; j>=0; j--) { int64_t cnt=0; for (int i=0; i<n; i++) { if (a[i]>=(1<<j)&&a[i]<(1<<(j+1))) { cnt++; } } ans+=cnt*(cnt-1)/2; } cout<<ans<<'\n'; } }

n! = n*(n-1)! (n!)^(mod-2)%mod = ((n * (n-1)!)^(mod-2))%mod = (n ^(mod -2) * (n-1)!^(mod-2))%mod (invf(n)*n)%mod = (n!^(mod-2) * n)%mod = (n^(mod-1) * (n-1)!^(mod-2))%mod = 1 * (n-1)!^(mod-2)%mod = invf(n-1)%mod

vector<pair<int, int> > v; int L[n], R[n]; int ans = 0, cnt = 0; for (int i = 0; i < n; i++) { cin >> L[i] >> R[i]; v.push_back({L[i], 1}); // starting time of lamp v.push_back({R[i] + 1, -1}); //ending time + 1 of lamp } sort(v.begin(), v.end()); for (int i = 0; i < 2 * n; i++) { int S = v[i].second; if (S == 1) { if (cnt >= k - 1) { ans += NCR(cnt, k - 1, MOD); ans = ans % MOD; } cnt++; } else { cnt--; } } cout << ans;

/* * For every endPoint, add the number of ways to choose k segments * Such that the segment(call it p) with that endpoint is one of the k segments * Which basically leaves you to choose the k-1 segments that pass through * that endpoint, other than p. */

Comments (203)

Show archived | Write comment?

gepardo

4 years ago, # |

Auto comment: topic has been translated by gepardo (original revision, translated revision, compare)

→ Reply

Wind_Eagle

← Rev. 2 →

+172

Actually, many people say that now there are practically no data structure or DP tasks, only ad-hoc and constructive. Well, this contest disproves this. There were no constructive tasks, but there were Segment Tree in C2 (there is a solution using it) and DP in E.

tfg

4 years ago, # ^ |

+245

One contest doesn't disprove a trend.

kcravuritg

I am actually better at ad-hoc and constructive currently and loved the trend :sobs:, But I do like these rounds where I learn a lot.

zukonit14

Please could you explain/hint about the Segment tree solution of C2. How will we merge answer of left and right for getting answer of current node ? Thanks in advance!

_Bishop_

+19

You can maintain four quantities in the segment tree node i.e maximum strength given that
1. first element has $$$plus$$$ sign before it and last one also has $$$plus$$$ sign
2. first element has $$$plus$$$ sign and last has $$$minus$$$
3. first element has $$$minus$$$ sign and last has $$$plus$$$ sign
4. first element has $$$minus$$$ sign and last has $$$minus$$$ sign
You can easily implement merge and update operations.
Here is my implementation https://mirror.codeforces.com/contest/1420/submission/93716090. There may be easier solutions but this one was more obvious and intutive

Got it bro. Thanks a lot _Bishop_.

snehal007

+18

Tch tch, "fenwick123" had to use segment tree. xD

The_mysterio

Hi Sir, I wanted to know is there any site that focusses more on data structures and less on ad-hoc?

Vedkribhu

I think atcoder contests are better in this matter.

spashal

I implemented segtree and received tle in C2. someone willing to help? https://mirror.codeforces.com/contest/1420/submission/93731337

Aggu_01000101

Massive output, use \n instead of endl.

Also, if you exepected segtree, the constraints could've been relaxed (:

Mohsina_Shaikh

2 years ago, # ^ |

Wind_Eagle Can you please Explain 1420B.Rock and Lever Not able to understand.From the editorial code

Why to need to check a[i]<(1<<j+1).Why the total number of pairs is cnt*(cnt-1)/2.Can you please explain.Tried to understand for so long still not clear

+12

a[i]>=(1<<j) is actually $$$a_i \geq 2^j$$$, and a[i]<(1<<(j+1)) is actually $$$a_i < 2^{j+1}$$$. If some number is greater or equal than $$$2^j$$$, and strictly smaller, than $$$2^{j+1}$$$, than it contains $$$j$$$-th bit as the greatest bit.

cnt*(cnt-1)/2 is actually $$$C_n^2$$$, where C is binomial coefficient. More easily, this is the number of ways to pick to different things from a set of $$$n$$$ things. In this task we want to choose the number of ways to pick two numbers from $$$n$$$ numbers.

It seems to me that you should first study some simple maths, like basic combinatorics, bit arithmetics and number theory. With mathematical base it is much easier to study CP.

Thank you for solving my contest!

boleyn.su

+62

For E, if we accept solution in O(n^5), I think it may be better to choose a smaller n. For me, I was too afraid that an O(n^5) solution will TL and just tried to think of a better solution.

+26

We initially thought of making $$$n \le 100$$$ in this problem, because my solution in $$$O(n^5)$$$ passes with such constraints. Anyway, the DP solution is pretty fast and $$$n^5$$$ can be easily divided by some constant factor, so the constraints may seem too big for $$$O(n^5)$$$.

+34

I can understand your choosing n<=80. Maybe I just need to accept the fact that computers are much faster nowadays.

BTW, my solution at https://paste.storage.boleyn.su/cfr672d2e.cpp should have a better time complexity. I have not analyzed its TC yet so this is just my guess.

Can you explain your approach? Is it true that you use a DP similar to the one described in the editorial, but you optimized the code by skipping invalid DP states?

← Rev. 4 →

For the last dimension, it is non-decreasing, so we can express it using a vector<pair<int,int>> and reduce some useless states.

My guess is that there may be a bound better than O(n^2) for it. It would be nice to prove or disprove my guess.

UPD: I tested it with larger constraints. It turns out that my guess is wrong. The size of vector<pair<int,int>> is still O(n^2).

rahulmysuru7

I felt like D had hard time limits. But, later realized after contest with precomputation of inverse of each factorial it could have been made faster.

chaluMan

Can you(someone) explain how to do that? Cause I got a TLE on pretest-11.

precompute fact 1 to n and mod inverse of those instead of computing them every time.

-6

Actually, I did that but still got a TLE(UPDATED: NEVER MIND I PASSED IT)

codephilic

you can calculate the mod inverse at run time also, pre-computing the factorials is enough

grimjaw

Precompute factorials from 0 to N (here $$$3*10^{5}) $$$ initially in an array, say $$$fact$$$.
Then calculate inverse factorial of N using formula $$$inv_f(N) = fact(N)^{MOD-2} \% MOD $$$
which takes $$$ O(\log_2 (MOD))$$$ time. Using this value, now fill the remaing values from N-1 to 0 using recurrence : $$$inv_f(i) = (i+1)*inv_f(i+1)$$$.
Thus we can precompute inverse factorials in $$$ O(N) $$$ complexity.
Solution for reference : 93729779

endless_summer

Hey, I precomputed the factorials and inverse factorials similar to your solution but I am getting TLE.

Use binary exponentiation in your powmod function to calculate $$$a^b$$$ in log(b) complexity, you are calculating it in O(b) time. You can read it here : cp-algo

Thanks, didn't know about binary exponentiation.

spookywooky

Note that precalculateing the inverse factorial does not give a better runtime. We never need to calculate more than $$$3*10^5$$$ inverses. So precalculating them just uses more memory, but not less time.

But to calculate $$$3*10^5$$$ inverse factorials you'd have to calculate them in $$$O(log_2(Mod))$$$ individually resulting in $$$O(N*log(Mod)$$$ time. But here I am precalculating them in $$$O(N)$$$, so I am unable to understand your logic behind why this won't give better runtime? Can you please elaborate.

You are right. I was irritated by the long runtime of your solution. But that is caused by the use of the map instead of sorting an array once.

kgargdun

I encountered this for the first time, tried to dig around but can't figure out why this works-

invf[i] = invf[i+1]*(i+1)%mod;

Can u help or tell me where to look.

← Rev. 3 →

Does this makes sense?

$$$(n+1)!^{−1} = (n!*(n+1))^{-1}$$$
$$$(n+1)!^{−1} = n!^{−1}*(n+1)^{−1}$$$
$$$(n+1)!^{−1}*(n+1) = n!^{−1}$$$

Watch from 16:00 onwards Lecture on Combinatorics By Kevin Charles Atienza

sasha4

Hi there. Please, tell me about inverse factorial of N. What is it?

aniervs

I precomputed the factorials and the inverses, thus my solution was linear, though I think $$$O(n\log mod)$$$ should fit in time.

Ebiarat

+53

Great problemset, like it . waiting for your next contest,thanks

sid9406

For some problem editorial is in english and for others it's in russian . Pls Fix.

The English versions will be loaded in several minutes.

shubham7999

Can you tell if the elements of array in C1 are not distinct what will happen??

I don't understand your question. The elements $$$a_i$$$ are distinct by problem statement. Or do you want to know if the intended solution works if equal elements are allowed?

Yess if we go by approach for maxima and minima approach for C1 problem then will it work??

It is unclear how we will treat an element if its neighbor is equal to it. Or even worse, when both neighbors are equal to our element. Will it be a local maximum or local minimum? Also it's important that local minima and local maxima must alternate (this property is used actively in the proof).

Anyway, it's an interesting question on how to fix the solution when equal elements are allowed.

If you solve C1 using dynamic programming, you should not be affected by equal elements.

tourist

+56

One may assume that out of equal neighboring elements, the left one is smaller. This is equivalent to adding $$$i \cdot \varepsilon$$$ to $$$a_i$$$, which makes all elements distinct while keeping the answer the same if $$$\varepsilon \rightarrow 0$$$.

clyring

+20

A trivial fix is to break ties by index, so that if $$$a_i = a_{i+1}$$$, we can simply pretend $$$a_i < a_{i+1}$$$. This is safe because the value of the optimal solution is continuous in the values of $$$a_i$$$ and we can imagine taking $$$a_i = \lim_{\varepsilon \to 0^+} a_i + i \cdot \varepsilon$$$ simultaneously at all indices. Alternatively, the solutions based on the $$$\sum_{i = 1}^n \max (0, a_i - a_{i-1})$$$ formula or on segment trees are completely unchanged, since they do not make use of the distinctness constraint in any way.

EDIT: I was beaten to the punch! So it goes.

codepasta

what is "Разбор"???? is it rated??

It is a Russian word for "parsing".

+21

It's better to translate "разбор" as "analysis" :)

meet2mky

Round was great. Kudos to author!!

.Mr.

My editorial is in Russian :(

Everything is OK now, except problem E.

UPD: Problem E is also OK now.

blobugh

Does "except problem E" mean this?

Yes.

Prateek_Jakhar

I could not solve a single problem in this contest, should I die? Why am I even alive.

Greedy_Optimizzer

+14

Don't feel demotivated you have given just few contest.just keep practicing and keep motivated :)

ustaritz

just take your meds : one tourist stream every day ! starting with : https://www.youtube.com/watch?v=97tieEKfvBs

streaker_rules

I don't understand the solution for problem D. This is my understanding so far :-

For example: if there is a starting point on 42 and our count of segments excluding this is 10 and the value of k is 3. then I found out the solution is to add ncr(10,k-1).Now the count is 11 and if we again reach to another starting position on 45 so again we will add ncr(11,k-1).Now count is 12.

So if I understood it correctly, then we will choose one out of k to be the current one and choose rest k-1 from the previous segments right?

If we fixed a starting point, we may pick any $$$k$$$ segments that contain it. But there must be at least one segment that start in this point (note that more than one one segment may start in the same point).

Actually, I was talking something like this solution.

I tried this concept after the contest.

I understood your solution. It's mostly the same as in the tutorial, but it considers the segments that start in the same point in a shorter and nicer way.

PurpleThinker

+15

How to solve C2 using segment tree?

Jesuscristo

+27

I used an ad-hoc segment tree with 4 values for each range:

The max value you can get by using an even number of elements in that range
The max value for an odd number of elements
The min value for an even number of elements
The min value for an odd number of elements

The value (max/even) of a node is the max from these values:

left child's (max/even) + right child's (max/even)
left child's (max/odd) – right child's (min/odd)

The rest of the values are obtained similarly. Updating the tree is the same as a regular segment tree.

Here is my code.

DumbCarlsen69

I am still a bit surprised to see the no. of submissions on the Problem D, Is it some sort of a standard algo or question? Because some were able to this in a matter of just few minutes, if it is then please do tell me

Prajwal_

I don't know weather it's standard or not but many interval related questions have similar approach.

decoder__

What is the approach?

My approach for D. Rescue Nibel : I will tell you the method i used to solve, i don't know weather its similar to tutorial's but mine one is easy.I will create a array and put inside L and R + 1 for all intervals. So array will have size = 2 * N. Then sort the array. We keep a variable cnt which will tell the number of lamps which are currently ON. Initially cnt = 0. Now we will iterate over all element of array (in chronological order). Let the current element be x. Now x will either be starting time or (ending time + 1) of some interval. If it is starting time and cnt >= k - 1 then we add NCR(cnt, k - 1) to the answer. We take k - 1 because the kth element will be the lamp whose starting time is x and remaining k - 1 we take from our cnt(lamps which are already in ON state). Then we increment our cnt by 1. If x is (ending time + 1) of some lamp then we just decrement our cnt. This is my code --->

harshitsoora

Can you please share the code for NCR(cnt, k, MOD)

I did follow the same approach but end up having a run time error which was certainly "Arithmetic Exception" while computing this Combination

It fails over 11th test case....could you please share some insights

My Submission

Please read this https://www.geeksforgeeks.org/compute-ncr-p-set-3-using-fermat-little-theorem/

Mooncrater

For anyone looking for his submission.

I think you should try ** Min Stations for trains ** problem on GFG.

SkroX_5

https://www.geeksforgeeks.org/number-of-ways-to-choose-k-intersecting-line-segments-on-x-axis/ I found it on geeks for geeks after the contest.

robinz62

I had a Java solution TLE on D which I fixed by just using a hashmap to cache the results of modular inverse :/

w/o cache: https://mirror.codeforces.com/contest/1420/submission/93688984

w/ cache: https://mirror.codeforces.com/contest/1420/submission/93691811

RisingWizard

Solution for Problem C2

Suddenly it turns out that a single request will change the state of no more than six elements

Didn't get intuition for this statement! Someone Please elaborate

+24

Only the neighbour elements affect local maximum and local minimum. So, if you swap $$$l$$$ and $$$r$$$, you will need to recalculate only $$$l$$$, $$$r$$$ and their neighbours (in the worst case, it will be $$$l-1$$$, $$$l$$$, $$$l+1$$$, $$$r-1$$$, $$$r$$$ and $$$r+1$$$, thus six elements).

pandastic

Whether an element is a local maximum (or minimum) or not solely depends on its relation to the neighbors, as the required condition for a local maximum is a_i > a_{i-1} and a_i > a_{i+1} (reverse signs for a minimum). Now when swapping elements at given positions l and r, this will affect the status of the elements at positions l-1, l, l+1, and r-1, r, r+1 (note that these triplets may partly overlap).

PalForce

https://medium.com/@harryjobz/pok%C3%A9mon-army-3f4348d57143

sharad1103

can't we solve problem b using tries

Venti_chai

+43

Why cut a apple using sword?

H4wk3ye

Yes, we can. https://mirror.codeforces.com/contest/1420/submission/93721816

naeem4265

I am looking details description of solution D or a video tutorial

Isn't the solution of problem D detailed enough?

may be..but i can't understand this..

learner7321

You can try looking at secondthread's youtube channel , I guess it may help.Youtube channel name is "SecondThread".

mehulcoder

This is the implementation. Hope it helps.

Sharon

+30

I've succesfully uphacked myself. Turns out I forgot to check if all ai are different and still passed systests :D

Devised_positive

Never knew that we can hack our own answers!

__StealthX__

Can anyone please explain why my code for Problem D got Accepted when I used fastI/O but was getting TLE without using fastI/O even after that fact that the question had just one testcase per testfile?

Aditya_Sharma

fast I/O is so powerful that it can even reduce your acceptance time to half. Even I was amazed when I first observed it, few days back.

minhcool

-49

It is because of the large input file, you have to read 6 * 10^5 numbers

Thanx a lot buddy!

-47

Your welcome my friend.

r1683n1j23

-63

C1 was already on internet. Make this contest unrated https://tkramesh.wordpress.com/2009/10/08/maximum-alternating-sum-subsequence/

Samarth12

+31

C1 is a standard dp problem, you can't argue on that basis to make the contest unrated. It is difficult to make new problems and very easy to make new accounts.

ACeed

I agree making a new problem is a difficult task but at least we shouldn't use the same language as has been used on other sources. Just googling "maximum alternating sum subsequence" gives you the solution (I googled after seeing this).

MohamedMagdy

Can C1 be solved w/o dynamic programming i.e greedy?

See solution for C2.

goingOn

← Rev. 7 →

can u help me find where am i going wrong in my C1 submission. 93713126
UPD : I didn't use long long .

Nick9

3 years ago, # ^ |

LOL. had the same problem. thank you for your UPD! kinda underestimated the max score that can be achieved.

yes you just need to select each local maxima and minima, and remove the last minima.

lightseba

+11

Is there a typo in the summation of $$$p(A)$$$ for E? Isn't it

$$$\displaystyle \sum_{1 \leq i,j \leq k} f_i \cdot f_j = \left( \sum_{i=1}^k f_i \right)^2 \neq \left( \sum_{i=1}^k f_i^2 \right)^2$$$

Yes, this is a typo. Thanks for finding it out.

StrawHatWess

Segment tree solution for C2

O_E

i think the time limit for D is really bad, i made this submission in the contest and i got a TLE after it ,i don't know why my code gets TLE is the problem with factor 2? https://mirror.codeforces.com/contest/1420/submission/93716363

Same happend with me bro see this, can you help why its happening in my code?

blue_m00n_eoy

I really like this kind of editorial where hints are provided rather than directly displaying the solution approach completely. I wish more than more authors adopt this strategy. Also, the problem statements were quite nice. Thank you for the good round and helpful editorials.

You are welcome :)

wttc

thanks for giving hints before solutions, would expect the same for future contests

ashanksingh1997

I heard that D was copied from gfg. I was surprised at number of submissions of D

An observation free solution for C1 and C2:

First let's try to solve this problem using divide and conquer. On every divide stage, we need the following information for each halve:

$$$1 = $$$ Maximum sum subsequence starting with a positive element and ending with a positive element
$$$2 = $$$ Maximum sum subsequence starting with a positive element and ending with a negative element
$$$3 = $$$ Maximum sum subsequence starting with a negative element and ending with a positive element
$$$4 = $$$ Maximum sum subsequence starting with a negative element and ending with a negative element

For the state where $$$R - L = 0$$$,

$$$1 = A_l$$$
$$$2 = -10^{16}$$$
$$$3 = -10^{16}$$$
$$$4 = -A_l$$$

In the conquer stage, we need to merge these properties. We can easily do this as follows:

$$$1 = $$$ max(1L, 1R, 1L + 3R, 2L + 1R)
$$$2 = $$$ max(2L, 2R, 1L + 4R, 2L + 2R)
$$$3 = $$$ max(3L, 3R, 3L + 3R, 4L + 1R)
$$$4 = $$$ max(4L, 4R, 4L + 2R, 3L + 4R)

Our final answer will be $$$max(1, 2)$$$ over the entire array. This divide and conquer can solve C1. The question is, how do we keep updating this?

Well, in this 'divide and conquer', the conquer part is only $$$O(1)$$$. Therefore, this is a decomposable search problem which can be handled fast by segment trees.

Therefore to AC C2 we can create a segment tree where each node in the segment tree stores the 4 properties mentioned above and the merge function performs the max operations as described above. After this, it's standard.

Time Complexity: $$$O(N + Q \cdot LogN)$$$

You can view my submission for clarity:

https://mirror.codeforces.com/contest/1420/submission/93674507

Note: There was a much easier solution using some problem specific observations, but it's always better to have a solution that works on a wider range of cases. I hope this helped :)

To make it cleaner change 2 and 3 to "Maximum sum subsequence of even size starting with a positive element" and "Maximum sum subsequence of even size starting with a negative element", this removes isolated 1L 1R 2L 2R from the max, because the base case turns into 0 (empty subsequence has even size) instead of -10^16.

Yeah GG

kshitijg9899

I liked the fact that hints are provided before the full solution. Thank you so much!

utr491

What is event method mentioned in D's solution?

NewRules

We intend to keep track of the number of lamps which are switched ON at a particular instance of time,let's call that variable x. Now consider turning ON and turning OFF of each lamps as seperate events.It's intuitive that the numbers of lamps switched ON (variable x) only changes at the time of an event(switching ON or switching OFF of some lamp). Hence we sort all the events in non-decreasing order of their occurence and process them one by one. (Note that for some time T if there are multiple events occuring ,we place the switching ON events before the switching OFF events because the question mentions that lamps are switched ON from li to ri inclusive).

As we process the events 1by1, if we encounter a switching ON event, x is incremented, otherwise x is decremented.(quite intuitive). Whenever there is a switching on event and x>=k-1 , there are X choose K-1 ways of selecting k-1 lamps ,which cobmined with the current lamp being switched on gives us a new solution of k lamps.(My Submission)

If you are new to this technique, https://www.geeksforgeeks.org/minimum-number-platforms-required-railwaybus-station/ is a good introductory problem.

Newton_01

Thanks for your explanation. I had this thought first but hesitated in doing it.

Thank you for the explanation. For anyone who did not get it, an alternate explanation could be that we just have to chronologically sort all that happens and iterate over it while maintaining a counter of the number of bulbs still switched on. Just keep in mind to switch off a bulb after switching on all the bulbs at that instant.

lmhoang

+25

I found an AC submission to problem B in my room:

https://mirror.codeforces.com/contest/1420/submission/93659129

It has O (n^2) complexity but runs surprisingly fast due to the use of intel vectorization.

Is this possible on AMD processors?

AMD processors have the same instruction set as Intel and support vectorization, so I think it's possible.

VishalGaurav

Does this mean that we can just pass simple n2 loops like these with this vectorization technique?

Sometimes it helps to fit into time limits, but these cases are not so often. Vectorization improves speed only in constant time, so if I make $$$200'000$$$ instead of $$$100'000$$$, then the $$$O(n)$$$ solution will be two times slower, but vectorized $$$O(n^2)$$$ will become four times slower, because it's still $$$O(n^2)$$$. Additionally, not every solution is good to vectorize.

What is intel vectorization and why this O(n^2) solution passing??

In short, it allow the core perform multiple instructions simutanenously.

frauddd

OMG!! I tried my luck with the exact same approach but got tle.....Can you plz explain what is Intel vectorization?

You need a new generation of intel cpu e.g. 9th or 10th gen. Older cpus do not support this feature.

So it means who has high processor laptop dont need to worry about time complexity??

+10

No, having a better algorithm is always better.

Vectorization is available for long time. According to Wiki, AVX2 is supported in Intel CPUs since 2013, and SSE is used for even longer period. Why do you think that only 9th or 10th generations are supported?

I guess so because this code runs on my computer in 6 seconds (intel xeon 6th gen).

1) What compiler are you using?

2) Did you see the generated assembly?

If your CPU doesn't support AVX2 instructions and the compiled binary has them, the program will simply crash.

I am using GCC 8.2 Windows, C++14. The compiler bypass all directives and generate the output exe as the same as normal

The compiler bypass all directives and generate the output exe as the same as normal

Then your compiler didn't produce such instructions. It doesn't mean that CPU is not capable of executing AVX2.

If you try to compile https://mirror.codeforces.com/contest/1420/submission/93659129 with -O2, everything should be optimized and AVX2 will be used.

Yash2000

So does this mean that someone with 9th/10th gen processor can get AC easily by doing 10^10 operations in like 1 second?

It's not true for every program. Vectorization improves the speed much only in the specific parts, like short loops which iterate over the array. You cannot just write #pragma GCC target("avx,avx2") and make any program two times faster. But there are some problems on Codeforces where solutions in $$$O(n^2)$$$ with vectorization passed with $$$n \le 100'000$$$.

Lelouch.Lamperouge

Clean solution to C1 using maxima minima approach

CODE

#include <bits/stdc++.h>
#define int long long 
using namespace std;
void solve(){
  int n,q ;cin >> n >> q ;
  vector<int>a(n) ;
  for(int &x:a)
    cin >> x ;
  a.insert(a.begin(),(int)-1e9) ;
  a.push_back((int)-1e9) ;
  int ans = 0 ;
  for(int i=1;i<n+1;i++){
    if(a[i]>a[i-1]&&a[i]>a[i+1])
      ans+=a[i] ;
    if(a[i]<a[i-1]&&a[i]<a[i+1])
      ans-=a[i] ;
	}
  cout << ans << '\n' ;
}
signed main() {
  int t ;cin >> t ;
  while(t--)
    solve() ;
}

Anachor

+35

For C2, the answer can be written as $$$\sum_{i=0}^{n}{max(0, a_i-a_{i+1})}$$$, where, $$$a_0 = a_{n+1} = 0$$$. Viewed this way, updates become very easy.

SCZ.Official

Actually , problem E has a solution in only $$$O(N^4)$$$ with simple convex hull trick . You can check my solution here .

CH_wzy

Fantastic!

doublevgp

Can anyone explain to me the problem b's code? why j from 29 to 0?

ankitsingh123

It is because 1<= Ai <=10^9. So Ai can have maximum length of 29 in binary form.

ok thx!

lucifer1004

For 1420C2 - Pokémon Army (hard version), instead of storing the local maxima/minima, we actually only need to store the current value of each element.

Why?

Because the best answer we can get can also be represented in

$$$\sum\max(0, a_i-a_{i+1})$$$

If we add an extra $$$0$$$ to the end of the array.

So for each query, we need to check at most 4 pairs: $$$(l-1,l)$$$, $$$(l,l+1)$$$, $$$(r-1,r)$$$, $$$(r,r+1)$$$. Note that there might be duplicates, so we need to use a set to deduplicate.

If the original pair is counted in the current sum, then we subtract it from the sum.

Then we do the swap and check the pairs again. If the new pair should be counted, we just add it to the sum.

I think this solution is clearer than the official one.

AC submission: 93743846

zendaya007

i am getting wrong answer for test case 44 in problem D. please help me out. i have done the same thing as mentioned in editorial most probably. submission : 93742154

Please help lucifer1004

This solution is really nice :) Do you have an AC submission?

I have added my submission.

I think it should be $$$\sum max(0, a_{i+1} - a_{i})$$$, $$$i = 0$$$ to $$$n$$$, $$$a_0 = 0$$$ and $$$a_{n+1} = 0$$$

The two forms are equivalent.

Note that I only add one $$$0$$$ at the end, while you add $$$0$$$ at both front and end.

EN_SA

I really like this editorial providing hints not only solutions.

sarthak.ag

+13

Actually DP solution with $$$O(n + q\sqrt{n})$$$ is also conceivable for C2 although because of hard constraints it gets a TLE. My approach, drawing inspiration from Mo's algorithm, was basically to split the pokemons into $$$\sqrt{n}$$$ segments of size $$$\sqrt{n}$$$ each. For each segment we can calculate following 4 things in $$$O(\sqrt{n})$$$ time (using the approach in C1). 1. Max +- sequence with even number of terms. 2. Max -+ sequence with even number of terms. 3. Max +- sequence with odd number of terms 4. Max -+ sequence with odd number of terms.

Finally we can run the dp on these $$$\sqrt{n}$$$ segments to get the answer in $$$O(\sqrt{n})$$$ time. Now, for every query just swap the pokemons and then recompute the (atmost) 2 segments in which the swapped pokemon lie and finally recompute the DP. For each query this can be done in $$$O(\sqrt{n})$$$ time.

My solution, TLE ofcourse - 93717594

rahul_2

Can someone Please explain how to do c2 with segment trees?

https://mirror.codeforces.com/blog/entry/82978?#comment-701497

Read the comments above before posting yours.

prakhargupta.pg

Why is my code getting runtime error on test 5 of problem c1. I have applied DP using 2d array where dp[i][j] tells answer when i elements are considered with j elements included in subsequence. https://mirror.codeforces.com/contest/1420/submission/93757697

I will be obliged if someone helps.

It can happen that $$$n = 3\cdot10^5$$$ in this task. And $$$3\cdot10^5\times3\cdot10^5$$$ is too large to fit into memory on stack, so the program crashes.

PremKapshi1234

When will ratings update?

Whats the reason behind the delay?

Overstars

For C2, we can use set to avoid categorical discussions.

set<int>now{l-1,l,l+1,r-1,r,r+1};
for(auto &x:now)
	erase(x);

Can anyone tell me about the "event metho" which author is referring in solution of D? "It should be noted that p(x) and s(x) can be easily supported using the event method. Then, the total runtime will be O(nlogn)."

The idea is as follows. You have segments $$$[l_i; r_i]$$$ and want to calculate (for example) the number of segments that contain specified points $$$x_i$$$. To do this, you create an array of events:

a segment is started (at time $$$l_i$$$)
a point is encountered (at time $$$x_i$$$)
a segment is finished (at time $$$r_i$$$)

Then, sort such events by time and process them in sorted order. Now it's easy to keep track of number of open segments in the specified points $$$x_i$$$.

tejas_warambhe

I still don't understand the logic behind A Cube Sorting

It is not difficult to see that the answer «NO» in this task is possible when and only when all ai are different and sorted in descending order.

How can one derive such logic

PS — newbie here, want to learn :D

Last_whisper

If you know something about bubble sorting, you can think of the worst time complexity of bubble sorting ($$$\frac{n-1 \cdot n}{2}$$$) and its reasons

Giyan

In 1420B- Rock and lever can you explain what this part of the code is doing?

for (int j=29; j>=0; j--) { int64_t cnt=0; for (int i=0; i<n; i++) { if (a[i]>=(1<<j)&&a[i]<(1<<(j+1))) { cnt++; } } ans+=cnt*(cnt-1)/2; }

PS-Newbie Here

if (a[i]>=(1<<j)&&a[i]<(1<<(j+1))) { cnt++; }

stores the no of ( numbers having their rightmost set bit e.g: 4(binary representation : 001) ,in 4 the rightmost set bit is at index 2(0 indexed) ) .

cnt*(cnt-1)/2 part counts the nos satisfying the condition mention in the question for that count.

sumit_gaurav

-9

Hello everyone !! I hope you are having a good day!

i tried to make editorial for this round :- https://www.youtube.com/watch?v=uc6mN7BZVbI&list=PLrT256hfFv5WH3U33DuKUzL9fYnVsTGAj

problem solved. -> A , B , C1 , C2 , D

language of communication -> Hindi

programming language -> C++

seniLogik

in problem c2 let y be the total number of local maximum and minimum we have now, can we change the parity of the y with an update(swapping two element)?

Believe

Can someone explain the dp which we have done in ques E. Thanks

DDDDDDDDOOOOOOOONNNNNNNN

Can anyone help me with the problem D. I got WA in test 11. Here is my code:93779028

If I am not wrong you have compressed the values but you cannot compress the values in this question. You have to keep them as it is.

Can you explain why I cannot compress the values? I think the values are not important in this problem.

For example, two segments [1,4] and [2,6] after being compressed will become [1,3] and [2,4], which will not affect the result.

galen_colin

Here are video solutions for all problems, including a description of the weird segtree sol for C2

aditya_sheth

Very nice Editorial for problem E. Thanks :). Has anyone solved it using flows?

Ertugrul_

problem D is a nice educational problem thanks :)

jahnvi_13

Why is my code for B failing on 2nd pretest? I have stored the count of occurrences of each MSBs in a map and then have iterated over it to calculate the final count. Link to submission: https://mirror.codeforces.com/contest/1420/submission/93796185

39dll

Because when you are adding m(m-1)/2s, m can be up to 10^6, so if you set m as an int, it can overflow.

rkguptagkp462

I today constucted a easy way to do the C2 problem or C1 problem. The ans is essentially the half of sum of absolute value of adjacent indices. (To work it, we have to add zero at start and end of the array). Now the update is really easy, we have to just remove the old absolute values and add the new ones. This can also handle updates like changing the strength of a indice. See the Ac code for implementation details- AC code

h_squared

Can someone spot why I am getting WA on test 11 in Problem D. Here is my submission. 93837255

breddreamer

I like this step by step style of the editorial. Thanks.

laoji

In the tutorial of problem E, should $$$z_i$$$ be $$$\sum^i_{j=1}{f_j}$$$? And could you please explain why it takes $$$|z_i-(s+h)|$$$ steps from $$$dp_i$$$ to $$$dp_{i+1}$$$?

In the tutorial of problem E, should $$$z_i$$$ be $$$\sum\limits^i_{j=1} f_j$$$?

Yes, there is a typo in the editorial.

And could you please explain why it takes $$$|z_i-(s+h)|$$$ steps from $$$dp_i$$$ to $$$dp_{i+1}$$$?

Remember the reasoning above on how many steps it takes to transform the array $$$a$$$ into $$$b$$$. It's $$$\sum\limits_{i=1}^n \left|z_i - \sum\limits_{j=1}^n b_j\right|$$$, as shown above. In our DP, we know the array $$$a$$$ and build an optimal array $$$b$$$ step by step. After deciding that the value of the $$$i$$$-th element of the array $$$b$$$ will be $$$h$$$, we say that $$$\sum\limits_{j=1}^i b_j = s + h$$$. So, we add $$$\left|z_i - \sum\limits_{j=1}^i b_j\right| = \left|z_i - (s + h)\right|$$$ to the number of steps.

Is it now clear for you?

Yes, thanks. I didn't realize $$$dp_i$$$ starts from 0 so I mistook $$$h$$$ as $$$b_{i+1}$$$.

dantrag

I think there is some inconsistency in the editorial text. $$$s + h$$$ should be of the same "size" as $$$z_i$$$, so there I would rather write $$$z_{i+1}$$$.

ShubhamAvasthi

Can someone please tell me why I am not getting TLE on my solution to problem E?

It seems like magic. I really don't understand how a not very well optimized code like this is working.

Here is the solution link: 93898124.

Isn't it just $$$O(n^5)$$$? You have very straightforward loops, so it should be pretty fast. Switching to plain arrays might speed it up, but otherwise this looks like an intended solution.

Yes, but $$$80^5 > 32 \times 10^8$$$ and therefore, should not work, unless it is optimized heavily.

There is a small constant in front, rep(i,1,cnt+1) rep(j,i,n+1) is $$$n^2/2$$$, rep(p,0,j) further gives $$$n^3/3$$$. $$$10^8$$$ is quite okay, it is also very simple operations inside the loops.

deep_tricks

Can anyone help me with the reason for this solution getting TLE? Problem c2 https://mirror.codeforces.com/contest/1420/submission/93901635

Use '\n' instead of endl and use fast I/O.

leonidas328

-7

In editorial of 1420C1 - Pokémon Army (easy version) and 1420C2 - Pokémon Army (hard version).

The outputs of the editorial codes themselves do not match for the same input test cases. For example,
Input
1
4 0
5 2 2 5

Output (C1 Editorial Code)
8

Output (C2 Editorial Code)
10

I know it's a bit late, but the editorial remains as is. So, this should be corrected.

Is the input correct? The numbers in the input must be distinct, according to the statement.

Oh, I'm really sorry for that. Thanks for correcting me.

prmondal

Can anybody explain how this formula is deduced in solution of problem E?

$$$p(A) = \sum_{1 \le i < j \le k} f_i\cdot f_j = \frac 12 \left( \sum_{1 \le i, j \le k} f_i\cdot f_j - \sum_{i=1}^k f_i^2 \right) = \frac 12 \left( \left(\sum_{i=1}^k f_i \right)^2 - \sum_{i=1}^k f_i^2 \right)$$$

We can get this equation in an alternative way:

$$$\displaystyle\left(\sum\limits_{1 \le i \le k} f_i\right)^2 = \sum\limits_{1 \le i \le k,\ 1 \le j \le k} f_i f_j = \sum\limits_{1 \le i \le k,\ 1 \le j \le k,\ i < j} f_i f_j + \sum\limits_{1 \le i \le k,\ 1 \le j \le k,\ i > j} f_i f_j + \sum\limits_{1 \le i \le k,\ 1 \le j \le k,\ i = j} f_i f_j = $$$

$$$\displaystyle = 2\sum\limits_{1 \le i < j \le k} f_if_j + \sum\limits_{i=1}^k f_i^2 \implies 2\sum\limits_{1 \le i < j \le k} f_if_j = \left(\sum\limits_{1 \le i \le k} f_i\right)^2 - \sum\limits_{i=1}^k f_i^2$$$.

In simple words, we take $$$\left(\sum\limits_{1 \le i \le k} f_k\right)^2$$$, which is a sum of pairwise products. Split these pairwise products onto three parts: where $$$i < j$$$, where $$$i > j$$$ and where $$$i = j$$$. The first two parts are equal to $$$\sum\limits_{1 \le i < j \le k} f_if_j$$$, and the last part is a sum of squares. From this we get that $$$2\sum\limits_{1 \le i < j \le k} f_if_j = \left(\sum\limits_{1 \le i \le k} f_k\right)^2 - \sum\limits_{i=1}^k f_i^2$$$, therefore the original formula is also true.

I did not get this part: "$$$\left(\sum\limits_{1 \le i \le k} f_i\right)^2$$$ = sum of pairwise products". Can you explain how?

Okay, do you agree that

$$$(f_1 + f_2 + f_3)\cdot(f_1 + f_2 + f_3) = f_1(f_1 + f_2 + f_3) + f_2(f_1 + f_2 + f_3) + f_3(f_1 + f_2 + f_3) =$$$

$$$= f_1f_1 + f_1f_2 + f_1f_3 + f_2f_1 + f_2f_2 + f_2f_3 + f_3f_1 + f_3f_2 + f_3f_3$$$?

(So, we got that $$$(f_1 + f_2 + f_3)^2$$$ is a sum of pairwise products of $$$f_1$$$, $$$f_2$$$ and $$$f_3$$$.

You can do this for arbitrary number of elements in the same manner.

Yeah got it. Thanks.

EDIEEee

I have an alternative DP solution to problem E. Observe that the answer is at most C(# zero , 2), in stead of counting the good pairs, I count the bad pairs. (The number of bad pairs is the sum of C(len[i] , 2) for each consecutive 0s with length len[i])

Firstly, I considered interval DP, however, we need also keep track with the number of orders, so we have a state at least DP[l <= 80][r <= 80][k <= 80 * 79 / 2], when merging interval, the time complexity is at least O(N^6), so this doesn't seem like a good approach. Then, I consider a left to right solution, let DP[i][j][k][l] = the minimal bad pairs for the first making length of i sequence such that, we have used j orders, and have a length k full of consecutive 0s suffix, and have used l ones. For transition, DP[i][j][k][l] -> DP[i + 1][j][k + 1][l] + k if we put 0 on the (i + 1) th position. DP[i][j][k][l] -> DP[i + 1][abs(pos[l + 1] — (i + 1))][0][l + 1]. where pos[i] is ith one's position in the initial sequence.

Time complexity O(N^5 / 8) Runtime is some thing like O(80 * 80 * 80 * 79 / 2 * 40 * 40) which should be ok for 2 seconds.

My code: https://mirror.codeforces.com/contest/1420/submission/95068002

eidan

I read this blog while upsolving problem E. I want to thank gepardo for such a clear, detailed and well-written editorial. The hints were accurate, the math formulas nice and understandable, the explanations concise, I can even say I enjoyed reading it.

Recently, I've stumbled upon editorials that aren't very good. They're written in a hurry, and sometimes not even clear or complete. This editorial was the exception, and for that, the author deserves recognition. It's a great example of how problem tutorials should be written.

Thanks :)

MuhammadSawalhy

15 months ago, # ^ |

So sice editorial... ❤️

tgp07

3 years ago, # |

My code for problem D is giving me a run time error.

153332275.

Could anyone please tell me what I'm doing wrong here?

2 years ago, # |

Can anyone please explain Solution B.Not able to understand it

vitosevskich

See my submission: 169305548.

Let's calculate the first bit that is 1 in every $$$a[i]$$$ from $$$1$$$ to $$$n-1$$$ inclusively going from highest to lowest bit. (function f) Let's call it $$$f(a[i])$$$.

Observation is: $$$a[i]$$$ & $$$a[j]$$$ $$$>=$$$ $$$a[i]$$$ ^ $$$a[j]$$$ if $$$f(a[i])$$$ = $$$f(a[j])$$$.

Why? If $$$f(a[i])$$$-th bit in $$$a[j]$$$ is 1 then value of that bit when calulating $$$a[i]$$$ & $$$a[j]$$$ will be 1 and when calculating $$$a[i]$$$ ^ $$$a[j]$$$ will be 0. If $$$f(a[i])$$$-th bit in $$$a[j]$$$ is 0 then value of that bit when calulating $$$a[i]$$$ & $$$a[j]$$$ will be 0 and when calculating $$$a[i]$$$ ^ $$$a[j]$$$ will be 1.

It's known that when comparing two numbers, first that have 1 bit is bigger. Because $$$2^n$$$ $$$>$$$ sum of all powers of 2 from $$$0$$$ to $$$n-1$$$ inclusively. (I know there is way to write it using LaTex but I don't know how to do it.)

Hope it helped.

gepardo's blog