Next: About this document ...
��������� I ���� No.12
������Ⱦ�Υơ���:
���� 12.1
�Ĵ��͡�����
�����������Ƥ���Ȥ��롣 ���ΤȤ�,
��ľ��ʬ��
�ǡ� ��ľ�°���
�������
��
�ΤߤǤ���褦�ʤ�Τ�¸�ߤ��롣
PID �򥤥ǥ���dz�ä��褦�ʴĤˤĤ��ƤϤ��Ǥˤ�äȾܤ����ͻҤ��ΤäƤ��롣
���� 12.2
Ǥ�դβĴ���
������,
��
�Υ��ǥ���Ǥ��롣
�⤷
������ʤ�С�
�������
��
�ΤߤǤ��롣
���� 12.3
����
�� PID �ǡ�
�ʤ顢
�ʤ�ľ��ʬ��ǡ�
���ΤǤ���褦�ʤ�Τ�¸�ߤ��롣
���� 12.4
�Ĵ���
��
��ľ��ʬ�򤵤�Ƥ����Ȥ����,
Ǥ�դ�
-�÷�
��,
-�÷��Ȥ��Ƥ�ľ��ʬ��
�ǡ�
��
�ؤκ��Ѥ�����ľ��ʬ��ˤ�ä��гѤη�
�˽񤱤�褦�ʤ�Τ��ġ�
PID
���Ф���, -�÷����ɤΤ褦�ʹ�¤���Ĥ�����
ñ������������Ǥ��ä���
��Ⱦ�Ǥ��ǰ���ʬ��θ�Ψ���ˤĤ��ƹͤ�������
���ǰ���ʬ�򤹤���ָ���Ū����ˡ�ϡ�
����
�ޤǤ�
�����ǽ缡��äƤ�����ˡ�Ǥ��롣
����������Ǥ�
�����㤹����֤������롣
(��ˡ1)
�ǽ��
�Ĥ������Í�¸Ĥο��ò·«¤ï¿½ï¿½Ö¤ï¿½ï¿½ï¿½ï¿½ï¿½
��Ϳ�����Ƥ���Ȥ��롣
���ΤȤ�,
�Ȥ�����,
�Ȥʤ�
��¸�ߤ��롣
����������¿�༰
����ơ�
�ʤ��������������������ˤ��μ�ˡ��Ŭ�Ѥ����,
����
�ʤ�
����Ȥޤä�,
���ǰ��Ҥ���ޤ��礬���롣����� Pollard ��
ˡ�Ȥ�����
���� 12.1
ˡ�ޤ��Ϥ���¾����ˡ���Ѥ���
���ǰ�������衣
���Ѥ�����ˡ��, �ץ�����ࡢ����Ӽ¹Ի��֤򤽤��뤳�ȡ�
(����)
ˡ�� UBASIC �ǽ񤯤ȼ��Τ褦�ˤʤ롣
(���������ץ�������ʣ���ˤ���Τ��򤱤뤿���
�㴳��Ψ�ΰ����ץ������ˤ��Ƥ��롣
��̣�Τ���ͤϲ��ɤ��뤫����ʬ�ǰ줫��񤤤Ƥߤ뤳�ȡ�)
10 n=607143768775207 'ʬ�򤹤��
30 a=1:b=1 ' a=x_i, b=x_{2i} �ΤĤ��Ǥ��롣
50 while 1 ' �������� wend (150 ��)�ޤ�̵�¥롼��
60 a=fnf(a)@n
70 b=fnf(fnf(b))@n
80 x=gcd(a-b,n)
90 if x<>1 then print x,n :end
150 wend
160 end
1000 def fnf(k) ' fnf() �ؿ������; fn... �ϥ桼������ؿ�
1010 return(k*k+1) ' �������ۤ��δؿ��Ǥ�褤��
ubasic �ʤɤΥץ����������Ȥä����Ȥ��ʤ��ҤȤϡ�
�ֵ��Υۡ���ڡ�������ubasic �����������ɤ�,
ubasic ����ष�����Ȥ��Υե�����˥��Ģ�ʤɤǾåµï¿½×¥ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½
�񤤤ơ�"rho.txt" �ʤɤ���¸��,
���Τ���ubasic ��ư����,
load "rho.txt"
�Ȥ���Ф褤��������
�ֵ��Υڡ����ˤϿ����������ե� mupad �⤪���Ƥ���Τ�,
�����ȤäƤ���������������,mupad �� ifactor �ؿ����Ѥ����
�������˴�ñ������Τǡ�ifactor �����Ϻ���϶ؤ���Ȥ��롣
Yoshifumi Tsuchimoto
2001-01-11