序
考试爆炸祭…暴力分那么高，打什么正解
结果
题面
T1
T2
T3
题解
T1
一眼二分答案，结果写错了爆5分
T2
一眼看到这个:
对于所有编号为奇数的数据:$a_i=0, i∈[0,7]$;
于是大力贪心骗到了$65pts$
正解差分约束？？？我还看到过？？？
T3
我要好好写一下这一题的题解暴力80分我还来想正解，SHIFT
最近学了一个叫做$Generating\;Function$的东西
于是对于一个区间$[l, r]$，有
$$G(x)=\prod_{i=l}^r(a_ix+1)$$
因为$k<=10$
所以可以维护母函数前$10$项系数
$pushup$很好写，就是一个多项式乘法
$pushdown$的取反很好写，只要考虑奇数项
区间加则不那么好写，(毒瘤线段树)
令$[x^k]$表示$G(x)$中$x^k$的系数
$$\therefore[x^k] = \sum_{i_1<i_2<i_3<…<i_k}\prod_{j=1}^ka_{i_j}$$
$$\because new[x^k] = \sum_{i_1<i_2<i_3<…<i_k}\prod_{j=1}^k(a_{i_j} + c)$$
$\therefore$对于每一个展开式中的$l$
$\prod_{j=1}^la_{i_j}$与$\prod_{j=1}^{k-l}a_{i_j^{,}}$一一对应
又$\because \prod_{j=1}^{k-l}a_{i_j^{,}}$有$\binom{len-l}{k-l}$个
$\therefore$这里贡献的答案为$c^{k-l}\binom{len-l}{k-l}\prod_{j=1}^la_{i_j}$
$$\therefore new[x^k]=\sum_{j=0}^{k-1}\binom{len-j}{k-j}[x^j]c^{k-j}$$
于是就可以$pushdown$了
注意指针的问题!!!
代码
T1
最后AC的代码让我自己都懵1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47#include<bits/stdc++.h>
#define RG register
#define file(x) freopen(#x".in", "r", stdin);freopen(#x".out", "w", stdout);
#define clear(x, y) memset(x, y, sizeof(x));
using namespace std;
inline int read()
{
int data=0, w=1;
char ch=getchar();
while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar();
if(ch=='-') w=-1, ch=getchar();
while(ch>='0'&&ch<='9') data=(data<<3)+(data<<1)+(ch^48), ch=getchar();
return data*w;
}
const int maxn(1e5+10);
int x[maxn], y[maxn], n, m, qx, qy;
inline bool check(int mid)
{
long long a = 1ll * y[mid] * qx;
long long b = 1ll * x[mid] * qy;
long long c = 1ll * x[mid] * y[mid];
return c <= (a + b);
}
int main()
{
n = read();
for(RG int i=1;i<=n;i++) x[i]=read();
for(RG int i=1;i<=n;i++) y[i]=read();
sort(x+1, x+n+1); sort(y+1, y+n+1);
m = read();
while(m--)
{
qx = read(); qy = read();
int l = 0, r = n;
while(l<r)
{
int mid(l+r>>1);
if(check(mid + 1)) l = mid + 1;
else r = mid;
}
printf("%d\n", l);
}
return 0;
}
T2
没有
T3
1 | #include<bits/stdc++.h> |
是不是很短