[联合集训6-14] 一道很好玩的几何题

设$k_{i,j}$表示$i,j$连线的斜率，$v_i$表示第$i$个关键点。
给出一个结论：$\forall i,f(v_i)=v_{i+1}$。假设该结论成立，那么就有$k_{v_{i-1},v_i}>k_{v_i,v_i+1}$。
考虑证明这个结论。假设$f(v_i)<v_{i+1}$，如果$k_{v_i,f(v_i)}>k_{v_{i+1},f(v_{i+1})}$，那么直接把$v_{i+1}$接到$f(v_i)$上就减少了长度而且答案不变；如果$k_{v_i,f(v_i)}\le k_{v_{i+1},f(v_{i+1})}$，因为$v_i$看不到$v_{i+1}$所以$k_{v_i,f(v_i)}>k_{f(v_i),v_{i+1}}$，又因为$f(v_i)$看不到$f(v_{i+1})$，所以$k_{f(v_i),v_{i+1}}>k_{v_{i+1},f(v_{i+1})}$，然后就矛盾了。
那么可以把高度先下降$a$，再持平$b-a$的一段看成非负分数$\frac{a}{b}$，那么问题转化成找到一个最长的递增分数序列使得$\sum b\le n$ 。可以发现一定存在一种最优方案使得$a,b\le \sqrt n$，所以只需要把其中的最简分数按分母从小到大排序然后贪心选择即可。（实测对于$n=1000000$，存在方案使得$a,b\le 200$）
代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define pii pair<int,int>
#define fs first
#define sc second
#define MP make_pair
#define N 1000010
using namespace std;
int n,m,h[N],top,tot;
pii a[N],ans[N];
int gcd(int x,int y)
{
    for(int t;y;t=y,y=x%y,x=t);
    return x;
}
bool cmp1(pii p,pii q)
{
    return p.sc<q.sc;
}
bool cmp2(pii p,pii q)
{
    return p.fs*q.sc>p.sc*q.fs;
}
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=min(n,200);i++)
        for(int j=1;j<=i;j++)
            if(gcd(i,j)==1) a[++m]=MP(j,i);
    sort(a+1,a+m+1,cmp1);
    for(int i=1;i<=m;i++)
        if(tot+a[i].sc<n) tot+=a[i].sc,ans[++top]=a[i];
    sort(ans+1,ans+top+1,cmp2); 
    int pnt=n;
    for(int i=1;i<=top;i++)
    {
        for(int j=1;j<=ans[i].sc-ans[i].fs;j++)
            pnt--,h[pnt]=0;
        for(int j=1;j<=ans[i].fs;j++)
            pnt--,h[pnt]=1; 
    }
    for(int i=n;i>=0;i--)
        h[i]+=h[i+1];
    printf("%d\n",top+1);
    for(int i=0;i<=n;i++)
        printf("%d ",h[i]);     
    return 0;
}