the large magellanic cloud ( lmc ) has revealed a very complex structure both in the stellar and in the gaseous component .
the elongation of the stellar disk in the direction of the galactic center , its substantial vertical thickness , the warp and the strong asymmetric bar are naturally predicted by numerical simulations as a result of the gravitational interaction between the lmc and the galaxy ( bekki & chiba 2005 , mastropietro et al .
2005 , hereafter m05 ) .
the old stellar distribution appears to be quite smooth in the outer parts of the disk , with no signs of spiral structures out to a radius of 10 kpc @xcite . within the same radius
the hi large scale structure reveals the presence of several asymmetric features that have no equivalent in the old stellar disk .
the gaseous disk is characterized by the presence of an elongated region located at the south - east of the galaxy and aligned with the border of the optical disk , where the column density distribution shows a steep increase @xcite .
since the lmc proper motion vector is directed to the east , it appears natural to associate this high density region with ram - pressure acting on the leading edge of the disk due to the orbital motion of the lmc and its consequent interaction with the diffuse hot gas in the halo of the milky way ( mw ) .
the presence of an extended hot halo surrounding galaxies and in hydrostatic equilibrium within the dark matter potential is expected by current models of hierarchical structure formation . in the mw ,
x - ray absorption lines produced by hot ( @xmath1 k ) gas are detected in the spectra of several bright agn @xcite .
some ionization features discovered in the magellanic stream and high velocity clouds indicate that this distribution of hot gas extends well beyond the galactic disk ( @xmath2 kpc ) .
constraints from dynamical and thermal arguments fix its density in a range between @xmath3 and @xmath4 cm @xmath5 at the lmc distance from the galactic center ( but kaufmann et al .
2009 suggest a value ten times higher ) .
@xcite have performed a detailed analysis of the lmc global star formation rate using asymptotic giant branch stars .
they find an irregular and patchy distribution in age , with the youngest carbon - rich systems located at the south - east of the disk .
the present star formation activity is rather clumpy and concentrated in stellar complexes characterized by intense hii emission and associated with bright h@xmath0 filamentary bubbles .
most of these very young structures lie on the south - east of the disk , in the proximity of 30 doradus , the largest star forming region of the lmc , some are located in the bar and the remainder form an asymmetric pattern that covers the entire disk with no apparent relation to the global geometry of the satellite .
it is not clear which is the overall physical mechanism responsible for triggering star formation with the observed asymmetric pattern and different models have been proposed in the past .
the stochastic self - propagating star formation ( sspsf ) model predicts a clear age gradient in the lmc s stellar complexes , with the edges being younger with respect to the center , in contradiction with observations @xcite .
@xcite proposed a scenario where the bow shock originated by the motion of the lmc through the hot galactic halo compresses the leading edge of the disk and induces star formation .
the pressure at the south - eastern edge of the lmc is indeed 10 times higher than the average in the rest of the lmc @xcite .
this model , which assumes the orbital motion vector lying in the plane of the disk , predicts increasing ages of the stellar complexes in the direction of the rotation , due to the fact that the material compressed at the front side of the disk moves , in time , away to the side .
the youngest systems would indeed lie at the south - east border of the disk , where the relative velocity between the corotating interstellar medium and the external diffuse gas is maximum .
several giant structures along the outer east and north edge of the lmc actually show a progression in age in a clockwise direction : moving from south - east to the north lmc 2 , 30 doradus and lmc 3 , lmc 4 , ngc1818 .
in particular the difference in age between 30 doradus and lmc4 is exactly their distance along the border of the disk divided by the satellite s rotational velocity ( harris , private communication ) .
studied the recent star formation history of the lmc using cepheids and other supergiant stars and found that although the majority of the star formation events in the last 30 myr are concentrated on the east border , others are distributed across the entire disk in partial contrast with the bow shock induced star formation model , that can not explain them .
in this work we use high resolution sph simulations to study the effects of the interaction between the lmc interstellar medium and the diffuse hot halo of the mw .
we want to investigate whether the ram - pressure acting on the leading edge of the lmc disk is responsible for the increase in density observed in the south - east and for triggering star formation .
the analytic model of @xcite assumes a pure edge - on model , but according to @xcite the present angle between the lmc disk and the orbital motion is nearly @xmath6 . even in the absence of precession and nutation , this angle is subjected to large variations during the orbital period in such a way that compression produced by the external hot gas can affect in time both edge - on and face - on .    moreover , the ram - pressure felt by the lmc is not constant and has a maximum when the satellite approaches the perigalacticon .
the motion of the lmc through the hot halo of the mw during the last 1 gyr is modeled using `` test wind tunnel '' simulations with increasing ram - pressure values .
the paper is structured as follows .
section 2 describes the models and the star formation criteria adopted , section 3 illustrates the results of simulations without star formation , focussed on the investigation of pure effects of compression on the lmc interstellar medium while section 4 describes the runs where star formation is activated .
several simulations have been performed , assuming different star formation models , disk inclinations and hot halo densities .
the initial conditions of the simulations are constructed using the technique described by @xcite .
our disk galaxy model is a multi - component system with a stellar and gaseous disk embedded in a spherical nfw @xcite dark matter halo .
the density profile of the dark matter halo is adiabatically contracted in response to baryonic infall @xcite .
the stellar disk follows an exponential surface density profile of the form : @xmath7 where @xmath8 and @xmath9 are the disk mass and radial scale length ( in cylindrical coordinates ) , respectively , while the thin vertical structure has a scale length @xmath10 : @xmath11 the gaseous disk is characterized by an exponential profile with the same radial and vertical scale length as the stellar component and by a constant density layer which extends up to 8@xmath9 .
the structural parameters of the disk and the halo are chosen so that the resulting rotation curve resembles that of a typical bulgeless late - type ( sc / sd ) disk galaxy .
they are similar to those adopted in m05 for the initial lmc model and reproduce quite well the peak of the rotation curve inferred by @xcite ( fig .
as seen in m05 , the interaction with the mw does not affect significantly the stellar and dark matter mass in the inner @xmath12 kpc of the lmc and consequently the global rotation curve within this radial range .
the choice of an extended gaseous component for the initial lmc model is motivated by the fact that spiral galaxies in the local universe are commonly observed to be embedded in extended disks of neutral hydrogen significantly larger than their stellar component . as seen in m05 ,
the combined effect of tidal interactions and ram - pressure stripping can remove a significant fraction of gas from a lmc disk orbiting within the hot halo of the mw , with a ram - pressure stripping radius which is a factor of three smaller than the initial radius of the gaseous disk .
also in the case of a lmc with orbital velocities significantly higher @xcite hydrodynamic and gravitational forces together are effective in resizing and reshaping the extended gaseous disk of the satellite beyond 8 kpc . in the present work we neglect the presence of gravitational forces focusing on the effects of pure ram - pressure .
therefore we do not expect to see a significant decrease in the radius of the gas distribution .
however , in order to take in account the loss of cold gas from the disk of the satellite and the star formation events , we assumed an initial amount of gas in the disk which is about @xmath13 times larger than the hi mass in the lmc ( @xmath14 @xmath15 according to putman et al .
2003 ) .
the mass within the virial radius is set equal to @xmath16 m@xmath17 and the fraction of mass in the disk is @xmath18 , equally distributed between the gaseous and stellar component . the contribution of the different components to the global rotation curve , assuming a disk scale length @xmath19 kpc and a dark halo concentration @xmath20 ( where @xmath21 is defined as @xmath22 , with @xmath23 and @xmath24 virial and scale radius of the nfw halo , respectively ) is plotted in fig .
[ rotcurve ] .
the halo spin parameter , which sets the disk scale length in our modeling , is @xmath25 , where @xmath26 relates the angular momentum @xmath27 and the total energy @xmath28 of a system with virial mass @xmath29 through the relation @xmath30    the initial stellar disk of the satellite galaxy has , within its scale radius @xmath9 , a central mass surface density of @xmath31 ( fig .
[ diskprofile ] ) , that corresponds to a b - band surface brightness of @xmath32 , assuming a mass to light ratio @xmath33 .
the central gas surface density is only @xmath34 since a significant fraction of gas is distributed in the external disk .
assuming @xmath35 hi abundance this value corresponds to an hydrogen column density of @xmath36 within @xmath9 , comparable with the values observed by @xcite with the lmc parkes multibeam hi survey .      in order to obtain a strongly stable disk against bar formation even in the presence of significant gas stripping and consequent perturbation of the satellite potential , the thickness of the stellar component is set such that the toomre s @xcite stability criterion is largely satisfied .
in particular the toomre s parameter for the stellar disk :    @xmath37    where @xmath38 is the radial velocity dispersion , @xmath39 is the local epicyclic frequency and @xmath40 the unperturbed stellar surface density , has a minimum at the disk scale length with @xmath41 . for a gaseous disk the stability of the disk
is expressed in terms of the gas sound speed @xmath42 and surface density @xmath43 through the relation : @xmath44 the gaseous disk has initially a constant temperature of 10000 k , which implies @xmath45 and @xmath46 .    according to and @xcite , the stability of a multicomponent disk
is not guaranteed by the individual stability of its single constituents , due to the mutual gravitational interaction between gas and stars .
stars are characterized by velocity dispersions 3 - 4 times larger than the typical sound velocities in the cold gaseous disk and even relatively small variations of the gaseous component parameters can significantly affect the stability of the whole disk .
therefore we choose a large value of @xmath47 to contrast the effects of ram - pressure . in the case of a two components - gaseous and stellar - disk , the stability condition
is expressed by @xmath48^{-1}>1 , \label{toomretotal}\ ] ] where @xmath49 and @xmath50 .
[ qparameter ] illustrates the dependence of @xmath51 on the dimensionless wavenumber of the perturbation @xmath52 within three different regions of the disk : at the disk scale length @xmath9 , at @xmath53 kpc and in the external region ( @xmath54 kpc ) , where the gas component predominates .
the criterium is always satisfied and the disk is stable against axisymmetric perturbations , independently of their wavelength . in order to check stability the disk
was initially evolved in isolation for 1 gyr .
= 8truecm    = 8truecm    = 8truecm    all the simulations we now discuss were carried out using gasoline , a parallel tree - code with multi - stepping @xcite which is an extension of the pure n - body gravity code pkdgrav developed by @xcite .
the code uses a spline kernel with compact support where the interaction distance for a particle @xmath55 is set equal to two times the smoothing length @xmath56 , defined as the @xmath57-th neighbour distance from the particle . in this paper @xmath58 .
the internal energy of the gas is integrated using the antisymmetric formulation of @xcite that conserves entropy closely .
dissipation in shocks is modeled using the quadratic term of the standard @xcite artificial viscosity .
the balsara @xcite correction term is used to suppress the viscosity in non - shocking , shearing environments .
the code includes radiative cooling for a primordial mixture of hydrogen and helium in collisional equilibrium . at temperatures below @xmath59
k the gas is entirely neutral and due to the lack of molecular cooling and metals , the efficiency of the cooling function drops rapidly to zero .
we used a star formation recipe that includes density and temperature criteria , while converging flow criterium is not required in most of the simulations .
gas particles are eligible to form stars only if the density of the star formation region has a minimum physical density corresponding to 0.1 hydrogen atoms per @xmath60 @xcite and an overdensity @xmath61 ( katz et al . 1996 ) , which basically restricts star formation to collapsed , virialized regions . the physical density threshold describes the steep drop in star formation rate observed in disk galaxies when the gas surface density is much lower than a critical value @xmath62 @xcite .
the density threshold @xmath63 cm@xmath5 is compatible with observational results . according to @xcite ,
the star forming region has to be part of a converging flow that implies a local negative divergence of the sph velocity field .
however , the converging flow criterium was introduced to describe star formation in cosmological simulations , where the geometry of the collapsing regions is approximatively spherical . in the case of star formation regions like 30 doradus , localized at the periphery of the lmc disk where gas particles relatively close in distance can have significantly different kinematics
, this criterium leads to underestimate of the star formation rate .
therefore the converging flow is not required in most of the simulations . a single run including
the converging flow requirement has been performed for comparison .
the star formation rate is assumed to be proportional to @xmath65 @xcite , where @xmath66 represents the volume density of the cold gas , and is given by the expression @xcite @xmath67 where the star formation timescale @xmath68 is the maximum between the local gas dynamical collapse time @xmath69 and the local cooling time .
if the gas is already cool enough to form stars i.e @xmath70 , then @xmath71 is used .
we assumed @xmath72k .
the constant star formation rate parameter @xmath73 is chosen such that we reproduce the global lmc star formation rate @xcite .
once a gas particle satisfies the above criteria , it spawns stars according to a probability distribution function .
in particular , the probability @xmath74 that a gas particle forms stars in a time @xmath75 is modeled as @xmath76 a random number is then drawn to determine whether the gas particle forms stars during @xmath75 .
for all the simulations in this paper @xmath77 myr .
the newly created collisionless particle has the same position , velocity and softening length as the original gas particle while its mass is a fixed fraction @xmath78 of the parent gas particle , whose mass is reduced accordingly . following @xcite we assumed for our favorite models a star formation efficiency @xmath79 .
up to six particles are then created for each gas particle in the disk . after its mass has decreased below @xmath80 of its initial value the gas particle is removed and its mass is re - allocated among the neighbouring particles .      in order to study the influence of pure ram pressure on a galaxy model orbiting in a milky way halo , we performed `` wind tunnel '' simulations where the ram - pressure value varies with time .
we represent the hot gas as a flux of particles moving along the major axis of an oblong of base equal to the diameter of the dark matter halo of the satellite and height @xmath81 , where @xmath82 is the velocity of the lmc at the perigalacticon and @xmath83 is the time scale of the simulation .
the hot particles have an initial random distribution and a temperature @xmath84 .
the box has periodic boundary conditions in order to restore the flow of hot gas that leaves the oblong .
the galaxy model is at rest at the center of the oblong .
kinematical data @xcite indicate that the lmc , presently located at @xmath85 kpc from the galactic center , is just past a perigalactic passage and has an orbital velocity of about 300 km s@xmath86 .
recent proper motion measurements by @xcite and @xcite suggest that the velocity of the satellite is substantially higher ( almost 100 km s@xmath86 ) than previously estimated and consistent with the hypothesis of a first passage about the mw @xcite . in both scenarios
the cloud is affected by the largest ram - pressure values during the last million years of its orbital evolution . indeed , while in the models proposed by @xcite the lmc does not enter the halo of the mw earlier than 1 gyr ago ( slightly different orbits are found in mastropietro 2008 ) , in m05 we have shown that the change in the orbital parameters due to dynamical friction strongly affects the ram - pressure stripping rate . even in the case of a `` low velocity '' model the largest ram - pressure on the satellite is expected during the last orbital semi - period ( about 1 gyr ) due to the increasing velocity and external gas density .
we followed ram - pressure acting on the lmc s igm during the past 1 gyr . the density of the hot external gas increases with time , in such a way that in our low velocity model the external pressure experienced by the cold disk varies from @xmath87 to @xmath88 at the time corresponding to the pericentric passage .
this is equivalent of assuming @xmath89 and @xmath90 , and a number density of the external gas that increases from @xmath91 to @xmath92 at @xmath85 kpc from the galactic center .
these density values are comparable with those provided by m05 , who modeled the mw hot halo assuming a spherical distribution of gas that traces the dark matter profile , with a mean number density of @xmath93 within 150 kpc .
we also consider the eventuality of a less dense galactic halo and performed runs where the gas density is a factor 10 lower .
models with higher velocities @xcite and orbital parameters similar to those suggested by @xcite are also explored . in details
, @xmath94 is the same as in the low velocity models since the higher orbital velocity at the beginning of the simulation ( about 250 km s@xmath86 ) is compensated by a lower external density ( @xmath95 , according to mastropietro 2008 . indeed 1 gyr ago the satellite has just passed through the virial radius of the mw ) .
the maximum pressure felt by the disk is @xmath96 , that corresponds to a cloud moving with @xmath97 through an external hot medium of density @xmath92 .
each galaxy model is simulated using 750000 particles , of which @xmath98 are in the dark halo and @xmath99 in the disk ( @xmath100 collisional and @xmath101 collisionless ) .
the hot gas in the `` wind tunnel '' has @xmath102 particles , in such a way that the mass ratio @xmath103 between hot particles and particles in the disk is close to the unity even when the halo density is the largest .
this choice permits to avoid the presence of scattering and numerical holes which artificially change the shape of the front edge and influence the morphology of the disk ( m05 ) .
the gravitational spline softening is set equal to 0.5 kpc for the dark halo and the hot gas in the oblong , while it is 0.2 kpc for stars and gas in the disk .
in order to study the effect of pure compression on the density distribution of cold gas in the lmc disk , we have run a first set of simulations where the gas cools radiatively but star formation was not activated .    according to @xcite and @xcite , the present angle between the lmc s disk and its proper motion vector
is roughly @xmath6 . even neglecting the effects of precession and nutation on the disk plane of the satellite @xcite this angle
is expected to vary significantly during an orbital period , especially in the proximity of a pericentric passage due to rapid changes in the velocity vector .
different relative orientations of the disk with respect to the orbital motion are therefore investigated .
the inclination angle @xmath55 is defined as the angle between the angular momentum vector of the disk and the flux of hot particles in the wind tube , so that a galaxy moving edge - on through the external medium is characterized by @xmath104 , while the observed lmc disk would have @xmath105 .
we explored cases with inclination angle @xmath55 of 90 , 45 and @xmath106 ( runs cool90 , cool45 and cool10 , respectively ) .    with a hot halo temperature of @xmath107 k
the relative velocity between the satellite and the external medium is supersonic ( sound speed @xmath108 km s@xmath86 and mach number @xmath109 and @xmath13 at the pericenter of the low and high velocity orbit , respectively ) and a bow shock forms in front of the disk ( fig .
[ shockimage ] ) .
since the cooling time of the post - shock gas is @xmath110 gyr , the shock can be considered adiabatic and hydrodynamical quantities at the two sides of the shock front are in first approximation related by the rankine - hugoniot equations for a stationary normal shock .
= 8truecm    for @xmath111 the jump conditions give @xmath112 and @xmath113 , where subscripts 1 and 2 denote upstream and downstream quantities .
the ram - pressure @xmath114 actually felt by the galaxy behind the shock front is therefore smaller than that it would suffer due to the upstream flux of hot particles , but conservation of momentum flux across the shock discontinuity implies that the reduction in dynamical pressure has to be balanced by an increase in thermal pressure ( see also rasmussen et al .
2006 ) .    fig .
[ shock ] illustrates the behavior of hydrodynamical quantities across the shock discontinuity for a snapshot corresponding to the perigalacticon of a low velocity orbit .
the disk inclination is @xmath115 .
the horizontal axis is centered on the lmc stellar disk and oriented perpendicularly to the bow shock nose , with the shock located at @xmath116 kpc and the satellite moving towards increasing values of @xmath117 .
the mach number derived by the temperature jump is @xmath118 , in good agreement with the theoretical value for a normal shock .
only hot halo particles are considered in computations but , due to the sph nature of the simulations , close to the border of the disk we observe a further density increase and a sharp drop in temperature .
the @xmath117-velocity profile is plotted in the system of reference where the pre - shock gas is at rest .
= 8truecm    pressure profiles across the shock are plotted in fig .
[ shockpressure ] , where @xmath119 is defined as @xmath120 .
the total pressure remains roughly constant until the edge of the disk .
the steep increase at @xmath121 kpc is due the rapid growth in density at the border of the disk , which is not immediately followed by a decrease in temperature .
for @xmath122 the shock wave is inclined with respect to the initial flow velocity and the rankine - hugoniot conditions apply to the normal components of the velocity across the shock discontinuity , while the component parallel to the shock front remain unchanged .
the flow is therefore deflected toward an oblique shock wave and the jump at the shock discontinuity is smaller .
an edge - on ( @xmath123 ) disk behaves like a wedge moving supersonically with the vertex facing upstream .
if the wedge angle is smaller than or equal to the maximum flow deflection angle , the oblique shock becomes attached to the vertex of the wedge and the flow is deflected so that the streamlines are parallel to the surfaces of the wedge .
the shock standoff distance thus depends on the external density profile of the collisional edge - on disk and on the mach number of the incident flow . in both low and high velocity edge - on models the shock results almost attached to the disk
.    figs .
[ densitymaps][densitymaps10 ] illustrate  for different values of the inclination angle @xmath55  changes in the disk gas density distribution as the satellite passes through increasing values of the external pressure , moving towards the perigalacticon .
each couple of panels illustrates the state of the disk at increasing ( from the top to the bottom ) times along the orbit .
panels on the left represent hi column density maps .
the density contrast is chosen in order to highlight the density gradient in the external disk , since the gas distribution in the central regions of the lmc is dominated by the presence of the bar and a direct comparison with pure ram - pressure simulations is not possible .
the color scale is logarithmic , with white corresponding to a density larger than @xmath124 @xmath125 and blue to values lower than @xmath126 @xmath125 .
hot gas particles flow on to the disk from the left to the right of each plot , with increasing ram - pressure values from the first to the fourth image .
the disk is seen face on and rotating clockwise . in the case of the edge - on run ( fig .
[ densitymaps ] ) cold gas particles lying in the left - bottom quadrant of the disk feel the largest ram - pressure , due to the fact that their relative velocity with respect to the external medium is maximum .
the rotational velocity of the external disk is @xmath127 km s@xmath86 , which implies a relative velocity at the pericenter of @xmath128 km s@xmath86 ( 450 km s@xmath86 in the case of a high velocity orbit ) .
panels on the right represent the change in mean density and radius of the gaseous external disk as a function of the azimuthal angle @xmath129 . referring to the geometry of the hi density images on the left
, @xmath130 corresponds to the bottom of the disk and increases clockwise in such a way that the disk moves in the direction of @xmath131 .
the gas density and the mean radius are both calculated within sections of a three dimensional annulus with internal and external radius equal to 7 and 15 kpc , respectively .
the initial azimuthal profiles ( not represented in the plots ) are flat since both these quantities have only radial dependence .
as soon as the satellite starts moving through the surrounding medium , the external gas density develops a peak centered on @xmath132 : disk particles localized in regions of maximum ram - pressure get compressed and move on inner orbits , while their circular velocity increases consequently .
after about a quarter of the orbital period the gas has reached its minimum radius and maximum local density .
the gaseous disk becomes strongly asymmetric : compression at the front edge produces a density increase along the left border of the disk , evident in the hi maps even at early times .
the high density region forms a thin ( @xmath133 kpc ) but continuous and well defined arc which has not an equivalent in the stellar distribution . at the perigalacticon this feature extends for almost @xmath134 with a density more than one order of magnitude higher than gas located at smaller radii .
its average thickness ( @xmath135 1.5 kpc ) and velocity dispersion along the line of sight are larger than the average values in the rest of the disk .    in the case of a satellite moving through the hot medium with an inclination angle different from @xmath136 the external pressure directed perpendicularly to the plane of the disk increases as cos@xmath55 while compression at the leading edge is much less pronounced .
[ densitymaps45 ] and [ densitymaps10 ] refer to runs with inclination angles @xmath137 and @xmath115 .
disks are shown face - on .
the increment in density along the leading edge is smaller ( cool45 ) than in the edge - on model and almost absent for @xmath115 ( cool10 ) , while compression perpendicular to the plane of the disk produces local gravitational instabilities in the external gaseous disk ( also mayer , mastropietro & tran in preparation ) .
this effect is more evident in the nearly face - on run cool10 where high density filaments delimitate regions where the local density is almost one order of magnitude lower . despite the absence of a peak in the azimuthal mean density profile ,
the integrated final density of cool10 is comparable to the other runs .
[ dennosf ] represents the azimuthally averaged hi column density profile of the final disk configuration for the three simulations . as a result of the increase in density along the edge of the disk ,
the mean column density shows a secondary peak at large radius .
a limb - brightened density profile has actually been observed by @xcite using the parkes multibeam hi survey of the lmc .
the outer profiles of cool90 and cool45 are very similar , while in run cool10 the increment in density with respect to the original profile of the disk ( long - dashed red curve in the plot ) is located at larger radii ( @xmath138 kpc ) . indeed ,
due to compression of the leading edge the final gas distribution of cool90 and cool45 is asymmetric , the dense front edge being much closer to the center than the opposite border of the disk .
therefore the hi peak in the external disk is located at relatively small radii , while the azimuthally averaged gas distribution is more extended than in the case of the run cool10 since ram pressure elongates the back side of the disk .
again we stress the fact that our gaseous disk is more extended than it would be in a fully self consistent simulation including both gravitational and hydrodynamical forces .
as seen in m05 and @xcite the combined effect of ram - pressure and tidal stripping is quite efficient in stripping gas from the outer satellite s disk , creating the tip of the magellanic stream already at large distances from the mw . the high density feature in cool90
would then form along the border of the disk at smaller radii and would not be easily subjected to further stripping due to the relatively ram - pressure values .
[ cols="<,^,^,^,^,^,^,^",options="header " , ]     [ sfruns ]    the compressive increase in hi density is naturally associated with excess star formation .
sph simulations can not follow the formation of molecular clouds but in first approximation the molecular gas fraction can be related to the density of atomic gas @xcite .    the main parameters of star formation simulations are summarized in table [ sfruns ] .    as we already mentioned in section 3 ,
our standard star formation model ( sf ) does not include the converging flow criterium and is characterized by an efficiency @xmath139 @xcite .
the star formation rate parameter @xmath73 is initially set equal to @xmath140 .
this model was adopted to run wind tube simulations with inclination angles @xmath141 ( sf90 , sf45 , sf10 ) ( table [ sfruns ] ) .
we also investigated different star formation recipes requiring converging flows and assuming different values of @xmath73 and @xmath78 .
we explored star formation rate parameter values in the range from 0.01 to 0.05 , that produce a star formation rate integrated over the entire disk comparable with the @xmath142 yr@xmath86 provided by @xcite .
an efficiency @xmath143 implies that whenever a gas particle satisfies the star formation requirements , it is immediately turned into a single star particle of the same mass @xcite .    in the last six simulations listed in table [ sfruns ] we used our standard star formation model sf to investigate the effects of different orbital parameters and gas halo densities
runs sfv400 are characterized by a maximum ram pressure value corresponding to a perigalactic velocity of 400 km s@xmath86 @xcite . such high velocity disks  when moving face - on through the external hot medium
are strongly affected by local instabilities and star formation is consequently enhanced .
the introduction of an artificial lower limit for the satellite gas temperature ( in runs sf10v400t12000 and sf10v400t15000 ) , higher than the cut - off in the cooling function , has the effect of reducing gravitational instabilities and fragmentation in the disk .
this temperature threshold can be justified in order to crudely model the effect of the uv background and stellar feedback @xcite .    finally , with simulations sf90ld and sf10ld in table [ sfruns ] we also consider the possibility of a galactic hot halo ten times less dense than our standard model .
[ inclstarform ] illustrates the state of the newly formed stellar disk at increasing times along the satellite orbit .
each couple of rows corresponds to one of the first three runs of table [ sfruns ] and is associated with a different inclination angle @xmath55 .
the first row of the pair represents the face - on projection of the disk with the galaxy moving towards the left of the page and the same geometry as in fig .
[ densitymaps ] .
time increases from the left to the right .
each small cross indicates a new star formation event at the time of the snapshot ( within a time interval of 40 myr ) while the circle delimitates the external disk ( @xmath144 kpc ) .
stars form in the central regions as soon as the star formation algorithm is activated , but here the star formation activity of the inner disk is not represented .
the second row represents the total mass @xmath145 of the newly formed stars in the external disk as a function of the disk azimuthal angle @xmath129 .    in the case of a galaxy moving edge - on through the external medium ( sf90 ) ,
stars form at the leading edge of the disk when the ram - pressure becomes larger than @xmath146 dyn @xmath125 , at time @xmath147 gyr .
the location of the star formation events initially corresponds to the hi column density peak observed in cool90 around @xmath148 ( fig .
[ densitymaps ] ) . later on it expands along the entire front edge , creating a thin stellar arc well distinct from the star formation events that characterize the central disk .
as soon as the satellite encounters ram pressure values comparable to those experienced by the lmc at the perigalacticon ( time @xmath149 gyr ) some episodes of star formation occur even on the back side of the disk ( last plot on the top right of fig .
[ inclstarform ] ) , although they are not relevant in terms of new stellar mass formed .
indeed @xmath145 shows a drastic drop at @xmath150
.    runs with inclination @xmath151 are characterized by significant star formation only for values of the external pressure larger than @xmath152 dyn @xmath125 . in the case of the nearly face - on run sf10 , at @xmath153 gyr star formation occurs in the entire external disk .
ram - pressure affects the plane of the disk almost perpendicularly and stars form along the delocalized and filamentary high density structures visible in fig .
[ densitymaps10 ] ( mayer et al . in prep . ) .
contrary to what has been found by @xcite who focused on higher ram - pressure values ( @xmath154 cm@xmath5 and @xmath155 km s@xmath86 ) typical of the outskirts of galaxy clusters , the newly formed stars are all located in the plane of the satellite s disk ( with the exception of the high velocity face - on run sf10v400 where about 10@xmath156 of the stars forms behind the disk ) .
the star formation events appear to be distributed nearly homogeneously along the azimuthal profile of the external disk , although a small peak in @xmath145 is observable near @xmath123 .
in fact , the orientation of the disk with respect to its orbital motion is not exactly face - on .
the case of the intermediate run sf45 is more complex . in a first phase , for low ram pressure values ,
star formation is produced by compression at the leading edge and a thin star formation front
although not so well defined as in the case of a pure edge - on model
 appears on the east side of the disk .
