[codeforces 804F. Fake bullions]
题目大意:
传送门。
给一个n个点的有向完全图(即任意两点有且仅有一条有向边)。
每一个点上有$S_i$个人,开始时其中有些人有真金块,有些人没有金块。当时刻$i$时,若$u$到$v$有边,若$u$中第$i%S_u$个人有金块(无论真假),且$v$中第$i%S_v$个人没有金块,则会给$v$中第$i%S_v$个人一个假金块。
假设这样传递了无数次。(即不会再满足上面的条件时)
最后,拥有真金块的人一定可以把金子买出去,而拥有假的人有一半的概率买出去,每买出去一个会给自己的点贡献1的价值。
问题是:对于每种所有人卖出金块的结果,卖出最多的前$a$个点中会等概率挑出$b$个点,每个不同的$b$的方案(即存在一个点x在方案一中出现而在方案二中没有出现,称为两个不同的方案)会给答案贡献1,求答案对$1e9+7$的余。
题解:
二柱子找到的神题……(它和dp有什么关系……
显然我们如果能求出最后每个点有多少假金子这就是一道组合裸体。
先不考虑真假。
part1:
考虑$u$对$v$的贡献。若u中第x个人有金子,则v中的y均会拥有金块,显然
。
即所有满足$y \equiv x \mod (s_v,s_u)$的y所有金子。这个推导很简单,把x移到左边然后贝祖定理就完了。
即相当与把x映射到$(s_v,s_u)$的剩余系中。我们将其称为$u$对$v$的贡献,用$F(gcd,x)$来表示$u$的元素$x$在$gcd$剩余系下的结果。
当$v$又向$w$连边时,我们会发现x对w的贡献是每一个y直接贡献的,即 。
类比$x$与$y$的关系我们知道$z\equiv y\mod (s_v,s_w)$。所有有贡献 。
再来思考$x$与$z$的关系。有一个显然的结果是$z\equiv x \mod \gcd (s_u,s_v,s_w)$。
这个也很好证明:
我们知道:设变量$a,b,c,d$满足$c|a,d|b$则会有$(c,d)|a+b$。
我们令$c=(s_u,s_v),d=(s_v,s_w),a=x-y,d=y-z$。自然就可以得到这个关系了。
同时我们再考虑是否所有满足这个关系的$z$一定会得到金块。
即:通过$z$与$x$的关系和$z$与$y$的关系,构造出一个$y$使得能满足$y$与$x$的关系:
这样我们就证明了逆定理的正确性。
这用归纳法我们可以证明:设$u$到$v$的一条路径上,设路径上所有点的最大公约数为$gcd$,$u$上$x$对$v$的贡献为$F(gcd,x)$。
part2:
可是当出现环呢?
我们发现上面的结论时刻成立的,两点间的贡献只和路径$gcd$有关。环内两点的路径可以遍历环内每一个元素,我们会发现最后环里任意两点$gcd$为所有元素最大公约数以后不再变化。那么我们就知道对于环内每个点上的x会对环上所有点产生$F(gcd,x)$的贡献。这样我们用状态$H$来表示在$gcd$的剩余系下都有那些位置被标记过,即$H=\cup_x F(gcd,x)$。
part3:
由于这是一个有向完全图,当我们缩点以后,会发现这还是一个有向完全图,不过这次是有向完全DAG。
我们已经计算出了缩点后每个环点的状态。
考虑在这个有向完全DAG的性质。
1.点的出度是递减的。
2.点的入度是递增的。
3.根据拓扑序遍历我们会得到一条链,且链上点标号递减。
以下为缩点后点数为5的例子(由于linux不太会用且校园网有限制只给出有向边,请读者手画一下……):
5->4,5->3,5->2,5->1;4->3,4->2,4->1;3->2,3->1;2->1
我们考虑此时环对环的贡献可以转化为$H$间的贡献,其贡献方式与点完全相同。由于缩点以后图的特殊性质,我们如果根据两点之间有边计算贡献无疑会T得飞起。以5对3的贡献为例,一种是5先给4贡献,4在计算此时$H_4$对3的贡献时再把5的贡献传给3,另一种是5对3直接计算贡献。显然我们发现这两种最后的结果都是一样的。那么前者无疑是跟为优秀的,我们按照拓扑后得到的链来传递贡献的话,那么就可以在只枚举一次每个环的计算出所有的贡献。
代码:
#include "bits/stdc++.h" using namespace std; inline int read(){
int s=,k=;char ch=getchar();
while (ch<''|ch>'') ch=='-'?k=-:,ch=getchar();
while (ch>&ch<='') s=s*+(ch^),ch=getchar();
return s*k;
} const int N=5e3+,M=2e6+,mod=1e9+;
typedef long long ll; struct edges {
int v;edges *last;
}edge[N*N],*head[N<<],*ecnt=edge; inline void push(int u,int v){
*ecnt=(edges){v,head[u]},head[u]=ecnt++;
} inline int powmod(int a,int b){
int ret=;
while (b) {
if(b&) ret=(ll)ret*a%mod;
a=(ll)a*a%mod;
b>>=;
}return ret;
} inline int gcd(int x,int y){
return y?gcd(y,x%y):x;
} int bcc_cnt,bccno[N],stk[N],dfn[N],low[N],idx,top,bgcd[N],s[N],n,a,b,instk[N]; inline void tarjan(int x) {
low[x]=dfn[x]=++idx;
stk[++top]=x;
instk[x]=true;
for (edges *i=head[x];i;i=i->last)
if(!dfn[i->v])
tarjan(i->v),low[x]=min(low[x],low[i->v]);
else if(instk[i->v]) low[x]=min(low[x],dfn[i->v]);
if(dfn[x]==low[x]) {
bcc_cnt++;int t;
do instk[t=stk[top--]]=false,bccno[t]=bcc_cnt,bgcd[bcc_cnt]=gcd(bgcd[bcc_cnt],s[t]);while (t!=x);
}
} char me[N],gold[M<<];
int fr[N<<],in[N],now[N],num[N],mx[N],mn[N]; int inv[N],fac[N]; #define C(i,j) ((ll)fac[(i)]*inv[(j)]%mod*inv[(i)-(j)]%mod) int main(){
//freopen("1.in","r",stdin);
//freopen("war.out","w",stdout);
n=read(),a=read(),b=read();
for (int i=;i<=n;++i) {
scanf("%s",me+);
for (int j=;j<=n;++j) if(me[j]^) push(i,j);
}
for (int i=;i<=n;++i) {
s[i]=read();
fr[i]=fr[i-]+s[i-];
scanf("%s",gold+fr[i]);
}
for (int i=;i<=n;++i) if(!dfn[i])
tarjan(i);
bgcd[]=s[n];
for (int i=;i<=bcc_cnt;++i)
fr[i+n]=fr[i-+n]+bgcd[i-];
fr[n+bcc_cnt+]=fr[n+bcc_cnt]+bgcd[bcc_cnt];
for (int i=;i<fr[n+];++i)
gold[i]^=;
for (int i=;i<=n;++i) {
for (edges *j=head[i];j;j=j->last)
if (bccno[i]!=bccno[j->v])
push(bccno[i]+n,bccno[j->v]+n),++in[bccno[j->v]];
for (int j=fr[i];j<fr[i+];++j)
if (gold[j])
gold[fr[bccno[i]+n]+(j-fr[i])%bgcd[bccno[i]]]=true;
}
now[bcc_cnt]=bgcd[bcc_cnt];
for (int i=bcc_cnt;i>;--i) {
now[i-]=gcd(bgcd[i],bgcd[i-]);
for (int j=;j<bgcd[i];++j) if(gold[fr[i+n]+j])
++num[i],gold[fr[i+n-]+j%now[i-]]=true;
}
for (int j=;j<bgcd[];++j) if(gold[fr[+n]+j])
++num[];
for (int i=;i<=n;++i) {
mx[i]=(ll)num[bccno[i]]*s[i]/bgcd[bccno[i]];
for (int j=fr[i];j<fr[i+];++j)
mn[i]+=gold[j];
}
fac[]=;
for (int i=;i<=n;++i) fac[i]=(ll)fac[i-]*i%mod;inv[]=inv[]=;
for (int i=;i<=n;++i) inv[i]=(ll)inv[mod%i]*(mod-mod/i)%mod;
for (int i=;i<=n;++i) inv[i]=(ll)inv[i]*inv[i-]%mod;
int ans=;
for (int i=;i<=n;++i) {
int c1=,c2=;
for (int j=;j<=n;++j) if (mn[j]>mx[i]) ++c1;
if (c1>=a) continue;
for (int j=;j<=n;++j)
if (mn[j]<=mx[i]&&(mx[j]>mx[i]||mx[j]==mx[i]&&j>i)) ++c2;
for (int j=min(b-,min(c2,a--c1));b-j-<=c1&&(~j);--j)
ans=(ans+C(c2,j)*C(c1,b-j-)%mod)%mod;
}
printf("%d\n",ans);
}
[codeforces 804F. Fake bullions]的更多相关文章
- 【Codeforces】【图论】【数量】【哈密顿路径】Fake bullions (CodeForces - 804F)
题意 有n个黑帮(gang),每个黑帮有siz[i]个人,黑帮与黑帮之间有有向边,并形成了一个竞赛完全图(即去除方向后正好为一个无向完全图).在很多年前,有一些人参与了一次大型抢劫,参与抢劫的人都获得 ...
- Solution -「CF 804F」Fake bullions
\(\mathcal{Description}\) Link. 给定 \(n\) 个点的竞赛图,第 \(i\) 个点代表了 \(s_i\) 个人,每个人(0-based)可能有真金条.此后在 ...
- Codeforces 802I Fake News (hard) (SA+单调栈) 或 SAM
原文链接http://www.cnblogs.com/zhouzhendong/p/9026184.html 题目传送门 - Codeforces 802I 题意 求一个串中,所有本质不同子串的出现次 ...
- Codeforces 802I Fake News (hard)
Codeforces 802I 题意:统计所有不同子串出现次数的平方的和. 想法:建一个SAM,$Ans=\sum (step[i]-step[fa[i]])*right[i]^2$ #include ...
- [CF804F]Fake bullions
Solution: 这题可以分为两个部分, 一个部分为处理出每个点最大的金条数与最小的金条数,记为 \([Min_i, Max_i]\) 第二部分为对于 \(n\) 个变量 \(x_i\i ...
- CodeForces 805A Fake NP
直觉. 一段区间中,肯定是$2$的倍数最多,因为区间长度除以$2$得到的数字最大.但只有$1$个数字的时候需要特判. #include <cstdio> #include <cmat ...
- codeforces411div.2
每日CF: 411div2 Solved A CodeForces 805A Fake NP Solved B CodeForces 805B 3-palindrome Solved C CodeFo ...
- Codeforces Round #310 (Div. 2) B. Case of Fake Numbers 水题
B. Case of Fake Numbers Time Limit: 20 Sec Memory Limit: 256 MB 题目连接 http://codeforces.com/contest/5 ...
- 构造 Codeforces Round #310 (Div. 2) B. Case of Fake Numbers
题目传送门 /* 题意:n个数字转盘,刚开始每个转盘指向一个数字(0~n-1,逆时针排序),然后每一次转动,奇数的+1,偶数的-1,问多少次使第i个数字转盘指向i-1 构造:先求出使第1个指向0要多少 ...
随机推荐
- Loader转换器
一.简介 webpack本身只能处理js模块,Loader可以理解为模块和资源的转换器,它本身是一个函数,接受文件作为参数,返回转换的结果.因此,我们就能通过require来加载任何类型的模块和文件. ...
- Integer 和int 比较
在jdk1.5的环境下,有如下4条语句: 1 2 3 4 Integer i01 = 59; int i02 = 59; Integer i03 =Integer.valueOf(59); Integ ...
- Javascript、CSS、HTML面试题
1 JS中的三种弹出式消息提醒(警告窗口.确认窗口.信息输入窗口)的命令是什么? alert confirm prompt 2声明一个已经存在一个CSS有几种方式? 1.导入一个已经存 ...
- Python-Flask框架之——图书管理系统 , 附详解源码和效果图 !
该图书管理系统要实现的功能: 1. 可以通过添加窗口添加书籍或作者, 如果要添加的作者和书籍已存在于书架上, 则给出相应的提示. 2. 如果要添加的作者存在, 而要添加的书籍书架上没有, 则将该书籍添 ...
- MongoDB与CouchDB 全方位对比
http://blog.nosqlfan.com/html/1519.html 本文见于MongoDB官方网站,MongoDB与CouchDB 很相似,他们都是文档型存储,数据存储格式都是JSON型的 ...
- jtds驱动更新对一个老问题的解决
07年年末的一篇blog: 以前网站做初期开发时,有一个问题:hibernate下text大字符串读取时出这个异常:JDBCExceptionReporter - The amount of data ...
- Interesting卡常数
C++ Interesting卡常数 作为一名OIer,在Noip中卡(kǎ 我就爱读kǎ)常数可以说是必备技巧.在此总结一下我所知卡常数的神奇手法: IO优化 fread 和 fwrite ,如果还 ...
- javaXML文件解析之DOM4J实操
既然前面说了DOM4J这里好那里好,大家都是在用这个,那咱就不得不写一个了. XML文件: <?xml version="1.0" encoding="UTF-8& ...
- node八-核心模块、包
学会查API,远比会几个API更重要. 核心模块意义 -如果只是在服务器运行javascript代码,并没有多大意义,因为无法实现任何功能>读写文件.访问网络 -Node的用处在于它本身还提供可 ...
- Android动态字符串拼接----%s
在开发经常遇到字符串中的某一数据或多个数据是动态变化的. 如下图 不要创建3个TextView,暂时不考虑颜色变化的情况,可以用以下做法. <string name="maintain ...