开始给了一个简单的题目,但我还是没有珍惜。

一个简简单单的树形 \(dp\),然而因为取模却不知道该如何比较大小。。

其实可以取 \(log\),然后我就梦中惊坐起,然后想到了魔法少女lbw

然后拿到了 \(15pts\)。

然而太虚真人因为拍掉了魔法少女lbw,然后码了一个高精。

然后。。。

过了?!

好吧,还是 \(nb\)

就是一个非常非常基础的树形 \(dp\)

用 \(f_{i,1/0} 来表示选择这个还是不选择这个\)

其实刚开始的时候会有一个非常 \(naive\) 的想法就是按照深度奇偶进行分开,然后分别进行统计。

但是这个显然是个假的。

就好比这个图:

后面是点权,然后发现那个想法就是假的。

#include<bits/stdc++.h>
using std::cout; using std::endl;
#define try(i,a,b) for(register signed i=a;i<=b;++i)
#define throw(i,a,b) for(register signed i=a;i>=b;--i)
#define asm(i,x) for(register signed i=head[x];i;i=edge[i].next)
namespace xin_io
{
#define sb(x) cout<<#x" = "<<x<<' '
#define jb(x) cout<<#x" = "<<x<<endl
#define debug cout<<"debug"<<endl
#define gc() p1 == p2 and (p2 = (p1 = buf) + fread(buf,1,1<<20,stdin),p1 == p2) ? EOF : *p1++
char buf[1<<20],*p1 = buf,*p2 = buf; int ak; typedef long long ll; typedef unsigned long long ull;
class xin_stream{public:template<typename type>inline xin_stream &operator >> (type &s)
{
register int f = 0;s = 0; register char ch = gc();
while(!isdigit(ch)) {f |= ch == '-'; ch = gc();}
while( isdigit(ch)) s = (s << 1) + (s << 3) + (ch xor 48),ch = gc(); return s = f ? -s : s,*this;
}}io;
}
#define int long long
using namespace xin_io; static const int maxn = 1e6+10,inf = 1e9+7,mod = 1e9+7; const ll llinf = 1e18+7;
namespace xin
{
int f[maxn][2];
class xin_edge{public:int next,ver;}edge[maxn];
int head[maxn],ced = 0;
inline void add(int x,int y) {edge[++ced].ver = y; edge[ced].next = head[x]; head[x] = ced;}
int w[maxn],n;
int ind[maxn];
int d[maxn];
long double sum[maxn][2];
void dfs(int x,int fa)
{
d[x] = d[fa] + 1;
asm(i,x)
{
register int y = edge[i].ver;
if(y == fa) continue;
dfs(y,x);
f[x][1] = f[x][1] % mod * f[y][0] % mod;
sum[x][1] += sum[y][0];
f[x][0] = (f[x][0] * (sum[y][1] > sum[y][0] ? f[y][1] : f[y][0])) % mod;
sum[x][0] += sum[y][1] > sum[y][0] ? sum[y][1] : sum[y][0];
}
}
bool sp1 = 1;
inline short main()
{
io >> n;
try(i,1,n) io >> w[i],f[i][0] = 1,sum[i][0] = log(1);
try(i,1,n - 1)
{
register int x,y; io >> x >> y;
add(x,y); add(y,x);
ind[x] ++; ind[y] ++;
if(!(x == 1) and !(y == 1)) sp1 = 0;
}
if(sp1)
{
int ans = 1;
try(i,2,n) ans = ans * w[i] % mod;
cout<<ans<<endl;
return 0;
}
try(i,1,n) if(ind[i] == 1) f[i][1] = w[i],sum[i][1] = log(w[i]); else f[i][1] = w[i],sum[i][1] = log(w[i]);
dfs(1,0);
// try(i,1,n) sb(i),sb(f[i][0]),jb(f[i][1]); cout<<(sum[1][0] > sum[1][1] ? f[1][0] : f[1][1])<<endl;
return 0;
}
}
signed main() {return xin::main();}

简单题

似乎并不是很简单。

是个组合数学的题目,然后使用 \(Lucas\)

然而我居然忘了调用 \(Lucas\) 只用 \(C\) 就过了???

细节很多,看代码吧

#include<bits/stdc++.h>
using std::cout; using std::endl;
#define try(i,a,b) for(register signed i=a;i<=b;++i)
#define throw(i,a,b) for(register signed i=a;i>=b;--i)
#define asm(i,x) for(register signed i=head[x];i;i=edge[i].next)
namespace xin_io
{
#define sb(x) cout<<#x" = "<<x<<' '
#define jb(x) cout<<#x" = "<<x<<endl
#define debug cout<<"debug"<<endl
#define gc() p1 == p2 and (p2 = (p1 = buf) + fread(buf,1,1<<20,stdin),p1 == p2) ? EOF : *p1++
char buf[1<<20],*p1 = buf,*p2 = buf; int ak; typedef long long ll; typedef unsigned long long ull;
class xin_stream{public:template<typename type>inline xin_stream &operator >> (type &s)
{
register int f = 0;s = 0; register char ch = gc();
while(!isdigit(ch)) {f |= ch == '-'; ch = gc();}
while( isdigit(ch)) s = (s << 1) + (s << 3) + (ch xor 48),ch = gc(); return s = f ? -s : s,*this;
}}io;
}
#define int long long
using namespace xin_io; static const int maxn = 1e6+10,inf = 1e9+7,mod = 1e7+19; const ll llinf = 1e18+7;
namespace xin
{
inline int ksm(int x,int y)
{
register int ret = 1;
while(y)
{
if(y & 1) ret = ret * x % mod;
x = x * x % mod; y >>= 1;
}
return ret;
}
int jc[mod+10];
inline void pre_work(int ms)
{
if(ms > mod) ms = mod;
jc[0] = 1;
try(i,1,ms) jc[i] = jc[i-1] * 1ll * i % mod;
}
inline int C(int n,int m)
{
if(n < m) return 0;
return jc[n] * ksm(jc[m] * jc[n-m] % mod,mod - 2) % mod;
}
int lucas(int n,int m)
{
if(!m) return 1;
return lucas(n/mod,m/mod) * C(n%mod,m%mod) % mod;
}
int n,qnum,base = 1,fang = 1,xu = 0,num,last;
inline short main()
{
io >> n >> qnum;
if(n > mod) {try(i,1,qnum) puts("0"); return 0;}
pre_work(n); last = n;
for(register int i=0;last;++i)
{
register int temp = last - (n / base);register int zhuan = (temp + (last & 1)) / 2;
xu += zhuan * (i / 2);
if(i & 1) (num += zhuan) %= mod; else (fang *= ksm(2,zhuan)) %= mod;
last = n / base; base *= 2;
}
try(que,1,qnum)
{
register int m; io >> m;
if(m - xu < 0 or m - xu > num) puts("0");
else printf("%lld\n",lucas(num,m-xu) * fang % mod);
// sb(num);sb(xu); jb(fang);
}
return 0;
}
}
signed main() {return xin::main();}