as soon as the external pressure reaches a critical level compression directed perpendicularly to the disk becomes the dominant mechanism driving star formation .
+ converting sfrs to h@xmath0 luminosities according to @xcite :    @xmath157    where sfr is the star formation rate averaged over the last 40 myr ( nearly two times the stellar age of 30 doradus ) , we obtain the h-@xmath0 maps illustrated in fig .
[ lalpha ] .
high emission regions are mainly concentrated in the external disk ( with the exception of run sfld90 where the ram - pressure exerted by the low density halo is not enough to induce star formation at the edge of the disk ) .
the continuous stellar arc forming along the leading side of the disk in edge - on runs breaks up into several distinct and very luminous h-@xmath0 regions that more closely resemble the star - forming complexes observed on the eastern border of the lmc . the inclusion of stellar and supernovae heating  which has been neglected in the present simulations
 could prevent further star - formation around highly emitting regions and consequently produce more compact and isolated star - forming complexes . nevertheless , modeling single star - formation complexes  whose linear extension is smaller than our softening length  is beyond the scope of this paper .
the present inclination of the lmc s disk with respect to the orbital motion is about @xmath158 ( according to the convention adopted in this paper ) .
since the satellite is currently near a perigalactic passage we expect the h-@xmath0 map at the leading border of the disk to be something in between pure edge - on runs and the run with inclination of @xmath159 . on the other side ,
it is very likely that the disk inclination during the phase of approach to the pericenter was different .
indeed in @xcite we have simulated the lmc s orbit according to the new proper motion measurements of @xcite and found that the cloud enters the mw halo face - on and moves almost face - on during most of the last 1 gyr .
it turns nearly edge - on only at the perigalacticon .
this would have a remarkable effect on the star - formation history of the external disk during the last 1 gyr and some impact also on the h-@xmath0 maps .
indeed , although the h-@xmath0 emission would be mostly concentrated on the eastern side of the disk due to the very recent edge - on motion , we expect to see some luminous clumps forming a patchy distribution on the entire disk , due to gravitational instabilities and subsequent star - formation induced by a nearly face - on compression of the gaseous disk before 30 myr ago .
the high velocity edge - on run sf90v400 presents a more elongated and thinner stellar arc along the leading border , with a geometry similar to that obtained increasing the star formation rate parameter to 0.05 ( sf90c0.05 ) .
the h-@xmath0 map of sf90v400 ( third panel of the third raw ) shows two distinct luminosity peaks .
one is located at the south - east region of the disk , roughly corresponding to the position of 30doradus and the two compact emission regions n159 and n160 .
= 9truecm    fig .
[ dmdt ] illustrates how the star formation rate of the external disk changes with time in the different models .
the three black curves refer to the standard star formation runs sf , characterized by low orbital velocities .
the edge - on disk starts forming stars earlier , but for large ram pressure values the star formation rate of sf10 and sf45 grows faster . at the perigalacticon passage sf45 has indeed a higher star formation rate than sf90 . in the case of an isolated lmc model star formation
is almost absent for @xmath160 kpc .
the remaining curves in fig . [ dmdt ] refer to different star formation recipes ( rows 4 - 9 of table [ sfruns ] ) .
the location of the star formation events in the external regions of the disk does not change significantly choosing different parameters in the star formation algorithm .
the star formation rate parameter @xmath73 determines the amount of new stars forming but does not affect the minimum threshold in ram - pressure neither the evolution of the star formation rate . in particular , in the case of edge - on runs after an initial steep increment the curve seems to converge to a constant value for increasing external pressures .
the consequences of an increased star formation efficiency ( @xmath161 ) are almost negligible ( but not in h@xmath0 maps where only the very recent star formation rate is taken in account : compare the first panels of the first and second row ) while including the convergency requirement ( sfconv ) has nearly the same effect than reducing the star formation rate parameter of a factor two .
= 9truecm    for convenience , the star formation rate of the last six runs of table [ sfruns ] is plotted separately ( fig .
[ dmdtalt ] ) . the high velocity edge - on model sf90v400
is characterized by a steeper increment in star formation at earlier times but later on the curve flattens and the star formation rate at the perigalacticon is similar to that of the low velocity case sf90 .
differences in h-@xmath0 maps are produced by a difference of @xmath162 m@xmath163 myr@xmath86 about time = 1 gyr .
on the contrary , the star formation generated by compression perpendicular to the disk increases with increasing ram - pressure values ( it shows a decrement only towards the end of the simulation ) and in the case of the high velocity run sf10v400 reaches a peak @xmath164 times higher than in sf10 .
the star formation rate of sf10v400 is strongly affected by the introduction of an artificial temperature threshold . with a temperature floor of 15000 k
we nearly suppress star formation in the external disk .
however , a threshold of 12000 k is already very high for a lmc model and more typical of luminous disk galaxies like the mw .
a temperature floor lower than 10000 k would not make sense since below this temperature the cooling function adopted in the present paper drops rapidly to zero . if , as pointed out by @xcite , the satellite is moving almost face - on until it gets very close to perigalacticon , we would not expect to see star formation before 0.6 gyr independently of the orbital velocity . differences between high and low velocity runs should be marginal also near the pericenter since star formation in edge - on runs  assuming our standard prescriptions for star formation
 seems to saturate around @xmath165 m@xmath163 myr@xmath86 .
