codeforces（D. Rescue Nibel!）（排序，二分）_codeforces 指令集 #pragma gcc target("avx2")-优快云博客

本文链接：https://blog.youkuaiyun.com/qq_43781431/article/details/114016902

这篇博客展示了如何使用C++进行高效的数学计算，包括模运算、快速幂和前缀和，并结合排序算法处理数组。文章通过定义结构体`node`，实现了对数据的排序，并使用`lower_bound`函数进行区间查询。最后，代码演示了如何求解特定数学问题并输出结果。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

原题链接
在这里插入图片描述

//#pragma GCC optimize("Ofast")
//#pragma GCC target("avx,avx2,fma")
//#pragma GCC optimization ("unroll-loops")
#include <bits/stdc++.h>
#define pb push_back
#define ll long long
#define IOS std::ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
inline int read() {
    char c = getchar(); int x = 0, f = 1;
    while (c < '0' || c > '9') {if (c == '-')f = -1; c = getchar();}
    while (c >= '0' && c <= '9') {x = x * 10 + c - '0'; c = getchar();}
    return x * f;
}
using namespace std;

ll n,k,ans;
ll mod=998244353;
struct node
{
    ll x,y;
}e[300005];
ll c[300005],a[300005];
bool cmp(node a,node b)
{
    if(a.x<b.x){return 1;}
    return 0;
}
ll apow(ll a,ll b)
{
    ll s=1;
    while(b)
    {
        if(b&1){s=s*a%mod;}
        a=a*a%mod;
        b>>=1;
    }
    return s;
}
void init()
{
    c[k-1]=1;
    ll val=1,base=1;
    for(ll i=1;i<=k-1;i++){val=val*i%mod;}
    for(ll i=1;i<=k-1;i++){base=base*apow(i,mod-2)%mod;}
    for(ll i=k;i<=n;i++)
    {
        val=val*i%mod;
        val=val*apow(i-(k-1),mod-2)%mod;
        c[i]=val*base%mod;
    }
}
ll query(ll key)
{
    ll pos=lower_bound(a+1,a+1+n,key)-a;
    return n-pos+1;
}
int main()
{
    scanf("%lld %lld",&n,&k);
    for(int i=1;i<=n;i++)
    {
        scanf("%lld %lld",&e[i].x,&e[i].y);
    }
    init();
    sort(e+1,e+1+n,cmp);
    for(int i=1;i<=n;i++){a[i]=e[i].y;}
    sort(a+1,a+1+n);
    for(int i=n;i>=1;i--)
    {
        ll cc=query(e[i].x);
        cc-=(n-i+1);
        if(cc>=(k-1))
        {
            ans=(ans+c[cc])%mod;
        }
    }
    printf("%lld\n",ans);
    return 0;
}