TopCoder SRM500 Div1 1000 其他
原文链接https://www.cnblogs.com/zhouzhendong/p/SRM500-1000.html
SRM500 Div1 1000
设 \(v_1,v_2,\cdots ,v_9\) 分别为一个数中 \(1-9\) 的出现次数。
那么可以列出如下方程组:
v_1+2v_2+3v_3+\cdots +9v_9 &=& S\\
v_2+2v_4+3v_8+v_6 &=& p_2\\
v_3+v_6+2v_9 &=& p_3\\
v_5 &=& p_5\\
v_7 &=& p_7
\end{eqnarray*}
\]
有 9 个变量,5 个方程,所以只需要枚举其中 4 个就可以得到所有的。于是问题被转化成了知道 \(v_1\cdots v_9\),求这些数对答案的贡献。先假设所有的 \(t=\sum_{1\leq i \leq 9}v_i\) 个数都互不相同,则任意一个数在任意一个位置都会产生 \((t-1)!\) 次贡献,所以总贡献为 \(a=(\sum_{0\leq i <t} 10 ^i)(\sum_{1\leq i\leq 9}iv_i)((t-1)!)\) 。由于值相同的数是等价的,所以最终得到的贡献为 \(a/(\prod _{1\leq i \leq 9}v_i!)\) 。
代码
static const int N=3005,mod=500500573;
int Pow(int x,int y){
int ans=1;
for (;y;y>>=1,x=1LL*x*x%mod)
if (y&1)
ans=1LL*ans*x%mod;
return ans;
}
int Fac[N],Inv[N],Iv[N],sum[N];
void prework(int n){
for (int i=Fac[0]=1;i<=n;i++)
Fac[i]=1LL*Fac[i-1]*i%mod;
Inv[n]=Pow(Fac[n],mod-2);
for (int i=n;i>=1;i--)
Inv[i-1]=1LL*Inv[i]*i%mod;
for (int i=1;i<=n;i++)
Iv[i]=1LL*Inv[i]*Fac[i-1]%mod;
sum[0]=0;
for (int i=1;i<=n;i++)
sum[i]=(10LL*sum[i-1]+1)%mod;
}
int v1,v2,v3,v4,v5,v6,v7,v8,v9;
int calc(){
static int t;
t=v1+v2+v3+v4+v5+v6+v7+v8+v9;
return 1LL*Fac[t-1]*sum[t]%mod
*(v1*1+v2*2+v3*3+v4*4+v5*5+v6*6+v7*7+v8*8+v9*9)%mod
*Inv[v1]%mod*Inv[v2]%mod*Inv[v3]%mod
*Inv[v4]%mod*Inv[v5]%mod*Inv[v6]%mod
*Inv[v7]%mod*Inv[v8]%mod*Inv[v9]%mod;
}
int getSum(int p2, int p3, int p5, int p7, int S){
prework(3000);
int ans=0;
S-=p5*5+p7*7;
v5=p5,v7=p7;
for (v2=0;v2<=p2;v2++)
for (v4=0;v2+v4*2<=p2;v4++)
for (v8=0;v8*3+v4*2+v2<=p2;v8++){
v6=p2-v2-v4*2-v8*3;
int s2=S-v2*2-v4*4-v8*8-v6*6;
if (s2<0||p3<v6)
continue;
for (v9=(p3-v6)/2;v9>=0;v9--){
v3=p3-v6-v9*2;
v1=s2-v9*9-v3*3;
if (v1<0)
continue;
ans=(ans+calc())%mod;
}
}
return ans;
}
TopCoder SRM500 Div1 1000 其他的更多相关文章
- TopCoder SRM500 Div1 250 其他
原文链接https://www.cnblogs.com/zhouzhendong/p/SRM500-250.html SRM500 Div1 250 题意 (看题用了半个小时--) 有 n 个人(编号 ...
- TopCoder SRM500 Div1 500 分治
原文链接https://www.cnblogs.com/zhouzhendong/p/SRM500-500.html SRM500 Div1 500 没想到 double 的精度居然没有爆-- 考虑以 ...
- TopCoder SRM502 Div1 1000 动态规划
原文链接https://www.cnblogs.com/zhouzhendong/p/SRM502-1000.html SRM502 Div1 1000 题意 从 [0,n-1] 中选择 k 个不同的 ...
- TopCoder 649 div1 & div2
最近一场TC,做得是在是烂,不过最后challenge阶段用一个随机数据cha了一个明显错误的代码,最后免于暴跌rating,还涨了一点.TC题目质量还是很高的,非常锻炼思维,拓展做题的视野,老老实实 ...
- TopCoder SRM502 Div1 500 贪心 01背包
原文链接https://www.cnblogs.com/zhouzhendong/p/SRM502-500.html SRM502 Div1 500 好题. 首先,如果已经确定了解决所有问题的优先级, ...
- topcoder SRM712 Div1 LR
题目: Problem Statement We have a cyclic array A of length n. For each valid i, element i-1 the l ...
- TopCoder 603 div1 & div2
div2 250pts MiddleCode 题意:s串长度为奇数时,将中间字符取掉并添加到t末尾:长度为偶数时,将中间两个较小的字符取掉并添加到末尾. 分析:直接做,学习了一下substr(s, p ...
- TopCoder SRM704 Div1 800 构造
原文链接https://www.cnblogs.com/zhouzhendong/p/SRM704-800.html 题解 考虑构造一个 $n = 20$ 的图. 先把所有 $i$ 都连向 $i-1$ ...
- topcoder SRM642 div1 hard WheelofFortune
题目链接:vjudge 大意:有两个人参加一场游戏,这个游戏在一个编号为\(0\text~n-1\)的轮盘上进行,一开始轮盘上的数字均为0:一共有\(m\)轮,每一轮都有一个操作参数\(s_i\),主 ...
随机推荐
- mysqlfrm初步使用
这个工具也就是读取frm文件生成建表语句默认的模式是再生个实例,使用--basedir选项或指定--server选项来连接到已经安装的实例.这种过程不会改变原始的.frm文件.该模式也需要指定--po ...
- dnsmasq详解&手册
Dnsmasq为小型网络提供网络基础设施:DNS,DHCP,路由器通告和网络引导.它被设计为轻量级且占用空间小,适用于资源受限的路由器和防火墙.它还被广泛用于智能手机和便携式热点的共享,并支持虚拟化框 ...
- FlashBack 闪回
[学习目标] Flashback Database 功能非常类似与RMAN的不完全恢复,它可以把整个数据库回退到 过去的某个时点的状态,这个功能依赖于Flashback log日志.比RMAN 更快速 ...
- Confluence 6 MySQL 输入你的数据库细节
Confluence 的安装向导将会指导你一步一步的在 Confluence 中配置安装 MySQL 数据库. 使用 JDBC 连接(默认) JDBC 是推荐的连接你的 Confluence 到数据库 ...
- Java之动手动脑(三)
日期:2018.10.12 星期五 博客期:017 这次留了两个动手动脑作业!我需要一个一个来说!先说第一个吧! Part 1 :随机生成1000个随机数 代码: //以下为 RandomMaker. ...
- D3.js 添加zoom缩放功能后dblclick双击也会放大的问题
svg.call(zoom).on("dblclick.zoom", null); https://stackoverflow.com/questions/25007466/d3- ...
- iframe内容自适应高度
一直觉得要用JS才能实现iframe高度的自适应,其实CSS也可以,而且实现的更好,只是需要给包裹iframe的DIV设置个高度,然后让irame高度设置成100%就可以自适应了. 完美版Iframe ...
- mysql之视图,触发器,事务等。。。
一.视图 视图是一个虚拟表(非真实存在),其本质是[根据SQL语句获取动态的数据集,并为其命名],用户使用时只需使用[名称]即可获取结果集,可以将该结果集当做表来使用. 使用视图我们可以把查询过程中的 ...
- BrupSuite渗透测试笔记(十一)
一.数据查找和拓展功能的使用 1.BrupSuite高级功能在界面布局上主要分成两个部分,一是菜单栏,另一个是Engagement tools,Brup菜单下包含的数据查找Search 组件状态存储. ...
- eclipse的安装及使用
1.安装 2工作区 3透视图添加透视图 关闭和显示各个子视图 点击视图右上角的关闭按钮可以关闭当前视图 可以选择Window-->Show View菜单项打开各个子视图 4创建项目 选择File ...