粉丝

这个题目要求一个范围的。

我们可以转化为两个大的范围相减。

然后开始数位 \(dp\)

#include<bits/stdc++.h>
using std::cout; using std::endl;
#define try(i,a,b) for(register signed i=a;i<=b;++i)
#define throw(i,a,b) for(register signed i=a;i>=b;--i)
#define asm(i,x) for(register signed i=head[x];i;i=edge[i].next)
namespace xin_io
{
#define sb(x) cout<<#x" = "<<' '
#define jb(x) cout<<#x" = "<<endl
#define debug cout<<"debug"<<endl
#define gc() p1 == p2 and (p2 = (p1 = buf) + fread(buf,1,1<<20,stdin),p1 == p2) ? EOF : *p1++
char buf[1<<20],*p1 = buf,*p2 = buf; int ak; typedef long long ll; typedef unsigned long long ull;
class xin_stream{public:template<typename type>inline xin_stream &operator >> (type &s)
{
register int f = 0;s = 0; register char ch = gc();
while(!isdigit(ch)) {f |= ch == '-'; ch = gc();}
while( isdigit(ch)) s = (s << 1) + (s << 3) + (ch xor 48),ch = gc(); return s = f ? -s : s,*this;
}}io;
}
#define int long long
using namespace xin_io; static const int maxn = 1e6+10,inf = 1e9+7; const ll llinf = 1e18+7;
namespace xin
{
int mod,x,y,n;
int f[maxn],g[maxn],s[maxn];
inline int solve(int x)
{
if(x > n) return 0; memset(f,0,sizeof(f)); memset(s,0,sizeof(s)); memset(g,0,sizeof(g));
register int blo = std::max(x,300ll);
f[0] = g[0] = s[0] = 1;
try(i,x,blo-1) try(j,i,n) (f[j] += f[j-i]) %= mod;
try(i,1,n/blo)
{
register int t = i * blo;
try(j,i,n-t) (g[j] += g[j-i]) %= mod;
try(j,0,n-t) (s[j+t] += g[j]) %= mod;
}
int ret = 0;
try(i,0,n) (ret += f[i] * s[n-i] % mod) %= mod;
return ret;
}
inline short main()
{
io >> x >> y >> n >> mod;
printf("%lld\n",(solve(x) - solve(y+1) + mod) % mod);
return 0;
}
}
signed main() {return xin::main();}

字符串

还没过,先鸽了。

最新文章

  1. javascript-观察者模式
  2. 整块div设置为超链接进行界面跳转
  3. HDU 1693 Eat the Trees(插头DP)
  4. Python Locust性能测试框架实践
  5. Socket Programming in C#--Multiple Sockets
  6. python学习笔记之module &amp;&amp; package
  7. shell 命令
  8. 【Android 界面效果40】Android4.0-Fragment框架实现方式剖析(一)
  9. Python 基础【第九篇】运算
  10. 使用资源监控工具 glances
  11. GEF的MVC体系结构
  12. 切割图像(五)主动轮廓模型Snake简要模型
  13. win7远程连接ubuntu14.04.1桌面
  14. Spring+SpringMVC+MyBatis+easyUI整合优化篇(五)结合MockMvc进行服务端的单元测试
  15. CentOS7下安装MySQL的安装与配置(yum) (转)
  16. ElasticSearch-6.3.2 linux 安装
  17. LINUX内核分析第五周学习总结——扒开系统调用的“三层皮”(下)
  18. PWDX查找程序执行路径
  19. Docker常用命令学习
  20. DBLinq (MySQL exactly) Linq To MySql

热门文章

  1. 小马哥的 Java 项目实战营学习笔记(1)
  2. centos linux下配置固定ip,方便xshell连接
  3. 电脑常用快捷键及常用的DOS命令
  4. Vulhub-Phpmyadmin 4.8.1远程文件包含漏洞
  5. RHCE_DAY01
  6. Alibaba-技术专区-Dubbo3总体技术体系介绍及技术指南(目录)
  7. 一文让你彻底掌握ArcGisJS地图管理的秘密
  8. Windows上的暗色调Puppet书籍翻译写作环境
  9. 【原创】Spring Data Redis &lt;=2.0.3反序列化漏洞
  10. NOIP 模拟 $14\; \text{抛硬币}$