Digit DP - Codeforces

#	User	Rating
1	tourist	3985
2	jiangly	3814
3	jqdai0815	3682
4	Benq	3529
5	orzdevinwang	3526
6	ksun48	3517
7	Radewoosh	3410
8	hos.lyric	3399
9	ecnerwala	3392
9	Um_nik	3392

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

Wrote this article a long ago but during solving a problem recently thought of sharing this article publicly. Hope it will help some contestants to understand the idea clearly.

Digit dp is a very easy technique and also useful to solve many dynamic programming problems. Seeing the name “Digit DP” it’s easy to guess that we are going to do something using the digits. Yes we are actually going to play with digits. Let’s explain the concept using a classical problem.

Problem

How many numbers x are there in the range a to b, where the digit d occurs exactly k times in x? There may have several solutions including number theory or combinatorics, but let’s see how we can solve this problem using digit dp.

Solve for range (zero to a)

Using digit dp we always focus on building a number satisfying all the conditions. If we finally manage to build that number then we say, yes we have got one ;-) But how we’ll build that number? For the time being let’s say a is zero. So we need to find the total numbers which are not greater than b and also satisfy the given conditions.

Building a sequence of digits

Let’s consider the number as a sequence of digits. Let’s name the sequence sq. Initially sq is empty. We’ll try to add new digits from left to right to build the sequence. In each recursive call we’ll place a digit in our current position and will call recursively to add a digit in the next position. But can we place any of the digits from 0 to 9 in our current position? Of course not, because we need to make sure that the number is not getting larger than b.

Information we need to place a digit at the current position

Let’s say during the building of the sequence, currently we are at position pos. We have already placed some digits in position from 1 to pos-1. So now we are trying to place a digit at current position pos. If we knew the whole sequence we have build so far till position pos-1 then we could easily find out which digits we can place now. But how?

You can see that, in the sequence sq the left most digit is actually the most significant digit. And the significance get decreased from left to right. So if there exist any position t (1<=t<pos) where sq[t] < b[t] then we can place any digit in our current position. Because the sequence has already become smaller than b no matter which digit we place in the later positions. Note, b[t] means the digit at position t at number b.

But if there was no t that satisfy that condition then at position pos, we can’t place any digit greater than b[pos]. Because then the number will become larger than b.

Do we really need the whole sequence?

Now imagine, do we really need that whole sequence to find if a valid t exist? If we placed any digit in our previous position which was smaller than its corresponding digit in b then couldn’t we just pass the information somehow so that we can use it later? Yes, using an extra parameter f1(true/false) in our function we can handle that. Whenever we place a digit at position t which is smaller than b[t] we can make f1 = 1 for the next recursive call. So whenever we are at any position later, we don’t actually need the whole sequence. Using the value of f1 we can know if the sequence have already become smaller than b.

Extra condition

So far we focused on building the sequence sq, but we have forgotten that there is an extra condition which is, digit d will have to occur exactly k times in sequence sq. We need another parameter cnt. cnt is basically the number of times we have placed digit d so far in our sequence sq. Whenever we place digit d in our sequence sq we just increment cnt in our next recursive call.

In the base case when we have built the whole sequence we just need to check if cnt is equal to k. If it is then we return 1, otherwise we return 0.

Final DP States

If we have understood everything so far then it's easy to see that we need total three states for DP memoization. At which position we are, if the number has already become smaller than b and the frequency of digit d till now.

Solve for range (a to b)

Using the above approach we can find the total valid numbers in the range 0 to b. But in the original problem the range was actually a to b. How to handle that? Well, first we can find the result for range 0 to b and then just remove the result for range 0 to a-1. Then what we are left off is actually the result from range a to b.

How to solve for range a to b in a single recursion?

In the above approach we used an extra parameter f1 which helped us to make sure the sequence is not getting larger than b. Can’t we do the similar thing so that the sequence does not become smaller than a? Yes of course. For that, we need to maintain an extra parameter f2 which will say if there exist a position t such that sq[t] > a[t]. Depending on the value of f2 we can select the digits in our current position so that the sequence does not become smaller than a. Note: We also have to maintain the condition for f1 parallely so that the sequence remains valid.

Please check this to find the sample code of our initial approach.

Problem List

Is there some other problems? Some suggestions are really welcomed :)

Comments (105)

Show archived | Write comment?

minimario

7 years ago, # |

← Rev. 2 →

I've had this problem for a while, and most probably it's solvable with Digit DP. Maybe it's not, but I haven't found a solution:

Given a integer X, find the number of integers i in [l, r] such that i ≤ X and rev(i) ≤ X, where rev(i) is the number formed by reversing the digits of i. For example, rev(1560) = 651 and rev(156) = 651 (not 6510, 65100, etc)

