Skip to main content

Find the power of each factor in a given number using matlab function

Problem

Write a function in matlab to find the exponent of each prime factor of a given number. For example if the input is 28 then the output should be


Power of 2 in the number 28 is 2
Power of 3 in the number 28 is 0
Power of 5 in the number 28 is 0
Power of 7 in the number 28 is 1

Solution

The solution is simple and one has to know the function definition, factor function, for loop, nnz function, isprime function, and the if statement to understand the program. Don't worry if you don't know the function definition. You can even program this code without a function. The difference will be that there will not be a function definition line, and the value of n is directly defined instead. But it would be useful if you know them. If you want to learn any of the above then head on to the links below.

Find the explanation for the program below the program

Program


function ret = countfacpowers(n)
%this function will calculate the exponent of each prime factor and return
%the result as a printed statement

vec = factor(n);
highest = vec(length(vec));

for i = 1:highest
    if isprime(i) == 1
        fprintf('Power of %d in the number %d is %d\n',i,n,nnz(vec == i));
    end
end
To get the matlab for the above code click here: countfacpowers.m

Explanation

During the starting of the function a vector is created to store the prime factors of the given numbers using the factor function. Now digging deep a highest variable is created in order to store the largest prime factor of the number. Here you can also use the max function but not to complicate this program more for the novice programmers, I just have used the last number of the vector which is the largest prime factor. Now let us understand the for loop iteration by iteration. I am considering the input number as 28 for this example

First iteration

During the first iteration the value of i is 1 and this doesn't qualify the if statement and doesn't execute. Thus first iteration is over and the second iteration starts

Second iteration

During the second iteration the value of i is changed from 1 to 2 and the value of 2 qualifies the if statement and the fprintf statement gets executed. Here we will understand what is happening in the nnz function. The nnz function in general gives the number of non zero elements in the vector. But here we have given the input to the nnz function as the vector whose values are 1 only where the value becomes equal to i. Let us simplify this:
First the value of the vec is [2 2 7]
Now we will call vec == 2 and this will give the output as [1 1 0] , 1 will only be created at the places where there is 2(Here in this case 2 and where there is i in general). Finally nnz([1 1 0]) will give 2 i.e. 1+1 = 2.
and the second iteration comes to halt

Similarly all the other iterations continue until the i value changes to 7.

I have tried to explain this program as easy as possible, but if you don't understand then please do contact me so that many others will get benefited 

The above code was highlighted using hilite.me

To read the above post in pdf format then you can do so by downloading from here: countfacpowers.pdf

Popular posts from this blog

Making a quiz web app with python and flask

Edit : When you are creating a web app with h tml templates, then y ou will have to sa ve the html file in templates folder in the Current Wor ki ng Directory( CWD). If you save the file in the C W D directl y you will get a TemplateNotFound error. Thank you Udhay for pointing it out.   In this post we will create a quiz website using python . I will be using the flask framework . After reading this tutorial you will learn form submission , flask templates , python code in flask templates , shuffling the questions and options with the random module and few others.  Please note that this tutorial is not big as it seems to be. Some of the code has been rewritten to maintain consistency and also font size is somewhat big so that your eyes won't get stressed reading this tutorial. Also the content has not occupied the full width of the page. In this tutorial I am assuming that you are having a very basic understanding of the flask framework . Please refer the documentation

Problem 11 Project Euler Solution with python

Largest product in a grid In the 20×20 grid below, four numbers along a diagonal line have been marked in red. 08 02 22 97 38 15 00 40 00 75 04 05 07 78 52 12 50 77 91 08 49 49 99 40 17 81 18 57 60 87 17 40 98 43 69 48 04 56 62 00 81 49 31 73 55 79 14 29 93 71 40 67 53 88 30 03 49 13 36 65 52 70 95 23 04 60 11 42 69 24 68 56 01 32 56 71 37 02 36 91 22 31 16 71 51 67 63 89 41 92 36 54 22 40 40 28 66 33 13 80 24 47 32 60 99 03 45 02 44 75 33 53 78 36 84 20 35 17 12 50 32 98 81 28 64 23 67 10 26 38 40 67 59 54 70 66 18 38 64 70 67 26 20 68 02 62 12 20 95 63 94 39 63 08 40 91 66 49 94 21 24 55 58 05 66 73 99 26 97 17 78 78 96 83 14 88 34 89 63 72 21 36 23 09 75 00 76 44 20 45 35 14 00 61 33 97 34 31 33 95 78 17 53 28 22 75 31 67 15 94 03 80 04 62 16 14 09 53 56 92 16 39 05 42 96 35 31 47 55 58 88 24 00 17 54 24 36 29 85 57 86 56 00 48 35 71 89 07 05 44 44 37 44 60 21 58 51 54 17 58 19 80 81 68 05 94 47 69 28 73 92 13 86 52 17 77 04 89 55 40 04 52 08 83 97 35 99 16 07

Problem 60 Project Euler Solution with python

Prime pair sets The primes 3, 7, 109, and 673, are quite remarkable. By taking any two primes and concatenating them in any order the result will always be prime. For example, taking 7 and 109, both 7109 and 1097 are prime. The sum of these four primes, 792, represents the lowest sum for a set of four primes with this property. Find the lowest sum for a set of five primes for which any two primes concatenate to produce another prime. This problem is j u st a brute force problem. If you have come here because you don't know the limit upto which you will h ave to gener ate the prime numbers t hen go ahe ad and t r y with 10,000 . When I first start ed solving the problem I chose 1 million(beca use most of the problem s on project E uler have this limit ), but it took very long for the computer to fin d the solution. After searching on the internet then I found many people choosing 10, 000 so I have changed my in put f rom 1 million to 10000 and the output was f ast. He