finally , a hot halo ten times less dense than the one assumed in our favorite model would reduce drastically the star formation in a face - on lmc and suppress completely the star formation on the leading edge of the disk .
= 9truecm    the star formation rate of the entire disk is illustrated in fig .
[ dmdttot ] .
this plot is only indicative since we neglect the presence of the bar and its influence on the star formation history of the satellite .
the star formation rate of a lmc model evolved in isolation is plotted for comparison .
the largest contribution to sfr is given by star formation events in the central region of the disk .
curves peak between 15 and 35 m@xmath163 myr@xmath86 , comparable with observations of the magellanic clouds photometric survey ( harris et al . in prep ) .
the total star formation rate is clearly not affected by ram - pressure before 0.4 gyr .
indeed , after an initial sharp increment curves are rather flat despite orientation and intensity of the external pressure . an initial burst in star formation is obtained only by increasing the star formation efficiency parameter to 0.05 . for time@xmath166 0.5 gyr
the largest deviations from an average star formation rate of @xmath167 m@xmath163 myr@xmath86 are produced in the edge - on low density run sf90ld , where dm / dt drops to 10 m@xmath163 myr@xmath86 at time = 1 gyr , and in the high velocity face - on model sf10v400 whose star formation rate increases up to @xmath168 m@xmath163 myr@xmath86 toward the end of the simulation .
differences among the other models are of the order of few m@xmath163 myr@xmath86 and vary with time so that it would be quite difficult to use the star formation history of the entire disk to test the cloud s orbital parameters and the hot halo density .
[ agephi ] represents the mean stellar age of the external disk versus the azimuthal angle @xmath129 for the same selected runs of fig .
[ dmdttot ] . the maximum increment in age in a clockwise direction
is associated with edge - on runs , where stars forming at the leading edge move , in time , away to the side , due to the clockwise rotation of the disk .
the youngest stars are located at @xmath169 .
clearly the gradient in age is much weaker in models with @xmath151 . sf10v400  with a mean stellar age of @xmath170 myr
forms stars earlier with respect to the other nearly face - on runs , while the difference between sf90v400 and the corresponding low velocity run sf90 is about 5 myr at the leading edge and not significant in the rest of the external disk .    for the same runs we also plotted the final radial gas density profile ( fig .
[ densf ] ) . in most of the cases
a secondary peak is still present .
= 8truecm    = 8truecm
we have performed high resolution `` wind tunnel '' simulations to study the effects of ram - pressure by a tenuous galactic hot halo on the hi morphology of the lmc s disk , its recent star formation history , and location of the youngest star forming regions .
we did not focus on the mass loss produced by ram - pressure stripping since this would also be affected by tidal interactions . for the same reason our galaxies do not form any bar and we actually start with a galaxy model very stable against bar formation so that the pure effects of external pressure are more clearly visible .
our lmc is a multi - component system with a spherical nfw halo , an exponential stellar disk and a gaseous disk that extends up to 8 times the stellar disk scale length . in each simulation
the external flux of hot particles increases with time as the satellite approaches the perigalacticon in such a way that the pressure experienced by the disk is consistent with the lmc s orbital velocity and an average hot halo density of @xmath93 cm@xmath5 within 150 kpc from the galactic center .
low velocity runs are characterized by a `` classic '' pericentric velocity of 300 km s@xmath86 @xcite while high velocity runs have velocities compatible with the new proper motion measurements of @xcite .
we expect the angle between the lmc s disk and its proper motion to vary significantly during the last billion years of orbital interaction .
we have defined the inclination angle @xmath55 as the angle between the angular momentum vector of the disk and the flux of hot particles in the wind tube , so that the observed lmc s disk would have @xmath171 .
@xcite performed self consistent nbody / sph simulations of the interacting system mw / lmc adopting orbital constraints from the last lmc s proper motion measurements and found that the cloud enters the mw halo face - on ( @xmath172 ) , moving nearly face - on for most of the last billion years .
it turns edge - on only about 30 myr ago .
this means that the lmc is moving nearly edge - on close to the perigalactic passage , corresponding to the maximum ram - pressure values , consistently with the actual disk inclination measured by @xcite .
we have performed several simulations varying the inclination angle of the disk , the star formation recipe and the intensity of the external pressure .
we have shown that :    * the compression of the leading border of an edge - on lmc disk can account for the high density hi region observed at the south east . in our simulations
this high density feature is well defined ( with a mean density one order of magnitude higher than the surrounding gas ) and localized within 1.5 kpc from the border of the disk .
its average thickness and velocity dispersion along the line of sight are larger than the average values in the rest of the disk . in cool90
it extends for almost @xmath134 and could also explain the origin of the spiral arm e described by @xcite , which does not have an equivalent in the stellar disk . * compression directed perpendicularly to the disk ( in runs with @xmath151 ) produces local instabilities in the gas distribution and a clumpy structure characterized by voids and high density filaments similar to those observed by the parkes multibeam hi survey ( see fig . 3 of staveley - smith et al .
if the satellite was moving nearly face - on in the past  and according to @xcite this is likely to happened during most of the lmc / mw orbital history  ram - pressure could be responsible for the general mottled appearance of the hi disk .
* as a result of the increase in density along the edge of the disk the mean hi column density shows a secondary peak at large radius , in agreement with observations . * the compression of the satellite s igm is naturally associated with induced star formation activity .
we focussed on the external regions of the disk since the central parts of the real lmc would be dominated by the bar .
edge - on disks start forming stars earlier , but for large ram - pressure values the star formation rate of runs with @xmath151 grows much faster .
the high velocity edge - on model sf90v400 is characterized by a steeper increment in star formation at earlier times but later on the curve flattens and the star formation rate at the perigalacticon is similar to that of the low velocity case sf90 .
on the other hand , the star formation generated by a compression perpendicular to the lmc s disk increases with increasing ram - pressure values and in the case of the high velocity run sf10v400 reaches a peak @xmath164 times higher than in sf10 .
if the satellite is moving almost face - on until it gets very close to the perigalacticon we would expect not to see star formation before 0.6 gyr independently of the orbital velocity .
differences between high and low velocity runs should be marginal also near the pericenter since star formation in edge - on runs  assuming our standard prescriptions  seems to saturate around @xmath165 m@xmath163 myr@xmath86 . * in edge - on models the star formation of the external disk is characterized by a thin stellar arc along the leading border , well distinct from the star formation events in the central disk . if sfr is converted in h@xmath0 luminosities , this arc breaks in several distinct and very luminous h@xmath0 regions that more closely resemble the star forming complexes observed on the eastern border of the disk .
although the h@xmath0 emission is mostly concentrated on the eastern side as a consequence of the very recent edge - on motion , we expect to see some luminous clumps forming a patchy distribution on the entire disk , due to gravitational instabilities and subsequent star formation induced by a nearly face - on compression of the disk before 30 myr ago . as observed by @xcite stellar complexes on the leading edge show a progression in age in the clockwise direction , but a face - on compression in the recent past of the lmc would circumscribe this trend to the youngest stellar regions , with age @xmath173 30 - 40 myr .
we would like to thank m - r .
cioni , t. kaufmann and n. kallivayalil for useful discussions .
the numerical simulations were performed on the zbox1 supercomputer at the university of zurich and on the local sgi - altix 3700 bx2 ( partly funded by the cluster of excellence `` origin and structure of the universe '' ) .
this work was partly supported by the dfg sonderforschungsbereich 375 `` astro - teilchenphysik '' .


code coverage is a common metric in software and hardware testing that measures the degree to which an implementation has been tested with respect to some criterion . in its simplest form ,
one starts with a model of the program , and a partition of the behaviors of the model into _ coverage goals _
test _ is a sequence of inputs that determines a behavior of the program .
the aim of testing is to explore as many coverage goals as possible , ideally as quickly as possible . in this paper , we give complexity results for several coverage problems .
the problems are very basic in nature : they consist in deciding whether a certain level of coverage can be attained in a given system .
it is thus somewhat surprising that the problems have not been considered previously in the literature .
finite - state directed graphs have been used as program models for test generation of reactive systems for a long time ( see @xcite for surveys ) .
a coverage goal is a partition of the states of the graph , and a test is a sequence of labels that determine a path in the graph .
the maximal coverage test generation problem is to hit as many partitions as possible using a minimum number of tests . in the special case
the partitions coincide with the states , the maximal coverage problem reduces to the chinese postman problem for which there are efficient ( polynomial time ) algorithms @xcite . in this paper
, we show that the maximal coverage problem becomes np - complete for graphs with general partitions .
we also distinguish between _ system complexity _ ( the complexity of the problem in terms of the size of the graph ) and the _ coverage complexity _ ( the complexity of the problem in terms of the number of coverage goals ) .
then , the problem is nlogspace in the size of the graph ( but that algorithm uses space polynomial in the number of propositions ) .
we consider the special case where the graph has a special `` reset '' action that takes it back to the initial state .
this corresponds in a testing setting to the case where the system can be re - initialized before running a test . in this case , the maximal coverage problem remains polynomial , even with general partitions .
directed graphs form a convenient representation for deterministic systems , in which all the choices are under the control of the tester .
testing of non - deterministic systems in which certain actions are controllable ( under the control of the tester ) and other actions are uncontrollable lead to _ game graphs _ @xcite .
a game graph is a directed labeled graph where the nodes are partitioned into tester - nodes and system - nodes , and while the tester can choose the next input at a tester node , the system non - deterministically chooses the next state at a system node .
then , the test generation problem is to generate a test set that achieves maximal coverage no matter how the system moves . for general game graphs ,
we show the complexity of the maximal coverage problem is pspace - complete .
however , there is an algorithm that runs in time linear in the size of the game graph but exponential in the number of coverage goals .
again , the re - initializability assumption reduces the complexity of coverage : in case there is a re - initialization strategy of the tester from any system state , the maximal coverage problem for games is co - np - complete .
dually , we show that the problem of whether it is possible to win a safety game while visiting fewer than a specified number of partitions is np - complete .    finally , we consider the coverage problem in bounded time , consisting in checking whether a specified number of partitions can be visited in a pre - established number of steps .
we show that the problem is np - complete for graphs , and is pspace - complete for game graphs .
optimization problems arising out of test generation have been studied before in the context of both graphs and games @xcite . however , to the best of our knowledge , the complexities of the coverage problems studied here have escaped attention so far .
while we develop our theory for the finite - state , discrete case , we can derive similar results for more general models , such as those incorporating incomplete information ( the tester can only observe part of the system state ) or timing . for timed systems modeled as timed automata ,
the maximal coverage problem is pspace - complete . for timed games as well as for ( finite state ) game graphs with incomplete information ,
the maximal coverage problem becomes exptime - complete .
in this section we define _ labeled graphs _ and _ labeled games _ , and then define the two decision problems of coverage , namely , _ maximal coverage _ problem and _ coverage with bounded time _ problem .
we start with definition of graphs and games .
a _ labeled graph _ @xmath0 consists of the following component :    1 .   a finite directed graph with vertex set @xmath1 and edge set @xmath2 ; 2 .   the initial vertex @xmath3
a finite set of atomic propositions @xmath4 ; 4 .   a labeling function @xmath5 that assigns to each vertex @xmath6 the set @xmath7 of atomic propositions true at @xmath8 .
for technical convenience we will assume that for all vertices @xmath9 , there exists @xmath10 such that @xmath11 , i.e. , each vertex has at least one out - going edge .    * paths in graphs and reachability .
* given a labeled graph @xmath12 , a _ path _
@xmath13 in @xmath12 is a infinite sequence of vertices @xmath14 starting from the initial vertex @xmath3 ( i.e. , @xmath15 ) such that for all @xmath16 we have @xmath17 .
a vertex @xmath18 is reachable from @xmath3 if there is a path @xmath19 in @xmath12 and @xmath20 such that the vertex @xmath21 in @xmath13 is the vertex @xmath18 .
a _ labeled game graph _ @xmath22 consists of the components of a labeled graph along with a partition of the finite vertex set @xmath1 into @xmath23
. the vertices in @xmath24 are player  1 vertices where player  1 chooses outgoing edges , and analogously , the vertices in @xmath25 are player  2 vertices where player  2 chooses outgoing edges . again for technical convenience we will assume that for all vertices @xmath9 , there exists @xmath10 such that @xmath11 , i.e. , each vertex has at least one out - going edge .    *
plays and strategies in games . *
a _ play _ in a game graph is a path in the underlying graph of the game .
a strategy for a player in a game is a recipe to specify how to extend the prefix of a play .
formally , a strategy @xmath26 for player  1 is a function @xmath27 that takes a finite sequence of vertices @xmath28 ending in a player  1 vertex @xmath6 , where @xmath29 and @xmath30 , representing the history of the play so far , and specifies the next vertex @xmath31 choosing an out - going edge ( i.e. , @xmath32 .
a strategy @xmath33 is defined analogously .
we denote by @xmath34 and @xmath35 the set of all strategies for player  1 and player  2 , respectively .
given strategies @xmath36 and @xmath37 for player  1 and player  2 , there is a unique play ( or a path ) @xmath38 such that ( a )  @xmath15 ; ( b )  for all @xmath39 , if @xmath40 , then @xmath41 ; and if @xmath42 , then @xmath43 .
* controllably recurrent graphs and games . * along with general labeled graphs and games , we will also consider graphs and games that are _ controllably recurrent_. a labeled graph @xmath12 is _ controllably recurrent _ if for every vertex @xmath18 that is reachable from @xmath3 , there is a path starting from @xmath18 that reaches @xmath3 .
a labeled game graph @xmath12 is _ controllably recurrent _ if for every vertex @xmath18 that is reachable from @xmath3 in the underlying graph , there is a strategy @xmath26 for player  1 such that against all player  2 strategies @xmath37 , the path starting from @xmath18 given the strategies @xmath26 and @xmath37 reaches @xmath3 .
controllable recurrence models the natural requirement that systems under test are _ re - initializable _ , that is , from any reachable state of the system , there is always a way to bring the system back to its initial state no matter how the system behaves .    * the maximal coverage problem . *
the _ maximal coverage problem _ asks whether at least @xmath44 different propositions can be visited .
we now define the problem formally for graphs and games . given a path @xmath45 , let @xmath46 be the set of propositions that appear in @xmath13 .
given a labeled graph @xmath12 and @xmath47 , the maximal coverage problem asks whether there is path @xmath13 such that @xmath48 . given a labeled game graph @xmath12 and @xmath47
, the maximal coverage problem asks whether player  1 can ensure that at least @xmath44 propositions are visited , i.e. , whether @xmath49 such that for all player  2 strategies @xmath50 we have @xmath51 .    the maximal _ state coverage
_ problem is the special case of the maximal coverage problem where @xmath52 and for each @xmath53 we have @xmath54 .
that is , each state has its own label , and there are @xmath55 singleton partitions .
* the coverage with bounded time problem . * the _ coverage with bounded time problem _ asks whether at least @xmath44 different propositions can be visited within @xmath56-steps .
we now define the problem formally for graphs and games . given a path @xmath57 and @xmath58 , we denote by @xmath59 the prefix of the path of length @xmath60 , i.e. , @xmath61 . given a path @xmath57 and @xmath58 , we denote by @xmath62 . given a labeled graph @xmath12 and @xmath47 and @xmath63 , the coverage with bounded time problem asks whether there is path @xmath13 such that @xmath64 . given a labeled game graph @xmath12 and @xmath47
, the maximal coverage problem asks whether player  1 can ensure that at least @xmath44 propositions are visited within @xmath56-steps , i.e. , whether @xmath49 such that for all player  2 strategies @xmath50 we have @xmath65 .
* system - tester game . * a _ system _ @xmath66 consists of the following components :    * a finite set @xmath67 of states with the starting state @xmath68 . * a finite alphabet @xmath69 of input letters .
* a transition relation @xmath70 . * a finite set of atomic propositions @xmath4 and a labeling function @xmath5 that assigns to each state @xmath71 the set of atomic propositions true at @xmath71 .
we consider _ total _ systems such that for all @xmath72 and @xmath73 , there exists @xmath74 such that @xmath75 .
a system is _ deterministic _ if for all @xmath72 and @xmath76 , there exists exactly one @xmath77 such that @xmath75 . the tester selects an input letter at every stage and the system resolves the non - determinism in transition to choose the successor state .
the goal of the tester is to visit as many different propositions as possible .
the interaction between the system and the tester can be reduced to a labeled game graph @xmath78 as follows :    * _ vertices and partition .
_ @xmath79 ; @xmath80 and @xmath81 ; and @xmath82 .
* _ edges .
_ @xmath83 . * _ labeling . _ @xmath84 and @xmath85 .
the coverage question for game between tester and system can be answered by answering the question in the game graph .
also observe that if the system is deterministic , then for all player  2 vertices in the game graph , there is exactly one out - going edge , and hence the game can be reduced to a labeled graph . in this paper
we will present all the results for the labeled graph and game model .
all the upper bounds we provide follow also for the game between tester and system .
all the lower bounds we present can also be easily adapted to the model of the game between system and tester .
in this section we study the complexity of the maximal coverage problem . in subsection
[ subsec : maxexpl : graphs ] we study the complexity for graphs , and in subsection  [ subsec : maxexpl : games ] we study the complexity for game graphs .
we first show that the maximal coverage problem for labeled graphs is np - complete .
[ thrm : maxexpl - graphs ] the maximal coverage problem for labeled graphs is np - complete .
the proof consists of two parts .
we present them below .    1 .   _ in np . _
the maximal coverage problem is in np can be proved as follows . given a labeled game graph @xmath12 , let @xmath86 .
we show first that if there is a path @xmath13 in @xmath12 such that @xmath48 , then there is a path @xmath87 in @xmath12 such that @xmath88 .
if @xmath13 visits at least @xmath44 propositions , and there is a cycle in @xmath13 that does not visit a new proposition that is already visited in the prefix , then the cycle segment can be removed from @xmath13 and still the resulting path visits @xmath44 propositions .
hence if the answer to the maximal coverage problem is `` yes '' , then there is a path @xmath87 of length at most @xmath89 that is a witness to the `` yes '' answer . since @xmath90
, it follows that the problem is in np .
_ np - hardness .
_ now we show that the maximal coverage problem is np - hard , and we present a reduction from the sat - problem .
consider a sat formula @xmath91 , and let @xmath92 be the set of variables and @xmath93 be the set of clauses .
for a variable @xmath94 , let 1 .
@xmath95 be the set of indices of the set of clauses @xmath96 that is satisfied if @xmath97 is set to be true ; and 2 .
@xmath98 be the set of indices of the set of clauses @xmath96 that is satisfied if @xmath97 is set to be false .
+ without loss of generality , we assume that @xmath99 and @xmath100 are non - empty for all @xmath101 ( this is because , for example , if @xmath102 , then we can set @xmath97 to be true and reduce the problem where the variable @xmath97 is not present ) . for a finite set @xmath103 of natural numbers , let @xmath104 and @xmath105 denote the maximum and minimum number of @xmath106 , respectively , for an element @xmath107 that is not the maximal element let @xmath108 denote the next highest element to @xmath109 that belongs to @xmath106 ; i.e. , ( a )  @xmath110 ; ( b )  @xmath111 ; and ( c )  if @xmath112 and @xmath113 , then @xmath114 .
we construct a labeled game graph @xmath115 as follows .
we first present an intuitive description : there are states labeled @xmath116 , and all of them are labeled by a single proposition .
the state @xmath117 is an absorbing state ( state with a self - loop only ) , and all other @xmath118 state has two successors .
the starting is @xmath119 . in every state @xmath118
given the right choice we visit in a line a set of states that are labeled by clauses that are true if @xmath118 is true ; and given the left choice we visit in a line a set of states that are labeled by clauses that are true if @xmath118 is false ; and then we move to state @xmath120 .
we now formally describe every component of the labeled graph @xmath121 .
the set of vertices is @xmath122 there is a vertex for every variable , and a vertex @xmath117 .
there is a vertex @xmath123 iff @xmath124 , and there is a vertex @xmath125 iff @xmath126 , 2 .
the set of edges is @xmath127 we now explain the role if each set of edges .
the first edge is the self - loop at @xmath117 .
the second set of edges specifies that from @xmath128 the next vertex is @xmath129 and similarly , from @xmath130 the next vertex is again @xmath129 .
the third set of edges specifies that from @xmath97 there are two successors that are @xmath123 and @xmath131 where @xmath132 and @xmath133 .
the final sets of edges specifies ( a ) to move in a line from @xmath134 to visit the clauses that are satisfied by setting @xmath97 as true , and ( b ) to move in a line from @xmath135 to visit the clauses that are satisfied by setting @xmath97 as false .
fig  [ fig : np - hardness ] gives a pictorial view of the reduction .
3 .   the initial
vertex is @xmath136 .
@xmath137 , i.e. , there is a proposition @xmath138 for each clause @xmath138 and there is a proposition @xmath139 for all variables ; 5 .
@xmath140 ; i.e. , every variable state is labeled by the proposition @xmath139 ; and we have @xmath141 and @xmath142 , i.e. , each state @xmath123 and @xmath125 is labeled by the corresponding clause that is indexes .
+ the number of states in @xmath115 is @xmath143 , and the reduction is polynomial in @xmath91 . in this graph
the maximal number of propositions visited is exactly equal to the maximal number of satisfiable clauses plus  1 ( since along with the propositions for clauses the proposition @xmath139 for all variables is always visited ) .
the proof of the above claim is as follows .
given a path @xmath13 in @xmath144 we construct an assignment @xmath145 for the variables as follows : if the choice at a vertex @xmath97 is @xmath134 , then we set @xmath97 as true in @xmath145 , else we set @xmath97 as false . hence if a path in @xmath115 visits a set @xmath146 of @xmath147 propositions , then the assignment @xmath145 satisfies @xmath148 clauses ( namely , @xmath149 )
conversely , given an assignment @xmath145 of the variables , we construct a path @xmath150 in @xmath115 as follows : if @xmath97 is true in the assignment @xmath145 , then the path @xmath150 chooses @xmath134 at @xmath97 , otherwise , it chooses @xmath151 at @xmath97 .
if @xmath145 satisfies a set @xmath67 of @xmath148 clauses , then @xmath150 visits @xmath152 propositions ( namely , the set @xmath153 of propositions ) .
hence @xmath91 is satisfiable iff the answer to the maximal coverage problem with input @xmath115 and @xmath154 is true .
the desired result follows .
* hardness of approximation .
* we note that from the proof theorem  [ thrm : maxexpl - graphs ] it follows that the max - sat problem ( i.e. , computing the maximal number of clauses satisfiable for a sat formula ) can be reduced to the problem of computing the exact number for the maximal coverage problem . from hardness of approximation of the max - sat problem  @xcite
, it follows that the maximal coverage problem for labeled graphs is hard to approximate .
the maximal coverage problem for labeled graphs that are controllably recurrent can be decided in ptime .    to solve the maximal coverage problem for labeled graphs that are controllably recurrent ,
we compute the maximal strongly connected component @xmath155 that @xmath3 belongs to . since the graph is controllably recurrent , all states that are reachable from @xmath3 belong to @xmath155 .
hence the answer to the maximal coverage problem is `` yes '' iff @xmath156 .
the result follows .
[ thrm : maxexpl - games ] the maximal coverage problem for labeled game graphs is pspace - complete .
the proof consists of two parts .
we present them below .    1 .
_ in pspace .
_ we argue that the maximal coverage problem for labeled game graph can be reduced to the coverage in bounded time problem .
the reason is as follows : in a labeled game graph with @xmath157 vertices , if player  1 can visit @xmath44 propositions , then player  1 can visit @xmath44 propositions within at most @xmath158 steps ; because player  1 can always play a strategy from the current position that visits a new proposition that is not visited and never needs to go through a cycle without visiting a new proposition unless the maximal coverage is achieved .
hence it follows that the maximal coverage problems for games reduces to the coverage in bounded time problem .
the pspace inclusion will follow from the result of theorem  [ thrm : pspace - bouexpl ] where we show that the coverage in bounded time problem is in pspace .
2 .   _ pspace - hardness . _
the maximal coverage problem for game graphs is pspace - complete , even if the underlying graph is strongly connected .
the proof is a reduction from qbf ( truth of quantified boolean formulas ) that is known to be pspace - complete  @xcite , and it is a modification of the reduction of theorem  [ thrm : maxexpl - graphs ] .
consider a qbf formula @xmath159 defined on the set @xmath92 of variables , and @xmath93 are the clauses of the formula .
we apply the reduction of theorem  [ thrm : maxexpl - graphs ] with the following modification to obtain the labeled game graph @xmath115 : the partition @xmath160 of @xmath161 is as follows . for a variable @xmath97
if the quantifier before @xmath97 is existential , then @xmath162 ( i.e. , for existentially quantified variable , player  1 chooses the out - going edges denoting whether to set the variable true or false ) ; and for a variable @xmath97 if the quantifier before @xmath97 is universal , then @xmath163 ( i.e. , for universally quantified variable , the opposing player  2 chooses the out - going edges denoting whether to set the variable true or false ) .
the state @xmath117 is a player  2 vertex , and all other vertex has an single out - going edges and can be player  1 state .
given this game graph we have @xmath91 is true iff player  1 can ensure that all the propositions can be visited in @xmath164 .
formally , let @xmath165 and @xmath166 denote the set of all strategies for player  1 and player  2 , respectively , in @xmath115 .
then @xmath91 is true iff latexmath:[$\sup_{{\pi}_1 \in { \pi}_1^\phi }    \inf_{{\pi}_2 \in { \pi}_2^\phi }     observe that since @xmath117 is a player  2 state if we add an edge from @xmath117 to @xmath119 , player  2 will never choose the edge @xmath117 to @xmath119 ( since the objective for player  2 is to minimize the coverage ) .
however , adding the edge from @xmath117 to @xmath119 makes the underlying graph strongly connected ( i.e. , the underlying graph of the game graph becomes controllably recurrent ; but player  1 does not have a strategy to ensure that @xmath119 is reached , so the game is not controllably recurrent ) .
the desired result follows .
* complexity of maximal coverage in controllably recurrent games .
* we will now consider maximal coverage in controllably recurrent games .
our analysis will use fixing memoryless _
randomized _ strategy for player  1 , and fixing a memoryless randomized strategy in labeled game graph we get a labeled markov decision process ( mdp ) .
a labeled mdp consists of the same components as a labeled game graph , and for vertices in @xmath24 ( which are randomized vertices in the mdp ) the successors are chosen uniformly at random ( i.e. , player  1 does not have a proper choice of the successor but chooses all of them uniformly at random ) .
given a labeled game graph @xmath168 we denote by @xmath169 the mdp interpretation of @xmath12 where player  1 vertices chooses all successors uniformly at random .
an _ end component _ in @xmath170 is a set @xmath171 of vertices such that ( i )  @xmath171 is strongly connected and ( ii )  @xmath171 is player  1 _ closed _
, i.e. , for all @xmath172 , for all @xmath173 such that @xmath174 we have @xmath175 ( in other words , for all player  1 vertices , all the out - going edges are contained in @xmath171 ) .
[ lemm - contr - mdp ] let @xmath12 be a labeled game graph and let @xmath169 be the mdp interpretation of @xmath12 . then the following assertions hold .    1 .
let @xmath171 be an end - component in @xmath169 with @xmath176 .
then @xmath177 .
there exists an end - component @xmath178 with @xmath176 such that @xmath179 .
we prove both the claims below .    1 .
if @xmath171 is an end - component in @xmath170 , then consider a memoryless strategy @xmath50 for player  2 , that for all vertices @xmath180 , chooses a successor @xmath175 ( such a successor exists since @xmath171 is strongly connected ) .
since @xmath171 is player  1 closed ( i.e. , for all player  1 out - going edges from @xmath171 , the end - point is in @xmath171 ) , it follows that for all strategies of player  1 , given the strategy @xmath50 for player  2 , the vertices visited in a play is contained in @xmath171 .
the desired result follows .
an optimal strategy @xmath181 for player  1 in @xmath12 is as follows :
1 .   let @xmath182 and @xmath183 ; 2 .   at iteration @xmath184 ,
let @xmath185 represents the set of propositions already visited .
at iteration @xmath184 , player  1 plays a strategy to reach a state in @xmath186 ( if such a strategy exists ) , and then reaches back @xmath3 ( a strategy to reach back @xmath3 always exists since the game is controllably recurrent ) .
3 .   if a new proposition @xmath187 is visited at iteration @xmath184 , then let @xmath188 .
goto step ( b ) for @xmath189 iteration with @xmath190 .
if no state in @xmath186 can be reached , then stop . + the strategy @xmath181 is optimal , and let the above iteration stop with @xmath191 .
let @xmath192 , and let @xmath193 be the set of vertices such that player  1 can reach @xmath139 .
let @xmath194 . then @xmath195 and player  2 can ensure that from @xmath3 the game can be confined to @xmath196 .
hence the following conditions must hold : ( a )  for all @xmath197 , there exists @xmath198 such that @xmath174 ; and ( b )  for all @xmath199 , for all @xmath200 such that @xmath201 we have @xmath202 .
consider the sub - graph @xmath203 where player  2 restricts itself to edges only in @xmath196 .
a bottom maximal strongly connected component @xmath204 in the sub - graph is an end - component in @xmath169 , and we have @xmath205 it follows that @xmath171 is a witness end - component to prove the result .
the desired result follows .
[ thrm : maxexpl - contrgames ] the maximal coverage problem for labeled game graphs that are controallably recurrent is conp - complete .
we prove the following two claims to establish the result .    1 .
_ in conp .
_ the fact that the problem is in conp can be proved using lemma  [ lemm - contr - mdp ] .
given a labeled game graph @xmath12 , if the answer to the maximal coverage problem ( i.e. , whether @xmath206 ) is no , then by lemma  [ lemm - contr - mdp ] , there exists an end - component @xmath171 in @xmath169 such that @xmath207 .
the witness end - component @xmath171 is a polynomial witness and it can be guessed and verified in polynomial time .
the verification that @xmath171 is the correct witness is as follows : we check ( a )  @xmath171 is strongly connected ; ( b )  for all @xmath172 and for all @xmath200 such that @xmath174 we have @xmath208 ; and ( c )  @xmath209 .
hence the result follows .
2 .   _ conp hardness . _
we prove hardness using a reduction from the complement of the _ vertex cover _ problem . given a graph @xmath210 , a set @xmath211 is a _ vertex cover _ if for all edges @xmath212 we have either @xmath213 or @xmath214 . given a graph @xmath215 whether there is a vertex cover @xmath171 of size at most @xmath44 ( i.e. , @xmath216 ) is np - complete  @xcite .
we now present a reduction of the complement of the vertex cover problem to the maximal coverage problem in controallably recurrent games . given a graph @xmath210
we construct a labeled game graph @xmath217 as follows .
let the set @xmath2 of edges be enumerated as @xmath218 , i.e. , there are @xmath219 edges .
the labeled game graph @xmath220 is as follows . 1 .
_ vertex set and partition .
_ the vertex set @xmath221 is as follows : @xmath222 all states in @xmath2 are player  2 states , and the other states are player  1 states , i.e. , @xmath223 , and @xmath224 .
the set @xmath225 of edges are as follows : @xmath226 intuitively , the edges in the game graph are as follows : from the initial vertex @xmath3 , player  1 can choose any of the edges @xmath227 . for a vertex @xmath228 in @xmath221 ,
player  2 can choose between two vertices @xmath229 and @xmath230 ( which will eventually represent the two end - points of the edge @xmath228 ) . from vertices of the form @xmath229 and @xmath230 , for @xmath231 , the next vertex is the initial vertex @xmath3 .
it follows that from all vertex the game always comes back to @xmath3 and hence we have controllably recurrent game .
3 .   _ propositions and labelling . _
@xmath232 , i.e. , there is a proposition for every vertex in @xmath1 and a special proposition @xmath233 .
the vertex @xmath3 and vertices in @xmath2 are labeled by the special proposition @xmath233 , i.e. , @xmath234 ; and for all @xmath227 we have @xmath235 . for a vertex @xmath236 ,
let @xmath237 , where @xmath238 are vertices in @xmath1 , then @xmath239 and @xmath240 .
note that the above proposition assignment ensures that at every vertex that represents an edge , player  2 has the choices of vertices that form the end - points of the edge .
+ the following case analysis completes the proof . * given a vertex cover @xmath171 , consider a player  2 strategy , that at a vertex @xmath241 , choose a successor @xmath236 such that @xmath242 .
the strategy for player  2 ensures that player  1 visits only propositions in @xmath243 , i.e. , at most @xmath244 propositions . *
consider a strategy for player  1 that from @xmath3 visits all states @xmath245 in order .
consider any counter - strategy for player  2 and let @xmath211 be the set of propositions other than @xmath233 visited .
since all the edges are chosen , it follows that @xmath171 is a vertex cover .
hence if all vertex cover in @xmath215 is of size at least @xmath44 , then player  1 can visit at least @xmath154 propositions .
+ hence there is a vertex cover in @xmath215 of size at most @xmath44 if and only if the answer to the maximal coverage problem in @xmath12 with @xmath154 is no .
it follows that the maximal coverage problem in controllably recurrent games is conp - hard .
the desired result follows .
* complexity of minimal safety games .
* as a corollary of the proof of theorem  [ thrm : maxexpl - contrgames ] we obtain a complexity result about _ minimal safety games_. given a labeled game graph @xmath12 and @xmath44 , the minimal safety game problem asks , whether there exists a set @xmath171 such that a player can confine the game in @xmath171 and @xmath171 contains at most @xmath44 propositions .
an easy consequence of the hardness proof of theorem  [ thrm : maxexpl - contrgames ] is minimal safety games are np - hard , and also it is easy to argue that minimal safety games are in np .
hence we obtain that the minimal safety game problem is np - complete .
in this section we study the complexity of the coverage in bounded time problem . in subsection
[ subsec : bouexpl : graphs ] we study the complexity for graphs , and in subsection  [ subsec : bouexpl : games ] we study the complexity for game graphs .
[ thrm : bouexpl - graphs ] the coverage in bounded time problem for both labeled graphs and controllably recurrent labeled graphs is np - complete .
we prove the completeness result in two parts below .    1 .
_ given a labeled graph with @xmath157 vertices , if there a path @xmath13 such that @xmath64 , then there is path @xmath87 such that @xmath246 .
the above claim follows since any cycle that does not visit any new proposition can be omitted .
hence a path of length @xmath247 can be guessed and it can be then checked in polynomial time if the path of length @xmath248 visits at least @xmath44 propositions . 2 .   _ in np - hard .
_ we reduce the _ hamiltonian - path ( ham - path ) _  @xcite problem to the coverage in bounded time problem for labeled graphs . given a directed graph @xmath210 and an initial vertex @xmath6 , we consider the labeled graph @xmath12 with the directed graph @xmath215 , with @xmath6 as the initial state and @xmath249 and @xmath250 for all @xmath251 , i.e.
, each vertex is labeled with an unique proposition .
the answer to the coverage is bounded time with @xmath252 and @xmath253 , for @xmath86 is `` yes '' iff there is a ham - path in @xmath215 starting from @xmath6 .
the desired result follows .    *
complexity in size of the graph .
* we now argue that the maximal coverage and the coverage in bounded time problem on labeled graphs can be solved in non - deterministic log - space in the size of the graph , and polynomial space in the size of the atomic propositions . given a labeled graph @xmath12 , with @xmath157 vertices
, we argued in theorem  [ thrm : maxexpl - graphs ] that if @xmath44 propositions can be visited , then there is a path of length at most @xmath158 , that visits @xmath44 propositions .
the path of length @xmath158 , can be visited , storing the current vertex , and guessing the next vertex , can checking the set of propositions already visited .
hence this can be achieved in non - deterministic log - space in the size of the graph , and polynomial space in the size of the proposition set .
a similar argument holds for the coverage in bounded time problem .
this gives us the following result .
[ thrm : nlogmaxexpl ] given a labeled graph @xmath254 , the maximal coverage problem and the coverage in bounded time problem can be decided in nlogspace in @xmath255 , and in pspace in @xmath256 .
[ thrm : pspace - bouexpl ] the coverage in bounded time problem for labeled game graphs is pspace - complete .    we prove the following two cases to prove the result .    1 .   _ pspace - hardness_. it follows from the proof of theorem  [ thrm : maxexpl - games ] that the maximal coverage problem for labeled game graphs reduces to the coverage in bounded time problem for labeled game graphs .
since the maximal coverage problem for labeled game graphs is pspace - hard ( theorem  [ thrm : maxexpl - games ] ) , the result follows .
2 .   _ in pspace . _
we say that an _ exploration game tree _ for a labeled game graph is a rooted , labeled tree which represents an unfolding of the graph .
every node @xmath257 of the tree is labeled with a pair @xmath258 , where @xmath6 is a node of the game graph , and @xmath259 is the set of propositions that have been visited in a branch leading from the root of the tree to @xmath257 .
the root of the tree is labeled with @xmath260 .
a tree with label @xmath258 has one descendant for each @xmath261 with @xmath262 ; the label of the descendant is @xmath263 . + in order to check if @xmath44 different propositions can be visited within @xmath56-steps , the pspace algorithm traverses the game tree in depth first order .
each branch is explored up to one of the two following conditions is met : ( i )  depth @xmath56 is reached , or ( ii )  a node is reached , which has the same label as an ancestor in the tree .
the bottom nodes , where conditions ( i ) or ( ii ) are met , are thus the leaves of the tree .
in the course of the traversal , the algorithm computes in bottom - up fashion the _ value _ of the tree nodes .
the value of a leaf node labeled @xmath258 is @xmath264 . for player-1 nodes ,
the value is the maximum of the values of the successors ; for player-2 nodes , the value is the minimum of the value of the successors .
thus , the value of a tree node @xmath257 represents the minimum number of propositions that player  1 can ensure are visited , in the course of a play of the game that has followed a path from the root of the tree to @xmath257 , and that can last at most @xmath56 steps .
the algorithm returns yes if the value at the root is at least @xmath44 , and no otherwise .
+ to obtain the pspace bound , notice that if a node with label @xmath265 is an ancestor of a node with label @xmath266 in the tree , we have @xmath267 : thus , along a branch , the set of propositions appearing in the labels increases monotonically . between two increases , there can be at most @xmath268 nodes , due to the termination condition ( ii ) .
thus , each branch needs to be traversed at most to depth @xmath269 , and the process requires only polynomial space .
the result follows .
[ thrm : bouexpl - contrgames ] the coverage in bounded time problem for labeled game graphs that are controllably recurrent is both np - hard and conp - hard , and can be decided in pspace .
it follows from the ( pspace - inclusion ) argument of theorem  [ thrm : maxexpl - games ] that the maximal coverage problem for labeled game graphs that are controllably recurrent can be reduced to the coverage in bounded time problem for labeled game graphs that are controllably recurrent .
hence the conp - hardness follows from theorem  [ thrm : maxexpl - contrgames ] , and the np - hardness follows from hardness in labeled graphs that are controllably recurrent ( theorem  [ thrm : bouexpl - graphs ] ) .
the pspace - inclusion follows from the general case of labeled game graphs ( theorem  [ thrm : pspace - bouexpl ] ) .
theorem  [ thrm : bouexpl - contrgames ] shows that for controllably recurrent game graphs , the coverage in bounded time problem is both np - hard and conp - hard , and can be decided in pspace .
a tight complexity bound remains an open problem .
* complexity in the size of the game . *
the maximal coverage problem can alternately be solved in time linear in the size of the game graph and exponential in the number of propositions . given a game graph @xmath270 , construct the game graph @xmath271 where @xmath272 , @xmath273 iff @xmath274 and @xmath275 , @xmath276 for @xmath277 , @xmath278 , and @xmath279 .
clearly , the size of the game graph @xmath203 is linear in @xmath215 and exponential in @xmath4 . now consider a reachability game on @xmath203 with the goal @xmath280 .
player-1 wins this game iff the maximal coverage problem is true for @xmath215 and @xmath44 propositions . since
a reachability game can be solved in time linear in the game , the result follows .
a similar construction , where we additionally track the length of the game so far , shows that the maximal coverage problem with bounded time can be solved in time linear in the size of the game graph and exponential in the number of propositions .
[ thrm : lintimemaxexpl ] given a labeled game graph @xmath168 the maximal coverage and the coverage in bounded time problem can be solved in linear - time in @xmath281 and in exponential time in @xmath256 .
somewhat surprisingly , despite the central importance of graph coverage in system verification , several basic complexity questions have remained open .
the basic setting of this paper on graphs and games can be extended in various directions , enabling the modeling of other system features .
we mention two such directions .
so far , we have assumed that at each step , the tester has complete information about the state of the system under test . in practice , this may not be true , and the tester might be able to observe only a part of the state .
this leads to graphs and games of _ imperfect information _ @xcite .
the maximal coverage and the coverage in bounded time problem for games of imperfect information can be solved in exptime .
the algorithm first constructs a perfect - information game graph by subset construction  @xcite , and then run the algorithm of theorem  [ thrm : lintimemaxexpl ] , that is linear in the size game graph and exponential in the number of propositions , on the perfect - information game graph .
thus , the complexity of this algorithm is exptime .
the reachability problem for imperfect - information games is already exptime - hard @xcite , hence we obtain an optimal exptime - complete complexity .
second , while we have studied the problem in the discrete , finite - state setting , similar questions can be studied for timed systems modeled as timed automata @xcite or timed game graphs @xcite .
such problems would arise in the testing of real - time systems .
we omit the standard definitions of timed automata and timed games .
the maximal coverage problem for timed automata ( respectively , timed games ) takes as input a timed automaton @xmath282 ( respectively , a timed game @xmath282 ) , with the locations labeled by a set @xmath4 of propositions , and a number @xmath44 , and asks whether @xmath44 different propositions can be visited .
an algorithm for the maximal coverage problem for timed automata constructs the region graph of the automaton  @xcite and runs the algorithm of theorem  [ thrm : nlogmaxexpl ] on the labeled region graph .
this gives us a pspace algorithm .
since the reachability problem for timed automata is pspace - hard , we obtain a pspace - complete complexity .
similar result holds for the coverage in bounded time problem for timed automata .
similarly , the maximal coverage and coverage in bounded time problem for timed games can be solved in exponential time by running the algorithm of theorem  [ thrm : lintimemaxexpl ] on the region game graph .
this gives an exponential time algorithm .
again , since game reachability on timed games is exptime - hard @xcite , we obtain that maximal coverage and coverage in bounded time in timed games is exptime - complete .
e.  brinksma and j.  tretmans .
testing transition systems : an annotated bibliography . in _
movep 00 : modeling and verification of parallel processes _ ,
volume 2067 of _ lecture notes in computer science _ , pages 187195 .
springer , 2000 .
o.  maler , a.  pnueli , and j.  sifakis . on the synthesis of discrete controllers for timed systems . in _ proc . of 12th annual symp . on theor .
asp . of comp .
_ , volume 900 of _ lect .
notes in comp .
_ springer - verlag , 1995 .