Unfortunately I can remember neither the source nor the limits, but probably X < 10¹⁸ or something.

→ Reply

flash_7

7 years ago, # ^ |

I have found a solution but not fully sure about it.

First of all if X < l then the answer is zero. And i must be <= L where L = min(r, X). Now let’s solve the problem for range (0, L). Let's Len = number of digits in L. Now we can divide the solution in two parts.

In the first part, for each k (1 <= k < Len) we'll try to find the number of valid i which has k digit in it's representation. We can do this easily because if i <= X then rev(i) <= will also hold. As total number of digits is smaller. So we can just ignore that rev(i) <= X condition for this part.

Let's sum1 is the summation of result for each k (1 <= k < Len). Now each of the number x in sum1 will contribute to the final result twice. Because for each x we can add some trailing zeroes in x to make it a Len digit number. As we are adding some trailing zeroes the conditions i <= X and rev(i) <= X will still hold.

Now we’ll solve the second part where we need to find the number of Len digit number which has no leading or trailing zeroes. I'll explain this with an example. Let's L = 372967524. Let's divide this L into 3 parts. P1 = 372, P2 = 9, P3 = 67524. Here P2 has only one digit.

Let's build a number i whose first 3 digit is fixed (same as P1). And at the 4th position we place a smaller digit than P2 (for example 8). So i now looks like this 3 7 2 8 _ _ _ _ _. So i has already become less than X. It doesn't matter what we place in the remaining 5 positions. Condition i <= X will always satisfy. So now we focus on selecting 5 digits for the remaining position such that rev(i) <= X satisfy.

The reverse of the trailing 5 digits of i is actually the 5 leading digits of rev(i). As we don't need to worry about i <= L condition any more, so now we'll do digit DP to find this first 5 leading digits of rev(i), comparing to the first 5 digits in X. We need to handle one more case here. As we have already placed the leading 4 digits in i whose reverse is actually the trailing 4 digits of rev(i), so after choosing the leading 5 digits in rev(i) we just need to check if the whole rev(i) satisfy the condition (<= X) or not.

As i said above, we’ll have to divide L in 3 parts. For each digit in L we’ll select it as P2. And the remaining two parts as P1 and P2. Then we repeat the same approach each time. The summation of all these results(including the one in part1) is our final answer.

VladaMG98

Check BOI(Baltic) 2013 Numbers, i really liked that one.

sam29

+11

I think you can add this problem. I solved it using Digit DP. I really liked this problem.

UPD: Another interesting problem on DIGIT DP https://www.hackerrank.com/contests/morgan-stanley-codeathon-2017/challenges/dreamplay-and-clubbing/problem

Nice problem. Thank you :)

DebSourav33

6 years ago, # ^ |

Beautiful problem! I used digit dp, bitmask, binary search- all in 1 problem!!

Auto comment: topic has been updated by flash_7 (previous revision, new revision, compare).

BishalG

One may find this interesting.

Solved it. A very nice problem. Thanks Vai :)

xxxcxxx

How to solve LIDS ?

The maximum length of LIDS can be at max 10. So we can check for each length k, in how many ways we can make a number whose LIDS is k. We can then print the result we found for the maximum k.

I’m describing the states we need to find the number of ways we can build an LIDS of length k.

At which position we are.
Has the sequence already become less than y? (0 or 1)
Has the sequence already become greater than x? (0 or 1)
Did we place at least one non zero digit so far? (To handle the first digit of LIDS)
Last digit we have placed in the sequence.
Number of total digit in the LIDS sequence.

Basically we build two sequence parallelly. One is definitely the whole sequence, another one is the LIDS sequence (Which is the sub sequence of the original sequence). When ever we select a digit we want to place, it becomes the last digit of the original sequence till now. But we can also choose if we want to consider that digit in our LIDS sequence or not.

sm247

Why do need the fourth state ??@flash_7

To make sure we don't place any leading zeroes. Because if we build a number of length k with some leading zeroes then that number is not really a number of length k.

TopGun__

3 years ago, # ^ |

Yes, I debugged my code and was getting a WA for

10000 10000

my answer was 1,1 because it did not consider zeros in LIDS.

Actual answer is 1,5.

HynDuf

5 years ago, # ^ |

Could you provide more details about how the transitions of this DP are ? Or the accepted code so I could try to understand ?

Thanks !

Dardy

4 years ago, # ^ |

You probably don't need this anymore, but here is a solution: https://ideone.com/QGXuo9

littlecoder_

2 years ago, # ^ |

