696 research outputs found
Flocking with discrete symmetry: the 2d Active Ising Model
We study in detail the active Ising model, a stochastic lattice gas where
collective motion emerges from the spontaneous breaking of a discrete symmetry.
On a 2d lattice, active particles undergo a diffusion biased in one of two
possible directions (left and right) and align ferromagnetically their
direction of motion, hence yielding a minimal flocking model with discrete
rotational symmetry. We show that the transition to collective motion amounts
in this model to a bona fide liquid-gas phase transition in the canonical
ensemble. The phase diagram in the density/velocity parameter plane has a
critical point at zero velocity which belongs to the Ising universality class.
In the density/temperature "canonical" ensemble, the usual critical point of
the equilibrium liquid-gas transition is sent to infinite density because the
different symmetries between liquid and gas phases preclude a supercritical
region. We build a continuum theory which reproduces qualitatively the behavior
of the microscopic model. In particular we predict analytically the shapes of
the phase diagrams in the vicinity of the critical points, the binodal and
spinodal densities at coexistence, and the speeds and shapes of the
phase-separated profiles.Comment: 20 pages, 25 figure
Active Brownian Particles and Run-and-Tumble Particles: a Comparative Study
Active Brownian particles (ABPs) and Run-and-Tumble particles (RTPs) both
self-propel at fixed speed along a body-axis that reorients
either through slow angular diffusion (ABPs) or sudden complete randomisation
(RTPs). We compare the physics of these two model systems both at microscopic
and macroscopic scales. Using exact results for their steady-state distribution
in the presence of external potentials, we show that they both admit the same
effective equilibrium regime perturbatively that breaks down for stronger
external potentials, in a model-dependent way. In the presence of collisional
repulsions such particles slow down at high density: their propulsive effort is
unchanged, but their average speed along becomes . A
fruitful avenue is then to construct a mean-field description in which
particles are ghost-like and have no collisions, but swim at a variable speed
that is an explicit function or functional of the density . We give
numerical evidence that the recently shown equivalence of the fluctuating
hydrodynamics of ABPs and RTPs in this case, which we detail here, extends to
microscopic models of ABPs and RTPs interacting with repulsive forces.Comment: 32 pages, 6 figure
Optimal paths on the road network as directed polymers
We analyze the statistics of the shortest and fastest paths on the road
network between randomly sampled end points. To a good approximation, these
optimal paths are found to be directed in that their lengths (at large scales)
are linearly proportional to the absolute distance between them. This motivates
comparisons to universal features of directed polymers in random media. There
are similarities in scalings of fluctuations in length/time and transverse
wanderings, but also important distinctions in the scaling exponents, likely
due to long-range correlations in geographic and man-made features. At short
scales the optimal paths are not directed due to circuitous excursions governed
by a fat-tailed (power-law) probability distribution.Comment: 5 pages, 7 figure
Sparse Suffix and LCP Array: Simple, Direct, Small, and Fast
Sparse suffix sorting is the problem of sorting suffixes of a string
of length . Efficient sparse suffix sorting algorithms have existed for more
than a decade. Despite the multitude of works and their justified claims for
applications in text indexing, the existing algorithms have not been employed
by practitioners. Arguably this is because there are no simple, direct, and
efficient algorithms for sparse suffix array construction. We provide two new
algorithms for constructing the sparse suffix and LCP arrays that are
simultaneously simple, direct, small, and fast. In particular, our algorithms
are: simple in the sense that they can be implemented using only basic data
structures; direct in the sense that the output arrays are not a byproduct of
constructing the sparse suffix tree or an LCE data structure; fast in the sense
that they run in time, in the worst case, or in
time, when the total number of suffixes with an LCP value
greater than is in
, matching the time of the optimal yet much more
complicated algorithms [Gawrychowski and Kociumaka, SODA 2017; Birenzwige et
al., SODA 2020]; and small in the sense that they can be implemented using only
machine words. Our algorithms are simplified, yet non-trivial,
space-efficient adaptations of the Monte Carlo algorithm by I et al. for
constructing the sparse suffix tree in time [STACS
2014]. We also provide proof-of-concept experiments to justify our claims on
simplicity and efficiency.Comment: 16 pages, 1 figur
Range Shortest Unique Substring queries
Let be a string of length n and be the substring of starting at position i and ending at position j. A substring of is a repeat if it occurs more than once in; otherwise, it is a unique substring of. Repeats and unique substrings are of great interest in computational biology and in information retrieval. Given string as input, the Shortest Unique Substring problem is to find a shortest substring of that does not occur elsewhere in. In this paper, we introduce the range variant of this problem, which we call the Range Shortest Unique Substring problem. The task is to construct a data structure over answering the following type of online queries efficiently. Given a range, return a shortest substring of with exactly one occurrence in. We present an -word data structure with query time, where is the word size. Our construction is based on a non-trivial reduction allowing us to apply a recently introduced optimal geometric data structure [Chan et al. ICALP 2018]
Longest common substring made fully dynamic
Given two strings S and T, each of length at most n, the longest common substring (LCS) problem is to find a longest substring common to S and T. This is a classical problem in computer science with an O(n)-time solution. In the fully dynamic setting, edit operations are allowed in either of the two strings, and the problem is to find an LCS after each edit. We present the first solution to this problem requiring sublinear time in n per edit operation. In particular, we show how to find an LCS after each edit operation in Ă(n2/3) time, after Ă(n)-time and space preprocessing. 1 This line of research has been recently initiated in a somewhat restricted dynamic variant by Amir et al. [SPIRE 2017]. More specifically, they presented an Ă(n)-sized data structure that returns an LCS of the two strings after a single edit operation (that is reverted afterwards) in Ă(1) time. At CPM 2018, three papers (Abedin et al., Funakoshi et al., and Urabe et al.) studied analogously restricted dynamic variants of problems on strings. We show that the techniques we develop can be applied to obtain fully dynamic algorithms for all of these variants. The only previously known sublinear-time dynamic algorithms for problems on strings were for maintaining a dynamic collection of strings for comparison queries and for pattern matching, with the most recent advances made by Gawrychowski et al. [SODA 2018] and by Clifford et al. [STACS 2018]. As an intermediate problem we consider computing the solution for a string with a given set of k edits, which leads us, in particular, to answering internal queries on a string. The input to such a query is specified by a substring (or substrings) of a given string. Data structures for answering internal string queries that were proposed by Kociumaka et al. [SODA 2015] and by Gagie et al. [CCCG 2013] are used, along with new ones, based on ingredients such as the suffix tree, heavy-path decomposition, orthogonal range queries, difference covers, and string periodicity
Constructing Antidictionaries of Long Texts in Output-Sensitive Space
A word x that is absent from a word y is called minimal if all its proper factors occur in y. Given a collection of k words y1, ⊠, yk over an alphabet ÎŁ, we are asked to compute the set M{y1,âŠ,yk}â of minimal absent words of length at most â of the collection {y1, ⊠, yk}. The set M{y1,âŠ,yk}â contains all the words x such that x is absent from all the words of the collection while there exist i,j, such that the maximal proper suffix of x is a factor of yi and the maximal proper prefix of x is a factor of yj. In data compression, this corresponds to computing the antidictionary of k documents. In bioinformatics, it corresponds to computing words that are absent from a genome of k chromosomes. Indeed, the set Myâ of minimal absent words of a word y is equal to M{y1,âŠ,yk}â for any decomposition of y into a collection of words y1, ⊠, yk such that there is an overlap of length at least â â 1 between any two consecutive words in the collection. This computation generally requires Ω(n) space for n = |y| using any of the plenty available O(n) -time algorithms. This is because an Ω(n)-sized text index is constructed over y which can be impractical for large n. We do the identical computation incrementally using output-sensitive space. This goal is reasonable when â„M{y1,âŠ,yN}ââ„=o(n), for all N â [1,k], where â„Sâ„ denotes the sum of the lengths of words in set S. For instance, in the human genome, n â 3 Ă 109 but â„M{y1,âŠ,yk}12â„â106. We consider a constant-sized alphabet for stating our results. We show that allMy1â,âŠ,M{y1,âŠ,yk}â can be computed in O(kn+âN=1kâ„M{y1,âŠ,yN}ââ„) total time using O(MaxIn+MaxOut) space, where MaxIn is the length of the longest word in {y1, ⊠, yk} and MaxOut=max{â„M{y1,âŠ,yN}ââ„:Nâ[1,k]}. Proof-of-concept experimental results are also provided confirming our theoretical findings and justifying our contribution
Equality of opportunity and educational achievement in Portugal
Portugal has one of the highest levels of income inequality in Europe, and low wages and unemployment are concentrated among low skill individuals. Education is an important determinant of inequality. However, there are large differences in the educational attainment of different individuals in the population, and the sources of these differences emerge early in the life-cycle when families play a central role in individual development. We estimate that most of the variance of school achievement at age 15 is explained by family characteristics. Observed school inputs explain very little of adolescent performance. Children from highly educated parents benefit of rich cultural environments in the home and become highly educated adults. Education policy needs to be innovative: (1) it needs to explicitly recognize the fundamental long run role of families on child development; (2) it needs to acknowledge the failure of traditional input based policies
- âŠ