试题编号： 202104-4
试题名称： 校门外的树
时间限制： 1.0s
内存限制： 512.0MB 

问题描述：

X 校最近打算美化一下校园环境。前段时间因为修地铁，X 校大门外种的行道树全部都被移走了。现在 X 校打算重新再种一些树，为校园增添一抹绿意。

X 校大门外的道路是东西走向的，我们可以将其看成一条数轴。在这条数轴上有n个障碍物，例如电线杆之类的。虽然障碍物会影响树的生长，但是障碍物不一定能被随便移走，所以 X 校规定在障碍物的位置上不能种树。n个障碍物的坐标都是整数；如果规定向东为正方向，则 n障碍物的坐标按照从西到东的顺序分别为a1,a2,…,an。X 校打算在[a1,an]之间种一些树，使得这些树看起来比较美观。

X 校希望，在一定范围内，树应该是等间隔的。更具体地说，如果把[a1,an)划分成一些区间[ap1,ap2),…,[apm-1,apm), （1=p1＜p2＜···＜pm=n），那么每个区间[api,api+1)内需要至少种一棵树，且该区间内种的树的坐标连同区间端点api,api+1应该构成一个等差数列。不同区间的公差，也就是树的间隔可以不相同。

例如，如果障碍物位于0,2,6这三处，那么我们可以选择在[0,2)和[2,6)分别种树，也可以选择在[0,6)等间隔种树。如果是分别在[0,2)和[2,6)种树，由于每个区间内至少要种一棵树，坐标1上必须种树；而[2,6)上的树可以按照1的间隔种下，也可以按照2的间隔种下。下图表示了这两种美观的种树方案，其中橙色的圆表示障碍物，绿色的圆表示需要在这个位置种树，箭头上的数字表示种下这棵树时对应的间隔为多少。
 

对区间[0,2)和[2,6)分别以1和2的间隔种树是美观的

 

对区间[0,2)和[2,6)分别以1的间隔种树也是美观的

而如果选择在[0,6)区间等间隔种树，我们只能以3的间隔种树，因为无论是选择间隔1或者间隔2，都需要在坐标2上种树，而这个位置已经有障碍物了。下图分别表示了间隔为3,2,1时的种树情况，红色箭头表示不能在这里种树。

对区间[0,6)以3的间隔种树是美观的

对区间[0,6)以2的间隔种树是不美观的

 

对区间[0,6)以1的间隔种树也是不美观的

一般地，给定一个区间[al,ar)，对于树的坐标的集合 T（al,ar)（TZ），归纳定义T在[al,ar)上是美观的：

1.如果T≠Ø，T∩{al,al+1,…,ar}=Ø，并且存在一个公差d≥1，使得T∪{al,ar}中的元素按照从小到大的顺序排序后，可以构成一个公差为d的等差数列（显然，这个等差数列的首项为al，末项为ar），则T在[al,ar)上是美观的；
2.如果T∩{al,al+1,…,ar}=Ø，并且存在一个下标m（l＜m＜r），使得T∩(al,am)在[al,am)是美观的，且T∩(am,ar)在[am,ar)上是美观的，则T在[al,ar)上是美观的。

根据这一定义，空集在任意区间上都不是美观的；另外，如果存在下标i使得ai∈T ，那么T一定不是美观的。

我们称两种种树的方案是本质不同的，当且仅当两种方案中，种树的坐标集合不同。请帮助 X 校对[a1,an)求出所有本质不同的美观的种树方案。当然，由于方案可能很多，你只需要输出总方案数对109+7取模的结果。

输入格式

输入的第一行包含一个正整数n，表示障碍物的数量。

输入的第二行包括n个非负整数a1,…,an，表示每个障碍物的坐标。

保证对i=1,2,…,n-1，ai＜ai+1。

输出格式

输出一个非负整数，表示本质不同的美观的种树方案的数量对109+7取模的结果。

样例输入

3
0 2 6

样例输出

3

样例说明

这组样例即为题面描述中提到的那组。

样例输入

11
0 10 20 30 40 50 60 70 80 90 100

样例输出

256507

样例输入