but how to count it can you please elaborate?

Megumin2006

For those who are struggling, here is my accepted implementation: https://ideone.com/gdavOn

Feel free to ask me anything. Glad if i can help.

hydra_cody

21 month(s) ago, # ^ |

Any help getting TLE on LIDS? https://pastebin.com/zLnPJiaW

jonnnnnsinaga

15 months ago, # ^ |

wouldn't it TLE ? Since the number of operation that we should do for all query is 10000*10*2*2*2*10 = 8000000

pseudonymous

Here is the solution to LIDS problem using 3 flags for people who are looking for it.

Code

const int N=11;
ll dp[N][N][N][2][2][2];
string a,b;
/*counts no. of subsequence of len==need using 
 remaining suffix starting with lastDigit set at (pos-1) + 3 flags*/
ll count(int pos,int lastDigit,int need,int flagA,int flagB,int flagZero)
{
        if(pos==b.size())
                return need==0;
        if(dp[pos][lastDigit+1][need][flagA][flagB][flagZero]!=-1)
                return dp[pos][lastDigit+1][need][flagA][flagB][flagZero];
        int l,r;
        if(flagA)
                l=0;
        else l=a[pos]-'0';
        if(flagB)
                r=9;
        else r=b[pos]-'0';
        ll res=0;
        for(int i=l;i<=r;i++)
        {
                res+=count(pos+1,lastDigit,need,flagA||i>l,flagB||i<r,flagZero||i);
                if(i>lastDigit && need && (flagZero || i))
                        res+=count(pos+1,i,need-1,flagA||i>l,flagB||i<r,flagZero||i);
        }
        return dp[pos][lastDigit+1][need][flagA][flagB][flagZero]=res;
}
void solve()
{
        int x,y;
        cin>>x>>y;
        a.clear(),b.clear();
        while(y)
        {
                b.pb(y%10+'0');
                y/=10;
        }
        while(x || a.length()!=b.length())
        {
                a.pb(x%10+'0');
                x/=10;
        }
        reverse(all(a)),reverse(all(b));
        for(int len=10;len>=1;len--)
        {
                memset(dp,-1,sizeof(dp));
                if(count(0,-1,len,0,0,0))
                {
                        cout<<len<<" "<<count(0,-1,len,0,0,0)<<endl;
                        return;
                }
        }
}

SourabhCh

dp[pos][lastDigit+1][need][flagA][flagB][flagZero]=res;

why lastDigit+1 ?

Cuz i'm theoretically placing -1 at position (-1) as a dummy digit as it automatically considers all cases. So, 0 points to -1 so offset of 1 is used.

coco_elon

← Rev. 3 →

Some additional problems:-

Chef and special numbers

Target2018

How to solve "palindromic numbers"?

egor.okhterov

Initially we have the space of all numbers S = {0, 1, 2, 3, ..., 10¹⁷}.
This space is too big, so we want to reduce it.
From all that space we only need to consider the numbers with a special property (being palindrome) P ⊂ S.

Now let's consider another space where all of our original numbers are mapped into a different number S₂ = {f(0), f(1), f(2), f(3), ..., f(10¹⁷)}.

There is no point in cosidering this new space S₂ instead of original S if it has the same size.
Can you think of any mapping f: S⟶S₂ that will make the size of S₂ manageable?

harshit2202

Can't think of solution. Any more suggestions?

SajidZakaria

Vaiya can you please sort the problems according to difficulty?

← Rev. 4 →

I didn't solve them all but here goes my order

Investigation
Sum of Digits
Digit Sum
Bomb
Ra-One Numbers
LUCIFER Number
G-One Numbers
Digit Count
Round Numbers
Fast Bit Calculations
Magic Numbers
Maximum Product
If you need help then you can look at this.

MotaSanyal

In the problem "Investigation" whats the logic behind taking the max remainder as 90 only? you have defined the dp states as dp[90][2][90][90]. why 90 ?Target2018

kenechi

Since the sum of any number less than 2^31 will be less than 83(approx) so any k > 83 wont be able to divide the sum of numbers . So for k > 83 answer should always be 0. That's why we can take max remainder value as 90.

My sol link :- Solution

Cyber_Wizard

I have the same question The sum of digits of any number less than 2^31 will be less than 90 but the number itself should also be divisible by k So I think the remainder can be greater than 90 and it should be less than 10000 Can you please explain?

Since we are not checking the answer for k > 83 in the question because the answer for k > 83 will always be 0 as I explained above so if while creating a digit we are taking its mod with k (which is always less than 83) the max value of that state won't cross 90.

Deba

You can add this problem

Merripium

