https://atcoder.jp/contests/abc339

F - Product Equality

n个数字a[n], 问有多少个有序对(i,j,k)满足

A[i]*A[j]=A[k]

n 1000

$a[i] \in [0,10^{10}]$

2s

1024mb

我的思路

光是 枚举 (i,j) 那么也是 n^2了，那么乘法和对比就要足够快, 然而最快的就算ntt，也是 (长度 log(长度)) 会tle

一个想法 是 通过减少比较范围，但如果正好 一半长度 [len/2] 一半长度[len]

第二个想法是 选取一些 质数p,然后通过类似hash的想法来完成

a[i] * a[j] == a[k]

那么 (a[i] % p)*(a[j] % p) % p== a[k] %p

没有任何确定性算法

然后就是 十进制压8位，但这样也只是降低到 int[125] x int[125]

n的个数没有减少

Kirchhoff’s theorem(基尔霍夫定律)

无自环，可重边无向图$G$的生成树个数 = G的Laplacian矩阵(基尔霍夫矩阵)的$n-1$阶余子式的行列式

https://atcoder.jp/contests/abc338

G - evall

给定字符串S,由123456789+*组成

最开始和最后的字符是数字，没有相邻的非数字

对于$(i\le j)$, 如果s[i],s[j]都是数字，则eval(i,j)为这一段字符串的表达式结果

否则 为0

求 $\displaystyle \sum_{i=1}^{|S|}\sum_{i=j}^{|S|} eval(i,j)\mod 998244353$

$|S| \le 10^6$

2s

1024mb

我的思路

这里虽然没有括号，但是有乘法的优先运算

如果固定前面

1
2
3

1234123*234*345+3245*145+324*45*345
[............]
fix          i->

那我们得到的是 [base + prodbase * (cur=34)]

每次遇到数字就是cur*=10,cur+=int(s[i])

而进入新的乘法则是prodbase*=cur, cur = 0

而进入新的加法则是base+=prodbase*cur,prodbase=1,cur=0

这个方法的问题在于 如果想用矩阵表示，那么在做prodbase * cur的时候 无法做这个乘法

再增加prodres东西

1
2
3
4
5
6

1234123*234*345+3245*145+324*45*345
                                 i
[base                  ]          
                         xxxxxx      prodbase 
                                xx   cur
                         xxxxxxxxx   prodres

ans = base + prodres

prodres = prodbase * cur, 而新的转移 发现cur不需要了

那么 遇到数字时,

1
2
3
4
5
6

			sumall     base prodres    prodbase    1
sumall        1 
base          1         1     
prodres      10                 10
prodbase   int(char)         int(char)     1  
1                                                  1

那么 遇到+时,

1
2
3
4
5
6

			sumall    base prodres    prodbase     1
sumall        1 
base                   1          
prodres                1                
prodbase            
1                                         1        1

那么 遇到*时,

1
2
3
4
5
6

			sumall    base prodres    prodbase    1
sumall        1 
base                   1          
prodres                                   1               
prodbase            				       
1                                                 1

所以，后缀矩阵乘法就秒了

前置的是[0,0,0,1,1]

Hall’s theorem, 霍尔定理

二分图 左侧$n$点,右侧$m$点, $n\le m$

二分图的最大匹配个数$=n$ 的充要条件, 左侧点$n$的任意大小$(=k)$的子集 连到右图的点的个数都满足$\ge k$

$f(a,b,c,n) = \sum_{x=0}^n \lfloor \frac{ax+b}{c} \rfloor$

代码

1
2
3
4
5
6
7

// \sum_{x=0}^n \lfloor \frac{ax+b}{c} \rfloor
ll floor_sum(ll a,ll b,ll c,ll n){
  if(a==0) return (b/c)*(n+1);
  if(a >= c or b >= c) return n*(n+1)/2*(a/c) + (n+1)*(b/c) + floor_sum(a%c,b%c,c,n);
  ll m = (a*n+b)/c;
  return m*n - floor_sum(c,c-b-1,a,m-1);
}