333
33 44 67 210 528 762 873 984 1234 1466 1739 2859 3421 4061 4598 5172 5201 5220 5261 5322 5389 5559 6670 7070 7898 8079 8129 8192 8616 8641 8806 9559 9585 9750 10263 10627 10674 10692 10903 11649 11885 12179 12307 12743 13173 13352 13389 13496 13611 15292 15321 16018 16327 16415 16959 16972 17499 17617 17786 18476 18966 19239 19498 19875 20312 20392 21603 21620 21730 21967 21972 21999 22015 22590 22775 23709 23839 24165 24408 24595 25160 25479 25812 26482 27328 28101 28297 28305 28342 28557 28986 29110 29401 29765 30292 30493 30739 31027 31201 31218 31414 32089 32759 32770 32777 32815 32877 32890 33297 33457 33603 33757 33866 34498 34525 34659 34679 34861 34870 34997 35311 35846 36411 36457 36738 36902 37940 38228 40156 40320 40705 40737 40803 41066 41443 41460 41954 41968 42040 42062 42099 43281 43320 43527 43537 43587 43729 44750 44822 45655 45769 46109 46525 47060 47128 47999 48635 48887 48981 49366 49424 49524 50546 50580 50689 51332 51861 51943 52097 52702 53009 53067 53397 53526 53901 54280 54399 54801 55535 55592 55740 55843 56110 56428 56552 56682 56848 57179 57688 57797 57847 57959 58330 58831 59553 59699 59884 59939 61233 61636 61732 61908 62145 62549 62649 62740 62912 62971 63053 64312 64322 64412 64816 64845 64873 64923 64976 65023 65166 65496 66065 66491 66803 66941 67081 68331 68336 68360 68476 69179 69719 69758 69948 70072 70544 70598 70990 71014 71454 71687 71743 71958 72282 72384 72456 72985 73327 74325 75046 75097 76647 77062 77088 77431 77553 77673 77753 78217 78518 78564 79565 79588 79686 80275 80939 81052 81348 81386 81440 81589 81610 81793 82408 82801 82836 83239 83466 83610 83867 83943 84441 84467 85248 85305 85554 85565 85758 86251 86603 86743 87323 87565 87824 87833 88265 88309 89178 89509 89618 89699 89708 90331 90359 90878 90902 91449 92284 92374 92549 92609 93609 94345 94934 95140 95475 95733 95985 95995 96270 96641 96807 97003 97632 98160 98677 98853 98943 99037 99055 99075 99185 99395 99592

样例输出

7094396

评测用例规模与约定

对于10%的数据，保证n=2；

对于30%的数据，保证n≤10；

对于60%的数据，保证n≤100，ai≤1000；

对于100%的数据，保证2≤n≤1000，0≤ai≤100,000，且至少存在一种美观的种树方案。
 

解题思路：

dp算法，依次推出每一个点的最大方案数，同时进行标记

首先把范围内所有数的约数进行初始化，这个在后面的遍历中需要用到(类似于打表)

两层for循环进行dp递推

最后结果就是二维最后一个位置的方案数

 

解题代码：

#include<bits/stdc++.h>
#define ll long long
using namespace std;
bool vis[100005];
ll f[1005];
ll a[1005];
vector<int>p[100005];
const ll MOD=1e9+7;
int main()
{
    ios::sync_with_stdio(0);
    cin.tie(0);
    cout.tie(0);

    for(int i=1;i<=100000;i++)
    {
        for(int j=i+i;j<=100000;j+=i)
        p[j].push_back(i);
    }
    int n;cin>>n;
    for(int i=0;i<n;i++) cin>>a[i];

    f[0]=1;
    for(int i=1;i<n;i++)
    {
        memset(vis,false,sizeof vis);
        for(int j=i-1;j>=0;j--)
        {
            ll d=a[i]-a[j];
            ll cnt=0;
            for(int k=0;k<p[d].size();k++)
            {
                if(!vis[p[d][k]])
                {
                    cnt++;
                    vis[p[d][k]]=true;
                }
            }
            vis[d]=true;
            f[i]=(f[i]+(f[j]*cnt)%MOD)%MOD;
        }
    }
    cout<<f[n-1]<<endl;
}