919B - Perfect Number also can be solved by digit DP and it's very amazing if n not greater 10^6

Edward_Haddad

5 months ago, # ^ |

Can you send the solution?

bounce_back20

how to solve http://mirror.codeforces.com/gym/100886/problem/G ? please tell me the state of the dp

Take a look

thanks !

-10

i have seen the solution , can you please explain the use of FLAG in dp state , how its helping?

BumbleBee

Which way is more efficient? Double recursion call with 3 states or single recursion call with 4 states?

ILoveBitches

Can Anyone please provide solution for this tutorial also ? I am noob. Sorry for that.

sinnersneversleep

I am having really hard time thinking about the time complexity of this approach can anyone please help

6 years ago, # |

what states are used in solving the problem(Cantor)

Gr8

Is Classy Numbers also comes in this category?

zoombee

yes

It was the simplest problem of DIGIT DP.

Testduk

could you explain the solution

dp states for Classy numbers I used: index,prev_sum

here prev_sum denotes the number of non zero digits. check out my submission : Solution

TheRedLegend

Here is my solution using digit DP 42675923, it may help you. There are quite a lot of cases to take care of.

The main idea is that for each interval [L, R] the answer is F(R) - F(L - 1). Where F(x) equals the number of classy integers from 0 to x.

How to you calculate F(x)? Just using Dynamic Programming.

For a number x, the DP looks something like this:

dp[i][j][k][0] = number of integers less than x that are i digits long, end in j and they have k digits bigger than 0.

dp[i][j][k][1] = number of integers equal to x that are i digits long, end in j and they have k digits bigger than 0.

j is going to be between 0..9 and k in the range from 0..3, because a classy number contains no more than 3 non-zero digits.

My dynamic has more states than the ILoveBitches's but I think it's more intuitive. Hope it's going to help you.

https://www.hackerrank.com/topics/digit-dp this is a good tutorial

dreamplayDaddy

Easy to understand solution using the same template as described in the post.
Hope it helps :)
solution

ay2306

Another really good question on digit dp: Link

prodipdatta7

Another good problem in digit dp Balanced Number

CuteCoder07

Why the for loop start from 0? let a = 100, d = 5, k = 1

Please correct me if I'm wrong, because one of the base cases is pos==num.size(), so there will will be state with, 050 chosen right? this is representing number 50? and is that 005 means it's number 5?, because we return 1 when pos == num.size() so we have to fill each digit, this zero a bit difficult to understand.

Yes, you are correct. 005 means 5 and 050 means 50. So the way i have written my code, it'll consider all the numbers <= N including the ones whose number of digits is less than the number of digits in N. So it'll work perfectly.

But when d = 0, my given code will produce an incorrect output. Because it'll count zero for the numbers who have leading zeroes. For example, for the number 005, it'll count zero two times. Which is wrong. Because we don't consider the leading zeroes in the decimal numbers. You have to handle this case separately.

How can we do that? We can use another boolean state, which will keep track if we have placed at least one non zero digits so far. If it's true in any state, then we'll count the zero in that state. Otherwise, we don't. Hope it helps you to understand.

whoami_shubham

Thanks bro :) .

flash_7 can you please tell me how to solve problem
6.Maximum Product
by digitDp as i am not able to get what should be the states though i have already solved this question using greedy approach solution but i want to know how it can be solved using digitDp or more precisely is it solvable using digitDp approach or not?

In question Investigation I have been using dp states as pos, flag, sum_of_digits, remainder_of_so_far_number.

So the dp size should be dp[32][2][10000][10000]
but this gives RTE and dp[32][2][101][101] passes.

But when k>100 this should give the wrong answer but it is passing My submission https://vjudge.net/solution/18823741 I can't understand why such dp size works?

SAKIB_ALAMIN

just think about , what could be maximum sum of digits & does the value of k matters if it's greater than the sum of digits ?

Sum of digits maximum can be 100 only I know But the remainder_of_so_far_number can be greater than 100 upto value of k i.e.10000

so , if sod( sum of digits ) is < 100 & k >= 100 , these numbers will never be divisible by k according to this problem , right ? so , if k is >= 100 , we can simply print 0 , otherwise use dp , thus state gets reduced.

Thanks:) Got it

great :-)

What should be the dp state in Palindromic Numbers?

vamaddur

Problem Editorial

ssanjeev2424

5 years ago, # |

Thanks a lot, it really helped !! @flash_7

lapmid

flash_7 How To Solve Palindromic Numbers ??

sabertooth

It May Help you

nilesh8757

Add ENCODING to the list.

sumit083

