You are not logged in. Please login at www.codechef.com to post your questions!

×

# SUGGESTION:

Before posting your question with asking for help try this test.
The answer should be 436746489.
If it is hard to debug this test for you, here is the helpful information.
It contains the number of unsuccessful tries for each cell
or -1 if the cell was hacked by the "Grid Hacking Mechanism".
If it still hard to debug try another smaller test that contains answer and similar debug information as above.
If you pass this test then read carefully the bold tip in the QUICK EXPLANATION section.
If you follow this tip but still have WA then post your question and we will help you.

# PROBLEM LINK:

Author: Khadar Basha
Tester: Anton Lunyov
Editorialist: Anton Lunyov

EASY

# PROBLEM:

You are given N × N grid filled by non-negative integers. We single out 3 types of numbers on the grid: primes, even non-primes and odd non-primes. For each number we define the cost as follows: for the prime number it is the id of this prime in the 0-based list of all primes, for even non-prime number X it is X / 2, and for odd non-prime number Y it is (Y + 3) / 2 (the cost of a number is the shortcut for the number of unsuccessful tries to crack the server secured by the password equals to this number; the mentioned formulas for cost in all three cases are clear from the problem statement).

Two even non-prime numbers that share a side on the grid are connected to each other. The same is true for odd non-prime numbers. But prime number is not connected to any other number. In this way all cells of the grid form a bidirectional graph that has several connected components (each cell having prime number is a separate component). The cost of the component is defined as the cost of the first cell in this component that we meet traversing the grid in row-major order. Now your task is to find the sum of costs over all components.

# QUICK EXPLANATION:

Actually, most of the job has already made while we reformulated the problem above. Now you simply need to loop over all cells of the grid in row-major order and if we meet the unvisited cell we add its cost to the answer and run DFS from this cell in the graph constructed above to mark other cells of its connected component as visited.

Tip: the total cost could overflow 32-bit integer type. Use 64-bit type instead.

To be able to find the cost of the prime cell fast enough we need to run Sieve of Eratosthenes in advance (even before input the very first test) in order to find all prime numbers. Then we need to create the array of size 107 that contains for each prime its id. Now the cost of each number can be found in O(1).

The overall complexity is O(K * log log K + T * N * N).

# EXPLANATION:

As mentioned above at first we run Sieve of Eratosthenes to identify all prime numbers:

isPrime[0] = false
isPrime[1] = false
for i = 2 to K do
isPrime[i] = true
for i = 2 to sqrt(K) do
if isPrime[i] then
for j = i * i to K with step i do
isPrime[j] = false

where K = 107 − 1.

Next we fill array of costs. We maintain variable prime_id which is equal to the number of primes we met so far:

prime_id = 0
for i = 0 to K do
if isPrime[i] then
cost[i] = prime_id
prime_id = prime_id + 1
else
cost[i] = i div 2 + (i mod 2) * 2

Note the formula for the cost for non-prime number.

Now we can input grids and traverse them. We use two-dimensional array A[][] for the initial grid. Also we use another two-dimensional array visited[][] to check visited cells, which could be filled by false while input the grid:

for i = 1 to N do
for j = 1 to N do
read A[i][j]
visited[i][j] = false


Now we can traverse the grid. If we meet visited cell we skip it. Otherwise we always add its cost to the answer and if it is non-prime then run DFS to fill its connected component:

total_cost = 0 // this should be a variable of 64-bit integer type!
for i = 1 to N do
for j = 1 to N do
if (visited[i][j]) then
continue
total_cost = total_cost + cost[A[i][j]]
if (not isPrime[A[i][j]]) then
DFS(i, j)


DFS(i, j) is the recursive routine that traverse the grid passing from the cell (i, j) to its neighbors in the graph constructed above:

DSF(i, j)
if (not visited[i][j]) then
visited[i][j] = true
for (x, y) in {(i+1, j), (i-1, j), (i, j-1), (i, j+1)} do
if (not isPrime[A[x][y]]) and (A[x][y] is of the same parity as A[i][j]) then
DFS(x, y)

Note that we pass to the cell (x, y) only if the number in it is non-prime and of the same parity as in the cell (i, j). Missing of any of these checks will definitely lead to WA. Also see tester's solution as a reference to one of the convenient ways how to implement loop:
        for (x, y) in {(i+1, j), (i-1, j), (i, j-1), (i, j+1)} do


