考えたこと
- ソートしてよさそう
- lを買う中で最も小さいお好み焼きとすると買える最大のお好み焼きのindexは単調性がある
- [l, r) を尺取みたいな感じでもとめていく
- 区間の要素数からその区間で買うときのパターン数は求まる
- lは必ず買うものとしてr-l-1要素からM-1要素を選ぶのでコンビネーションでok
- 提出したら通った
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef vector<int> VI;
typedef vector<VI> VVI;
typedef vector<ll> VL;
typedef vector<VL> VVL;
typedef pair<int, int> PII;
#define FOR(i, a, n) for (ll i = (ll)a; i < (ll)n; ++i)
#define REP(i, n) FOR(i, 0, n)
#define ALL(x) x.begin(), x.end()
#define IN(a, b, x) (a<=x&&x<b)
#define MP make_pair
#define PB push_back
const int INF = (1LL<<30);
const ll LLINF = (1LL<<60);
const double PI = 3.14159265359;
const double EPS = 1e-12;
const int MOD = 1000000007;
template <typename T> T &chmin(T &a, const T &b) { return a = min(a, b); }
template <typename T> T &chmax(T &a, const T &b) { return a = max(a, b); }
int dx[] = {0, 1, 0, -1}, dy[] = {1, 0, -1, 0};
template<unsigned MOD>
class ModInt {
public:
unsigned x;
ModInt(): x(0) { }
ModInt(signed y) : x(y >= 0 ? y % MOD : MOD - (-y) % MOD) {}
unsigned get() const { return x; }
ModInt inv() const {
ll a = 1, p = x, e = MOD-2;
while(e > 0) {
if(e%2 == 0) {p = (p*p) % MOD; e /= 2;}
else {a = (a*p) % MOD; e--;}
}
a %= MOD;
return ModInt(a);
}
ModInt pow(ll e) {
ll a = 1, p = x;
while(e > 0) {
if(e%2 == 0) {p = (p*p) % MOD; e /= 2;}
else {a = (a*p) % MOD; e--;}
}
a %= MOD;
return ModInt(a);
}
ModInt pow2() {
ll a = 1, p = 2, e = x;
while(e > 0) {
if(e%2 == 0) {p = (p*p) % MOD; e /= 2;}
else {a = (a*p) % MOD; e--;}
}
a %= MOD;
return ModInt(a);
}
bool operator <(ModInt b) { return x < b.x; }
bool operator >(ModInt b) { return x > b.x; }
bool operator<=(ModInt b) { return x <= b.x; }
bool operator>=(ModInt b) { return x >= b.x; }
bool operator!=(ModInt b) { return x != b.x; }
bool operator==(ModInt b) { return x == b.x; }
ModInt operator++() { x++; return *this; }
ModInt operator--() { x--; return *this; }
ModInt &operator+=(ModInt that) {
x = ((ll)x+that.x)%MOD;
return *this;
}
ModInt &operator-=(ModInt that) {
x = ((((ll)x-that.x)%MOD)+MOD)%MOD;
return *this;
}
ModInt &operator*=(ModInt that) {
x = (ll)x * that.x % MOD;
return *this;
}
ModInt &operator/=(ModInt that) {
x = (ll)x * that.inv() % MOD;
return *this;
}
ModInt &operator%=(ModInt that) {
x = (ll)x % that.x;
return *this;
}
ModInt operator+(ModInt that)const{return ModInt(*this) += that;}
ModInt operator-(ModInt that)const{return ModInt(*this) -= that;}
ModInt operator*(ModInt that)const{return ModInt(*this) *= that;}
ModInt operator/(ModInt that)const{return ModInt(*this) /= that;}
ModInt operator%(ModInt that)const{return ModInt(*this) %= that;}
};
typedef ModInt<1000000007> mint;
ostream &operator<<(ostream& os, mint a) { return os << a.x; }
istream &operator>>(istream& is, mint &a) { return is >> a.x; }
ll combi(ll N_, ll C_) {
const int NUM_=400001;
static ll fact[NUM_+1],factr[NUM_+1],inv[NUM_+1];
if (fact[0]==0) {
inv[1]=fact[0]=factr[0]=1;
for (int i=2;i<=NUM_;++i) inv[i] = inv[MOD % i] * (MOD - MOD / i) % MOD;
for (int i=1;i<=NUM_;++i) fact[i]=fact[i-1]*i%MOD, factr[i]=factr[i-1]*inv[i]%MOD;
}
if(C_<0 || C_>N_) return 0;
return factr[C_]*fact[N_]%MOD*factr[N_-C_]%MOD;
}
class OkonomiyakiParty {
public:
int count(vector <int> osize, int M, int K)
{
int n = osize.size();
sort(ALL(osize));
int l = 0, r = 0;
mint ret = 0;
while(l <= r && l < n) {
while(r < n && osize[r] - osize[l] <= K) r++;
ret += combi(r-l-1, M-1);
l++;
}
return ret.x;
}
};