【USACO】牛的学术圈 I 一题三解

题目

由于对计算机科学的热爱，以及有朝一日成为 「Bessie 博士」的诱惑，奶牛 Bessie 开始攻读计算机科学博士学位。

经过一段时间的学术研究，她已经发表了 N 篇论文，并且她的第 i 篇论文得到了来自其他研究文献的 ci 次引用。

Bessie 听说学术成就可以用 h 指数来衡量。

h 指数等于使得研究员有至少 h 篇引用次数不少于 h 的论文的最大整数 h。

例如，如果一名研究员有 4 篇论文，引用次数分别为 (1,100,2,3)，则 h 指数为 2，然而若引用次数为 (1,100,3,3) 则 h 指数将会是 3。

为了提升她的 h 指数，Bessie 计划写一篇综述，并引用一些她曾经写过的论文。

由于页数限制，她至多可以在这篇综述中引用 L 篇论文，并且她只能引用每篇她的论文至多一次。

请帮助 Bessie 求出在写完这篇综述后她可以达到的最大 h 指数。

注意 Bessie 的导师可能会告知她纯粹为了提升 h 指数而写综述存在违反学术道德的嫌疑；我们不建议其他学者模仿 Bessie 的行为。

输入格式
输入的第一行包含 N 和 L。

第二行包含 N 个空格分隔的整数 c1,…,cN。

输出格式
输出写完综述后 Bessie 可以达到的最大 h 指数。

数据范围
1≤N≤105,
0≤ci≤105,
0≤L≤105
输入样例1：
4 0
1 100 2 3
输出样例1：
2
样例1解释
Bessie 不能引用任何她曾经写过的论文。上文中提到，(1,100,2,3) 的 h 指数为 2。

输入样例2：
4 1
1 100 2 3
输出样例2：
3
如果 Bessie 引用她的第三篇论文，引用数会变为 (1,100,3,3)。上文中提到，这一引用数的 h 指数为 3。

[解法一： 贪心]

最容易想到的策略，注意证明，这里证明略。

#include <iostream>
#include <cstring>
#include <algorithm>
#include <cmath>

using namespace std;
const int N = 1e5 + 10;
int a[N];

int main()
{
	int n, L;
	cin >> n >> L;
	memset(a, 0, sizeof a);
	for (int i = 1; i <= n; i++) cin >> a[i];
	sort(a + 1, a + n + 1, greater<int>());
	int i = 1, cnt = 0;
	
	while(i <= n && a[i] >= i) i++; //证明，若答案小于i则i一定更优，若大于i + 1则不满足要求，所以i或i + 1一定是答案
	
	if(a[i] < i - 1) cout << i - 1 << endl;
	else{
	    for(int j = 1; j <= i; j++) if(a[j] == i - 1) cnt++;
	    if(L >= cnt) cout << i << endl;
	    else cout << i - 1 << endl;
	}

	return 0;
}

[解法二：二分]

h是二分答案

#include<bits/stdc++.h>
using namespace std;
const int N=1e5+10;
int c[N];
int n,k;
bool check(int x)
{
    int ans=0,bns=0;
    for(int i=1;i<=n;i++)
    {
        if(c[i]>=x)ans++; //记录已经满足条件的论文个数
        else if(x-c[i]==1)bns++; 
        //记录可以通过引用之后，满足条件的论文数

        //只有和x相差1次，才可以通过引用使得总数增加，因为每篇论文只能最多被引用一次
    }
    bns=min(k,bns);  //要和k取min，因为有可能bns>k,这样我k的次数不够改变所有的bns，所以要取两者较小的哪个
    ans+=bns;
    if(ans>=x)return true;
    else return false;
}
int main()
{
    cin>>n>>k;
    for(int i=1;i<=n;i++)scanf("%d",&c[i]);
    //考虑二分
    int l=0,r=N;
    while(l<r)
    {
        int mid=l+r+1>>1;
        if(check(mid))l=mid;
        else r=mid-1;
    }
    cout<<l<<endl;
    return 0;
}

[解法三]双指针

类似于贪心，不是很好理解

#include <iostream>
#include <cstring>
#include <algorithm>

using namespace std;

const int N = 100010;

int n, L;
int q[N];

int main()
{
    scanf("%d%d", &n, &L);
    for (int i = 1; i <= n; i ++ ) scanf("%d", &q[i]);
    sort(q + 1, q + n + 1, greater<int>());

    int res = 0;
    for (int i = 1, j = n; i <= n; i ++ )
    {
        while (j && q[j] < i) j -- ;
        if (q[i] >= i - 1 && i - j <= L)
            res = i;
    }

    printf("%d\n", res);
    return 0;
}