链接:https://codeforces.ml/contest/1436/problem/C
题意:输入n,x,pos.(n表示全排列长度,x在数组pos位置上(从0开始))
问多少个排列满足经过以下代码能输出为真:
题解:根据上面的代码模拟就好,mid<=pos,则l=mid+1,if(mid<pos) v[mid]=-1(即应该小于x)。否则v[mid]=1.
然后就是一个组合数,统计v为1的数为z,-1为f。c=n-z-f-1;
我的错误是,没有考虑z>n-x||f>x-1。这样的话应该输出为0。
我考虑到了,但是却没有去实施(高四之后习惯性放弃emmmm)。须知少时凌云志,曾许人间第一流。
代码:
#include <bits/stdc++.h>
#define ll long long
#define pi acos(-1)
#define pb push_back
#define mst(a, i) memset(a, i, sizeof(a))
#define pll pair<ll, ll>
#define fi first
#define se second
#define mp(x,y) make_pair(x,y)
#define dbg(x) cout << #x << "===" << x << endl
using namespace std;
template<class T>void read(T &x){T res=0,f=1;char c=getchar();while(!isdigit(c)){if(c=='-')f=-1;c=getchar();}while(isdigit(c)){res=(res<<3)+(res<<1)+c-'0';c=getchar();}x=res*f;}
void print(ll x){if(x<0){putchar('-');x=-x;}if(x>9)print(x/10);putchar(x%10+'0');}
const ll maxn = 1e6 + 10;
const ll mod = 1e9 + 7;
ll v[maxn],n,x,pos;
ll inv(ll a){return a==1?1:(ll)(mod-mod/a)*inv(mod%a)%mod;}
ll C(ll n,ll m){
if(m>n) return 0;
if(n==0) return 1;
ll a=1,b=1,c=1;
for(ll i=1;i<=n;i++) a=a*i%mod;
for(ll i=1;i<=m;i++) b=b*i%mod;
for(ll i=1;i<=n-m;i++) c=c*i%mod;
return inv(b)%mod*a%mod*inv(c)%mod;
}
//ll gcd(ll a,ll b){return (b==0)?a:gcd(b,a%b);}
ll qpow(ll a,ll p,ll mod){ll ans=1;a=a%mod;while(p){if(p&1)ans=(ans*a)%mod;p>>=1;a=(a*a)%mod;}return ans;}
int main() {
ll _s = 1;
//read(_s);
for (ll _=1;_<=_s;_++) {
read(n),read(x),read(pos);
ll l=0,r=n;
while(l<r){
ll mid=(l+r)/2;//?所到之处
if(mid<=pos){
if(mid<pos) v[mid]=-1;
l=mid+1;
}
else{
//if(mid!=n)
v[mid]=1;
r=mid;
}
}//?所以上面应该没啥问题
//dbg(l-1);
ll z=0,f=0,c=0;
for(ll i=0;i<n;i++){
if(v[i]==1) z++;
if(v[i]==-1) f++;
}
c=n-z-f-1;//?x-1个选z个然后排
//dbg(z),dbg(f),dbg(c);
ll ans=1;
ans=ans*C(x-1,f)%mod;
ans=ans*C(n-x,z)%mod;
for(ll i=1;i<=z;i++) ans=ans*i%mod;
for(ll i=1;i<=f;i++) ans=ans*i%mod;
for(ll i=1;i<=c;i++) ans=ans*i%mod;
print(ans);
}
return 0;
}
/*
input:::
output:::
*/