  The phase transition in the q-state Potts model with homogeneous
ferromagnetic couplings is strongly first order for large q, while is rounded
in the presence of quenched disorder. Here we study this phenomenon on
different two-dimensional lattices by using the fact that the partition
function of the model is dominated by a single diagram of the high-temperature
expansion, which is calculated by an efficient combinatorial optimization
algorithm. For a given finite sample with discrete randomness the free energy
is a pice-wise linear function of the temperature, which is rounded after
averaging, however the discontinuity of the internal energy at the transition
point (i.e. the latent heat) stays finite even in the thermodynamic limit. For
a continuous disorder, instead, the latent heat vanishes. At the phase
transition point the dominant diagram percolates and the total magnetic moment
is related to the size of the percolating cluster. Its fractal dimension is
found d_f=(5+\sqrt{5})/4 and it is independent of the type of the lattice and
the form of disorder. We argue that the critical behavior is exclusively
determined by disorder and the corresponding fixed point is the isotropic
version of the so called infinite randomness fixed point, which is realized in
random quantum spin chains. From this mapping we conjecture the values of the
critical exponents as \beta=2-d_f, \beta_s=1/2 and \nu=1.
