@property
def r1(P):
    """
        Number of real roots with multiplicity
        """
    assert P.isreal()
    ans = 0
    s = P.sturm()
    while s:
        ans += s
        P = P.gcd(P.prime())
        s = P.sturm()
    return ans
