| 
	John Cremona on Wed, 06 Apr 2011 18:33:32 +0200
	 | 
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
	
	| 
        Re: Polynomial divisibility
	 | 
 
- To: Charles Greathouse <charles.greathouse@case.edu>
 
- Subject: Re: Polynomial divisibility
 
- From: John Cremona <john.cremona@gmail.com>
 
- Date: Wed, 6 Apr 2011 09:26:44 -0700
 
- Cc: pari-users@list.cr.yp.to
 
- Delivery-date: Wed, 06 Apr 2011 18:33:32 +0200
 
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed;        d=gmail.com; s=gamma;        h=domainkey-signature:mime-version:in-reply-to:references:date         :message-id:subject:from:to:cc:content-type         :content-transfer-encoding;        bh=tGMooy/giFm9eEbIQAKX62a2Tq9EQMAY76zv/qUTKH0=;        b=pvBluLR04B9YrDGF+6DFGVwatdzNrX1myVmy1ElL+1iwK4T6NDZ0dH1IlV1LXgq5D/         Ov5tsaaH9wv8fTdYD+OVapWyQ4ee1KMSRxhisiKkzLjSdaykReKPEywmJt0OcZW9+3ui         ElUYdZWRTeSTaTD8eRpl9HoN7SWUnbTnyfLtA=
 
- Domainkey-signature: a=rsa-sha1; c=nofws;        d=gmail.com; s=gamma;        h=mime-version:in-reply-to:references:date:message-id:subject:from:to         :cc:content-type:content-transfer-encoding;        b=o62Kw95r6vz5AfiA1YxLwciLghmuD3uFCYHFFAaVcFEbuxP87tgEp1H2JRHAxJs76J         vc9t1ncWWW+8dt15juzUBL4JJIE9Fe9dCDsmS7/SHUAE0sM5wWqqzlpxwD/IMV0kZV7f         9ewDEJ5dIvCHDlJF+0VvtSsqPzf4JbTVzYpv8=
 
- In-reply-to: <BANLkTi=t7FBTR=LHWUAjhTznCkA746O2vQ@mail.gmail.com>
 
- Mailing-list: contact pari-users-help@list.cr.yp.to; run by ezmlm
 
- References: <BANLkTi=t7FBTR=LHWUAjhTznCkA746O2vQ@mail.gmail.com>
 
Do you mean to count the number of roots P has modulo p, not counting
multiplicity?
If so then it is the degree of gcd(P,x^p-x).
John
On Wed, Apr 6, 2011 at 8:56 AM, Charles Greathouse
<charles.greathouse@case.edu> wrote:
> I wanted to know if there's an efficient way to count the number of
> residue classes (mod p) for which the polynomial is divisible by p.
> The straightforward approach
>
> sum(n=1,p,substpol(P,x,Mod(n,p))==0)
>
> is slow.
>
> In my case the polynomial is reducible and of degree 62 with
> 'reasonable' coefficients (wordsize on a 64-bit machine, the largest
> is 44 bits).  I could test the smaller polynomials first but I think
> the overhead would be more expensive than the benefit -- it's rare
> that p will divide any given value of the polynomial.
>
> I'm willing to code in PARI if GP does not suffice.  It's possible
> that there's a different approach that doesn't enumerate cases but I
> don't know of one.
>
> Charles Greathouse
> Analyst/Programmer
> Case Western Reserve University
>