Can someone help in the question that was asked in the long challenge of CodeChef in August 2019. Problem Code -ENCODING LINK — https://www.codechef.com/AUG19B/problems/ENCODING

babin

Hello! There is a very detailed tutorial for this problem on Codechef. If you have further questions, feel free to ask(I was the setter).

siisback

A good problem for Digit DP: https://mirror.codeforces.com/contest/288/problem/E

k_8

Thanks for this Article. This following problem is good to be added. Playing with Digits(Hackerearth)

ronaldocr7

can anyone help me with DIGIT COUNT I am not getting states for my dp array please help

Haunted_Cpp

dp[index][last_digit_placed][flag]

xodeChemist

thank you so much this is really helpful my digitdp journy has started

mike51

Really good article about digit DP.

ritik_patel05

E — Almost Everywhere Zero, AtcoderBegCon — 154

Similar question to the one discussed in the article. Thank you a lot for the very well explained article!

zarif_2002

This is also a nice problem as far as I think. Solved now after learning this topic right now. https://mirror.codeforces.com/problemset/problem/1036/C

cenation092

I have made a video on digit DP. Digit Dp Link

jiteshrao

4 years ago, # |

How to solve this problem : calculate the no of times digit d occurs in l to r? sample input : l = 1, r = 22, d = 2 sample output : 6 (2, 12, 20, 21, 22) frequency of 2 is 6. I know a solution having complexty (len(n) * len(n)).

same_lame_name

Can you send the problem link? I want to try my idea.

Based on just what you have mentioned, I have a sketch of a solution. Here

I'm curious whether this is correct. It could require some tweaks, as for case when d = 0, It also counts the leading zeros.

Thanks for your reply.I also did in the same way.

rogernadal

If someone prefers videos, here are few problems and their solutions :

Problems (approx 13-14) with solution : https://www.youtube.com/watch?v=Hsnbaixm1A4

Theory : https://www.youtube.com/watch?v=5CLtoXus7O4

cmlivb

Are questions given in the end are sorted according to the increasing level of difficulty ??

ChitreshApte

Thank You for the informative article flash_7

I_Killed_Jason_Todd

Investigation
Can somebody tell me the dp state for the above problem?

My dp state is dp[pos][rem][f][rem2] = count of numbers when we are at position pos with a remainder rem(sum of digits mod K) and with a remainder rem2(the actual number mod K) The f in my dp state tells us whether the number has already become smaller than the original number or not. Java Code

Although the approach is right, my solution got Memory Limit Exceeded. Any ideas for optimizations?

someguy

Can anyone tell what I'm doing wrong in the question 1-Investigation, code link- https://pastebin.com/sd44xj20

the0dd1out

someguy you initialized 'tight' with 0 instead of 1 in your dp function!

ND_

Is any editorial available for these problems? If yes please paste the link here. Sorry for Spamming.

kushb30

3 years ago, # |

flash_7 How to solve maximum product using digit dp .Here is my solution but don't know where it is going wrong

hardikkapoor

2 years ago, # |

what is dp storing really? I'm having trouble understanding that.

vantaablackk

11.Sum of Digits and 12. Digit Sum are basically the same problem with different input style. I suggest removing Sum of Digits as it's statement is also wrong (asks for range [1,n] but it's actually [a,b]) and it has weird input format. Without changing the ordering, you can substitute CSES Counting Numbers in its place.

MrCatsworth

Palindromic numbers. Few things to mention.

5 states in dp,

dp[pos][len][last_digit][smaller][zerolead].

1)pos — where we are? MAXN = 2e1.

2)len — current maximum len without leading zeros.Suppose we have some cases like 00121,where len = 3. MAXN = 2e1.

3)last_digit- what kind of digit does store the result? MAXM = 10.

4)smaller — has the sequence already become smaller than x? where x is i or j. MAXF = 2.

5)zerolead — have we already placed our first digit that is not 0? MAXF = 2.

6 states in rec.Let's say rec is

ll rec(int pos,int len,int last_digit,ll cur_num,bool smaller,bool zerolead).

6)cur_num to continue building our initial number after (pos + 1) > (len + 1) / 2* with just trying to give symmetrical value for O(1),

but before that,we could give any value for O(10).

ll pw[],where we store power of 10 to get the value from our cur_num for O(1) and put it in order to follow the rules after (len + 1)/2.

(P.s.*pos starts from zero).

Good luck,have fun.

_harsh_upadhayay_

21 month(s) ago, # |

Yet another digit DP problem : Unlucky Numbers.

RISHI_SHARMA

can anyone please help with the base condition in the sample code. i am unable to understand if(pos==num.size()) { (if(cnt==k) return 1; return 0; }