# ALTERNATIVE SOLUTION:

For some languages (like Java or Python) the recursive implementation of DFS could lead to runtime error due to stack overflow. In this case you need to implement either non-recursive DFS or BFS (breadth-first search) to mark cells of each connected component. Another alternative is to use disjoint-set data structure to do the same job. But this way requires some additional thinking.

The alternative very fast method to find all primes is to use Atkin Sieve invented in 2004 by Atkin and Bernstein. It allows to find all primes up to K using O(K / log log K) operations. Check it out first 6 related problems to train yourself at fast implementation of sieve.

The remaining 3 problems involves flood fill algorithm.

# AUTHOR'S AND TESTER'S SOLUTIONS:

Author's solution can be found here.
Tester's solution can be found here.

# RELATED PROBLEMS:

This question is marked "community wiki".

asked 15 Jan '13, 15:00

6.7k62119138
accept rate: 12%

fantastic editorial. loved it!

(15 Jan '13, 16:25)
3

@editorialist : great work to give additional links to similar problems :) thanks. Things are really changing nicely in new year :)

(15 Jan '13, 20:20)

Its Nice.....

(16 Jan '13, 15:04)

25 Answers:
 7 only mistake was... i used 32bit integer for total cost instead of 64 feels like kill-me-now answered 15 Jan '13, 16:42 2★charchit 126●1●4●9 accept rate: 0% I'm in the same page :(! (15 Jan '13, 17:11) tyrant2★ same here.. :( (15 Jan '13, 18:36) kshitij83★
 1 Hi, Actually i did the same solution described above in Python but it didn't work (i got Time Limit Exceeded), you can check my attempt here and here. But actually, i was able to find another algorithm for this problem that got accepted, the idea was that i only check for immediate neighbors of any server to know if it should be included or not, of course this may lead to count a server while it shouldn't (i.e. it should have been already cracked), but to be aware of this case i created a Python dictionary that hold information about which server cracked which one (cracked_by) and then each time i visit a server i check if i should update the cracked_by of the immediate neighbors to know if i added wrongly a server if yes i discard the value from the result. You can check my solution here, feel free to post any feedback. HTH, answered 16 Jan '13, 02:52 2★mouad 31●2 accept rate: 0% Hm... Quite neat. I was stuck for a while thinking that you process this test 3 4 1 4 4 4 4 1 1 1 incorrectly. But are you sure that this solution has complexity O(N * N) per test? What about test of the structure: 10 4 1 4 1 4 1 4 1 4 1 1 1 1 1 1 1 1 1 1 1 4 1 4 1 4 1 4 1 4 1 1 1 1 1 1 1 1 1 1 1 4 1 4 1 4 1 4 1 4 1 1 1 1 1 1 1 1 1 1 1 4 1 4 1 4 1 4 1 4 1 1 1 1 1 1 1 1 1 1 1 4 1 4 1 4 1 4 1 4 1 1 1 1 1 1 1 1 1 1 1 It seems to me that finding roots of paths here will be O(N) in average making your solution to be O(N3). (17 Jan '13, 11:54) If I'm correct, it would only mean that our test data are weak :) Well, it's impossible to predict with what approaches contestants can come up with to have a test failing each approach. (17 Jan '13, 12:02) @anton_lunyov: First of all thanks for your feedback, as for the two test cases that you gave, AFAIK my solution will return 6 and 52 respectively which is correct, and as for the complexity well i am not sure what is the correct answer ;-) specially i don't know what is the worst case to get roots of the two paths, but what i can tell you is that it will only take O(n**2) for the last test case that you gave, finding the roots will be only O(1) (2 extra iteration at max). (19 Jan '13, 02:44) mouad2★
 0 Can you that what is my problem ID is :1681905 answered 15 Jan '13, 15:59 43●6●10●14 accept rate: 0% 2 Did you test you program at least somehow? They produce wrong result on any possible test. I don't know. Consider for example this one: 1 1 3 The correct answer is 1, while you return -0 (yeah, with this crappy minus sign) One suggestion is that you should use long long type for unsec and print it as printf("%lld ",unsec). (15 Jan '13, 16:15) are you sure this test case my program give wrong answer because in my computer answer is 1 not -0.are you sure that you check correct submit my program id is:1681905 (15 Jan '13, 17:00) 1 I look here http://campus.codechef.com/files/stderr/1681905/ Only admins, setter and tester have access to this. Your answer for 8th test which is very small test is DATASET NUMBER: 8 -0 -268156158598... -0 -427197407184... -125542034707... -0 -0 -0 where ... is some other digits (like hundreds of digits more) Try to run your code on ideone.com using GNU C++ compiler against this test case. (15 Jan '13, 17:22)
 0 The editorial is pretty nice. I wonder what was wrong with my implementation. Can you please have a look and let me know my error. It would be great. I'm unable to figure out. http://www.codechef.com/viewsolution/1718060 answered 15 Jan '13, 16:33 2★coolmind 16 accept rate: 0% 1 Try to use usual array instead of queue. You could simply have a loop for(int k=0;k
 0 Nice editorial. I have an unrelated question though. I used BFS instead of DFS. Everything else was same as you suggested in your editorial. But I was getting Wrong answer somehow. I tried it with many test cases generated by hand and it was working fine. How should I approach in finding problems in my code in such scenarios. Any help will be highly appreciated. BTW, my problem submission was http://www.codechef.com/viewsolution/1720890 If anyone can point to what test cases my program was failing, it will be really helpful. Thanks in advance. answered 15 Jan '13, 16:34 46●2 accept rate: 0% I used BFS as well. I got TLE. :( (15 Jan '13, 18:00) coolmind2★ @vikram535 You have repeated the same mistake as hundreds of contestants and tens of contestants who asking for help after the contest. You should use 64-bit integer type long long for the answer and output it either using cout or %lld specifier. Regarding the bfs, see my answer above in the response http://discuss.codechef.com/answer_link/5181/ In short, you should use usual array instead of queue. (17 Jan '13, 02:52)
 0 I first used an integer array of size 10e7 for the sieve, but it ran out of memory though...had to half the size of the array and cancel all the even numbers. answered 15 Jan '13, 18:28 2★jbupit 6●1●1●1 accept rate: 0% Actually the main your issue is declaring all arrays inside main(). They are stored in a stack which is quite limited and hence you received runtime error. The displayed memory 1500M is just the bug in codechef system. The stack size seems to be 32M hence initially you have used 40M and get RE but then half the size to 20M and all became fine. (17 Jan '13, 11:28) @anton_lunyov you mean to say if I declare my array globally I can use upto 1500M memory? (17 Jan '13, 19:33) charchit2★ Yes. Even 1536M :) (17 Jan '13, 19:36)
 0 What is problem in my solution.. I am getting WA http://www.codechef.com/viewsolution/1723138 answered 15 Jan '13, 18:42 -1●1●1●2 accept rate: 0% Your output for this test 4 0 4 8 12 0 1 9 24 0 1 9 12 0 8 16 18 is 16 while should be 2. (17 Jan '13, 11:29) How is it coming out to be 2. Can you please explain? (21 Jan '13, 13:13) Its just for testing '1' at place (1,1)... Got my mistake... its a silly one in function kr1(int,int)....wrote pass[i][j] instead of typ[i]][j]... damn shit.... (21 Jan '13, 18:18)
 0 To implement grid hacking in java, if i use an ArrayDeque or a LinkedList, i get TLE( http://www.codechef.com/viewsolution/1724620 ), while if i change the data structure to ArrayList( http://www.codechef.com/viewsolution/1724632 ), i get one of the fastest execution times in java. Any reason as to why this is happening? answered 15 Jan '13, 21:09 5★karan173 685●9●11●19 accept rate: 0%
 0 Thanks for such a detailed editorial. I learnt a lot. But still my program is not getting submitted. If I am not asking too much, can I get a sample test case for which my code is failing. My submission id is 1724659. Thanks. answered 15 Jan '13, 21:35 1 accept rate: 0% 1 @sarfarajey Read the bold tip in the beginning of the editorial. (17 Jan '13, 01:48)
 0 I had the same solution in Java (except I used the Sieve of Atkin and instead of a boolean[][] I changed the original input array directly), but I kept getting WA :/. Can somebody check my solution to see what's wrong? I tried several test cases, but all of them printed out the right answer... http://www.codechef.com/viewsolution/1715077 answered 16 Jan '13, 01:50 2★kullalok 1.5k●11●22●36 accept rate: 14% You have repeated the same mistake as hundreds of contestants and tens of contestants who asking for help after the contest. You should use 64-bit integer type for the answer. Also your solution has stack overflow error: Exception in thread "main" java.lang.StackOverflowError (17 Jan '13, 11:33)
 0 I got AC after using Sieve of Eratosthenes but I want to implement it by another method by which I couldn't. We can find all prime numbers less than 10^7 by simply making another prime number program and copy the output and assign it to array in our program arr[]={all prime numbers less than 10^7 that we already have.}. Now just apply a for loop for(i=0;i
 0 Can the admins tell what's wrong with http://www.codechef.com/viewsolution/1703688? In other note :- The tester solution is not getting complied on gcc. Need to include stdio.h. answered 21 Jan '13, 12:49 31●1●5 accept rate: 0% 1 Your code is very long. Delete at least all this annoying debugging printf before asking for help. (21 Jan '13, 15:04) New id is 1742822. Thanks. (21 Jan '13, 16:05) Try this test: http://pastebin.com/raw.php?i=j6UWCAZ5 The correct answer is 436746489. While your answer is 458527862. (21 Jan '13, 23:06) I have updated the editorial. Please read and follow the SUGGESTION at the beginning of editorial. (22 Jan '13, 02:38) That's great anton_lunyvon. Thanks a lot. (22 Jan '13, 09:59)
 0 I'm getting WA my solution is this http://www.codechef.com/viewsolution/1742829. answered 21 Jan '13, 16:15 1●1●1●1 accept rate: 0% Hi There, your solution did not work on the following test case. 1 6 1034 5036 5076 8099 6543 7853 0 2 4 6 8 10 1 5 7 9 13 15 1034 5036 5076 8099 6543 7853 1034 5036 5076 8099 6543 7853 1034 5036 5076 8099 6543 7853 The answer is 9076. It is giving 5087. (21 Jan '13, 16:27)
 0 i have checked all the test cases which u have provided in the editorial and it gives right answer.I am unable to figure out the problem. http://www.codechef.com/viewsolution/1744261 answered 22 Jan '13, 09:59 1 accept rate: 0% In routines oddserver and evenserver you should check before accessing valid[k][l] (or valid[m][n]) that cell (k,l) (or (m,n)) lies inside the grid. Probably this is the reason. (22 Jan '13, 12:36) @anton_lunyov (valid[i][j] is itself a array which has 1 at (i,j) if i,j lies inside the grid otherwise valid[i][j] contains 0) it is also workin well .please help to sort out this problem ,still unable to fiure out the problem (23 Jan '13, 15:53) In evenserver() instead of  if(!valid[k][l]) return 0;  you should write  if(k<0 || k>=n || l<0 || l>=n || !valid[k][l]) return 0;  where n is the size of the grid that should be defined globally. The same should be applied to oddserver() but be careful you are using variable n there as y-coordinate of the cell. If you will have troubles with fixing it you can refer to my edit of your code that got AC: http://www.codechef.com/viewsolution/1745787 (23 Jan '13, 16:04)
 0 I habe a problem in my solution http://www.codechef.com/viewsolution/1745646 based on the editorial which seems to be more a compiler/c++ problem then an algorithmic problem. I got the error "Access violation reading location 0x01A4B7A4" when trying to read from my big prime array. Can anyone help? answered 23 Jan '13, 13:30 1 accept rate: 0% 1 Your array for primes is too small. Add one zero to K. Like here: http://www.codechef.com/viewsolution/1745675 (23 Jan '13, 14:15) Greatly appreciated, Anton! I should have noticed it by myself! (23 Jan '13, 22:48)
 0 Please have a look at the following solution : http://www.codechef.com/viewsolution/1746818 It gave correct results for the test cases mentioned in the editorial here, those in the problem statement and some cases generated by myself. But still it displays 'Wrong Result' here. Please help! answered 24 Jan '13, 13:19 2★avinrox 207●5 accept rate: 10% Use printf("%lld\n", tries). I should add this tip to the editorial. Several guys who asked for help all have this bug: they use long long but print it as %d (24 Jan '13, 14:27) thank you for the correction. i tried another approach.. facing a similar problem.. http://www.codechef.com/viewsolution/1747532 (24 Jan '13, 23:35) avinrox2★ Wow! Your solution fails for only one grid among 80 grids we have in official test data. It is a grid having the following structure: N=350, a[i][j] are all even non-primes, except very small percent of odd non-primes. Each group of odd passwords is generated as a boundary of a small square. Your solution is quite hard to follow so I can only suggest you to generate test of the above structure for smaller size and compare your answer with the answer generated by the correct program. (25 Jan '13, 01:08) approach is as : 1. generate prime numbers index. 2. read data (simultaneously add attempts for primes). store odd/even flag for remaining. 3. get_neighbors() will store 'valid' neighbors for each data. 'valid' means - (i) it is either [i+1][j] or [i][j+1]. (ii) it has to be of same type (odd/even) non-primes. 4. 'tags[i][j]' for each member is assigned as 'i*N+j'. 5. arrange_and_compute() will do traversal over all members and update 'valid' neighbors with lowest tags among them. this is done in loop until there is no change in an iteration. this was an attempt on amortization. will see! (25 Jan '13, 15:11) avinrox2★
 0 Can you please tell whats wrong with my code...its giving correct output for test cases.. http://www.codechef.com/viewsolution/1752815 answered 27 Jan '13, 15:54 1 accept rate: 0% You should not pass through the cells containing 2 in mark_check_even. Try this test: 4 0 4 2 12 2 1 9 24 0 1 9 12 2 8 2 18 ` The answer is 12 while your answer is 2. (27 Jan '13, 16:02) Thanks a lot... (27 Jan '13, 16:30)
 0 my code is giving a runtime error(NZEC).Can anyone plz tell what is the problem? answered 05 Feb '13, 10:55 1●1●2 accept rate: 0%
 0 Can anybody look at my solution and explain why am i getting Runtime error(segmentation fault). i cant see why..:( http://www.codechef.com/viewsolution/4125365 answered 22 Jun '14, 19:28 4★dcod 358●1●5●11 accept rate: 27%
 0 my code runs fine on my machine but shows wrong answer on submission Please help my solution is here http://www.codechef.com/viewsolution/5203370 answered 23 Oct '14, 23:06 1 accept rate: 0%
 0 @anton_lunyov @i_am_what_i_am Hi, I just studied DFS and then searched this problem through question tags. I passed the given test cases and also the cases given in editorial. But still couldn't get AC. I made these solutions: 1) https://www.codechef.com/viewsolution/8801554 ---> It gave me WA 2) https://www.codechef.com/viewsolution/8801563 ---> It gave me RE The only difference between both these solutions is that i have changed my counter variable which is needed to be printed from long to long long. I would be so glad if anyone could explain me the bug. answered 17 Nov '15, 18:26 6★torque 427●1●1●11 accept rate: 17%
 0 struggling from a day with my code here it gives correct answer to all cases ever in discussion or editorial but still getting a WA. help would be greatly appreciated answered 23 May '16, 19:02 1★sd_7 1 accept rate: 0%
 0 Hey I solved the problem with same approach as given above but it's showing wrong answer. Please can anyone help https://www.codechef.com/viewsolution/14628909 answered 21 Jul '17, 08:16 1 accept rate: 0%
 0 Hey I solved the problem with same approach as given above but it's showing wrong answer. Please can anyone help https://www.codechef.com/viewsolution/14628909 answered 21 Jul '17, 08:17 1 accept rate: 0%
 -1 The test data, http://pastebin.com/raw.php?i=j6UWCAZ5 is too large for me to debug the program, so I am unable to find why is it giving 421211088 instead of correct result. Can you please help. My solution is http://www.codechef.com/viewsolution/1696242 answered 21 Jan '13, 23:38 1★raman92 214●4●8●15 accept rate: 0% I have updated the editorial. Please read and follow the SUGGESTION at the beginning of editorial. (22 Jan '13, 02:41)
 toggle preview community wiki:
Preview

### Follow this question

By Email:

Once you sign in you will be able to subscribe for any updates here

By RSS:

Answers

Answers and Comments

Markdown Basics

• *italic* or _italic_
• **bold** or __bold__
• link:[text](http://url.com/ "title")
• image?![alt text](/path/img.jpg "title")
• numbered list: 1. Foo 2. Bar
• to add a line break simply add two spaces to where you would like the new line to be.
• basic HTML tags are also supported
• mathemetical formulas in Latex between \$ symbol

Question tags:

×15,852
×3,820
×733
×303
×22
×19

question asked: 15 Jan '13, 15:00

question was seen: 9,944 times

last updated: 21 Jul '17, 08:17