Aviso: Esta ainda não é uma resposta completa. Se os argumentos de plausibilidade o deixarem desconfortável, pare de ler.
Vou considerar uma variante em que queremos multiplicar (x - a_1) ... (x - a_n) pelos números complexos.
O problema é duplo na avaliação de um polinômio em n pontos. Sabemos que isso pode ser feito de maneira inteligente no tempo O (n log n), quando os pontos são a enésima raiz da unidade. Isso tira vantagem essencial das simetrias dos polígonos regulares subjacentes à Transformação Rápida de Fourier. Essa transformação ocorre em duas formas, convencionalmente denominadas decimação no tempo e decimação na frequência. Na raiz dois, eles se apóiam em um par de simetrias de polígonos regulares de lados pares: a simetria entrelaçada (um hexágono regular consiste em dois triângulos equilaterais entrelaçados) e a simetria de desdobramento do ventilador (corta um hexágono regular ao meio e desdobra as peças como ventiladores em triângulos equilaterais).
Nesta perspectiva, parece altamente implausível que um algoritmo O (n log n) exista para um conjunto arbitrário de n pontos sem simetrias especiais. Isso implicaria que não há nada de excepcional em termos de algoritmo nos polígonos regulares em comparação com conjuntos aleatórios de pontos no plano complexo.