%I#52 2023年10月8日09:56:49
%S 0,1,3,2,5,4,8,11,6,13,12,7,22,16,17,21,9,14,10,27,18,15,19,28,31,43,
%电话:34,38,23,39,25,36,41,20,55,63,42,30,24,33,26,32,65,66,51,59,29,56,35,
%U 62,85,81,37,49,46,48,74,78,88,48,44,75,40,47,97,76,93,79,92,54,58100,61107,52
%N词典学上最早的非负整数无限序列,{a(N)}表示N>=0,这样方程y=a(N。
%我们所说的“一般位置”是指没有两条线平行或重合,也没有三条线同时存在。
%H Luc Rousseau,插图,前五行</a>
%e第0行是y=0*x+0,因为除了默认的约束(词典学上最早的)之外,还没有其他约束适用,所以a(0)=0;
%e线#1不能是y=0*x+1,因为重复使用先前定义的线的斜率意味着平行;但它可以是y=1*x+1,所以a(1)=1;
%e第2行不能是y=2*x+2,因为第0、#1和#2行将在(-1,0)处并发;但第2行可以是y=3*x+2,所以a(2)=3。。。
%o(PARI)
%o g(n,an,p)=我的(b=1);对于(i=1,#p,如果(p[i][2]==an*p[i][1]+n,b=0;中断);b条
%o c(n,s,p)=我的(an=0);while((setsearch(s,an,0)=0)||!g(n,an,p),an++);一个
%o w(n,a,an,p)=对于(m=0,#a-1,my(am=a[1+m]);列表输入(p,[(m-n)/(an-am),(am*n-an*m)/(am-an)]);对
%o list_a(N)=我的(a=列表(),s=列表(”),p=列表(“),N=0,an);而(n<=n,an=c(n,s,p);p=w(n,a,an,p);列表(a,an);listinsert(s,an,setsearch(s,an,1));n++);车辆(a)
%o列表a(75)
%Y参考A000217(创建线#n后要避免的交点数量)。
%K nonn公司
%0、3
%吕克·卢梭,2023年9月22日