雯舞

Love three things int he world -- Manchery

[计数] 美团 CodeM 复赛 排列

这里写图片描述

#include<cstdio>
#include<cstdlib>
#include<algorithm>
using namespace std;
typedef pair<int,int> abcd;
typedef long long ll;

#define read(x) scanf("%d",&(x))

const int P=1e9+7;

inline ll Pow(ll a,int b){
  ll ret=1;
  for (;b;b>>=1,a=a*a%P)
    if (b&1)
      ret=ret*a%P;
  return ret;
}
inline ll Inv(ll a){
  return Pow(a,P-2);
}

const int N=100005;

int n;
int a[N],b[N],f[N];
int pos[N];

int c[N],maxn;
inline void add(int x,int r){
  for (int i=x;i<=maxn;i+=i&-i)
    (c[i]+=r)%=P;
}
inline int sum(int x){
  int ret=0;
  for (int i=x;i;i-=i&-i)
    (ret+=c[i])%=P;
  return ret;
}
inline int sum(int l,int r){
  if (l>r) return 0;
  return (sum(r)+P-sum(l-1))%P;
}

int main(){
  int x,y;
  freopen("t.in","r",stdin);
  freopen("t.out","w",stdout);
  ll fac=1;
  read(n); for (int i=1;i<=n;i++) read(x),read(y),a[x]=y,fac=fac*i%P,pos[y]=i;
  maxn=n;
  for (int i=n;i;i--)
    b[i]=sum(a[i],n),add(a[i],1),b[i]=Inv(b[i]);
  maxn=n+1; for (int i=1;i<=n;i++) c[i]=0;
  add(0+1,Inv(n));
  for (int i=1;i<=n;i++){
    ll ret=sum(0+1,a[i]-1+1);
    if (!b[i])
      f[pos[a[i]]]=ret*fac%P;
    else
      add(a[i]+1,ret*b[i]%P);
  }
  for (int i=1;i<=n;i++)
    printf("%d\n",f[i]);
  return 0;
}
阅读更多
版权声明:本文为博主原创文章,未经博主允许随意转载。 https://blog.csdn.net/u014609452/article/details/74999915
个人分类: 组合数学 杂题
上一篇[数论] LOJ #510. 「LibreOJ NOI Round #1」北校门外的回忆
下一篇[DP 分块] UOJ #300. 【CTSC2017】吉夫特
想对作者说点什么? 我来说一句

没有更多推荐了,返回首页

关闭
关闭