bzoj 3622 已经没有什么好害怕的了
2015-11-16 10:39
288 查看
我好害怕这种题 φ(≧ω≦*)?两个数组排序后,处理出next[i]表示满足tang[i]>yao[j]的最大的j。f[i][j]前i种糖果,有j个糖比药多,不考虑剩余情况的方案数//也就是剩余n-j个糖果的放法算一种,最后乘上阶乘。f[i][j]=f[i-1][j]+f[i-1][j-1]*max(next[i]-j+1,0) 开始忘记max了。因为f[i][j]不一定都是合法状态,所以要判一下。套容斥:ans[j]=f
[j]*jie[n-j]-sigma(ans[k]*c[k][j])总结一下,容斥原理对于等于k个的可以使用上面的公式,而不是+1-1.spring我算明白了?吐槽一下,这题样例也太强了吧,随便怎么写都能过。
[j]*jie[n-j]-sigma(ans[k]*c[k][j])总结一下,容斥原理对于等于k个的可以使用上面的公式,而不是+1-1.spring我算明白了?吐槽一下,这题样例也太强了吧,随便怎么写都能过。
#include<cstdio> #include<cstdlib> #include<cstring> #include<algorithm> #include<cmath> #define md 1000000009 #define ll long long #define inf (int) 1e9 #define eps 1e-8 #define N 2010 using namespace std; int tang ,yao ,next ; ll f ,c ,jie ; int main() { int n,k; scanf("%d%d",&n,&k); if ((n+k)&1) { printf("0\n"); return 0; } k=(n+k)>>1; for (int i=1;i<=n;i++) scanf("%d",&tang[i]); for (int i=1;i<=n;i++) scanf("%d",&yao[i]); sort(tang+1,tang+n+1); sort(yao+1,yao+n+1); int w=1; for (int i=1;i<=n;i++) { while (tang[i]>yao[w]&&w<=n) w++; next[i]=w-1; } //for (int i=1;i<=n;i++) printf("%d ",next[i]); printf("\n%d\n",k); f[0][0]=1; for (int i=1;i<=n;i++) for (int j=0;j<=n;j++) f[i][j]=(f[i-1][j]+f[i-1][j-1]*max(next[i]-j+1,0))%md; c[0][0]=1; for (int i=1;i<=n;i++) for (int j=0;j<=i;j++) c[i][j]=(c[i-1][j]+c[i-1][j-1])%md; jie[0]=1; for (int i=1;i<=n;i++) jie[i]=(jie[i-1]*i)%md; for (int i=n-1;i>=k;i--) { f [i]=(f [i]*jie[n-i])%md; ll sum=0; for (int j=i+1;j<=n;j++) sum=(sum+f [j]*c[j][i])%md; f [i]=(f [i]-sum+md)%md; } printf("%lld\n",f [k]); return 0; }
相关文章推荐
- 几种线程池的实现算法分析【转载】
- OC中的特有语法
- 搭建jenkins环境(linux操作系统)
- Codeforces Round #331 (Div. 2)
- java代码封装
- 客户端登录也危险?--防范SQL注入漏洞
- ft5x06获取PROJECT_ID/ FW_VER /VENDOR_ID/ glass_vendor
- 看了zzb的安排
- bzoj 3198 spring
- bzoj 1042 硬币购物
- bzoj 2563 阿狸和桃子的游戏
- bzoj 1853 幸运数字
- bzoj 3620 似乎在梦中见过的样子
- bzoj 3916 friends | 字符串hash
- bzoj 1047 理想的正方形
- bzoj 3918 postman
- bzoj 3098 hash killer
- bzoj 1875 HH去散步
- bzoj 3679 数字之积
- bzoj 屯题计划 完结