<output id="qn6qe"></output>

    1. <output id="qn6qe"><tt id="qn6qe"></tt></output>
    2. <strike id="qn6qe"></strike>

      亚洲 日本 欧洲 欧美 视频,日韩中文字幕有码av,一本一道av中文字幕无码,国产线播放免费人成视频播放,人妻少妇偷人无码视频,日夜啪啪一区二区三区,国产尤物精品自在拍视频首页,久热这里只有精品12

      分式分解——2025.6.8 鮮花

      分式分解——2025.6.8 鮮花

      余命2:30
      あ 生まれた
      何をして生きようか
      無限の可能性が広がる まだ8秒だ
      あんなことできるかな
      何にでもなれるかな
      來るべき不幸から
      目を背けながら
      真面目に生きようか
      小狡く生きようか
      君を助けたい 蹴落としたい
      どっちも自分だ
      夢や戀はファストで
      努力はつらいよね
      本や映畫は
      あらすじだけでいいかな
      余命2:30
      余命2:30
      「ありがとう」「さよなら」
      「またね」をくり返して
      もう45秒の春
      余命2:30
      使い捨てられる青
      3分に満たない
      替えが利く ぼくの命を
      涙で消費しないでね
      夢ってこんなもんか
      戀ってこんなもんか
      人生ってまるで
      よくある歌みたいだな
      1分以上過ぎて
      「自分」を思い知って
      だけど 悟るには まだまだ早いね
      チクタク チクタク
      もういない君を想う
      1分前 描いた夢とは ほど遠いが
      余命2:30
      余命2:30
      「よかった」「もうだめだ」
      「まあいっか」をくり返して
      1分45秒の秋
      余命2:30
      また一つ消える魔法
      君のため?ぼくのため?
      大袈裟に嘆く命を
      どうか美化しすぎないでね
      あの日 読み飛ばした
      本や映畫みたいに
      ぼくの命も 短い歌になって
      余命2:30でさよならでも
      幸せだったと言わせて欲しいの
      余命2:30
      余命2:30
      3分に満たない
      もう終わる ぼくの命を
      人事って思わないでね
      あ
      

      比較易。

      未加說明則 \(x\) 是多項式的變量。

      設分式形如 \(F = \frac 1{\prod\limits_i Q_i}\),我們想把他分解成 \(\sum\limits_i \frac{R_i}{Q_i}\),其中 \(\deg R_i < \deg Q_i\)

      有:

      \[\frac 1{\prod_i Q_i}=\sum_i \frac{R_i}{Q_i} \]

      \[1=\sum_i \frac{\prod_j Q_j R_i}{Q_i} \]

      考慮 \(\deg\) 的限制,我們對 \(Q_i\) 取模。

      \[R_i=(\prod_{j != i} Q_j)^{-1} \bmod {Q_i} \]

      這個東西用處只在于其形式,用來做類似 EI 營業日志 2020.5.20

      正常的實現這個感覺都 \(n^2\log^2 n\sim n^3\) 了。

      Upd:jijidawang 說 EI 已經 \(n\log^2 n\) 了,不過我不會

      我們在 \(Q_i = 1 - a_ix\) 時給出一個更正常的 \(n ^ 2\) 做法。

      考慮取模的意義,實際上是去掉其他項的貢獻,我們用一個類似拉差的做法做到這個。

      設原式

      \[\sum\limits_i \frac{R_i}{1 - a_ix} \]

      通分后 \(R_i\) 的限制

      \[\sum\limits_i R_i \prod\limits_{j != i} (1 - a_jx) = 1 \]

      我們取 \(x = \frac 1{a_k}\),容易發現此時所有 \(i \not= k\) 的項全等于 \(0\),所以此時有

      \[R_i \prod\limits_{j != i} (1 - \frac{a_j}{a_i}) = 1 \]

      即可解出

      \[R_i = (\prod\limits_{j != i} (1 - \frac{a_j}{a_i})) ^ {-1} \]

      例題:AT_abc241_h [ABC241Ex] Card Deck Score

      直接整出生成函數

      \[[x^m] \prod_{i = 1} ^ n \frac{1 - (a_ix_i)^{b_i}}{1 - a_ix_i} \]

      考慮分子的系數直接 \(\mathcal{O}(n2^n)\) 枚舉即可,重點在于分母的系數。

      我們直接套用分式分解,將

      \[\prod_{i = 1} ^ n \frac 1{1 - a_ix_i} \]

      改寫成

      \[\sum\limits_i \frac{R_i}{1 - a_ix} \]

      則可以 \(\mathcal{O}(n\log P)\)\(P = 998244353\) 為模數)求單次系數,將枚舉分子的 \(2^n\) 個系數帶入即可。

      總復雜度是 \(\mathcal{O}(n2^n\log P)\) 的。

      Code
      /*
      clear && g++ % -o %< -O2 -std=c++14 -DLOCAL -Wall -Wextra && time ./%< && size %<
      clear && g++ % -o %< -O2 -std=c++14 -DLOCAL -Wall -Wextra -fsanitize=address,undefined -g && time ./%< && size %<
      echo && cat in_out/out.out && echo
      */
      #include <bits/stdc++.h>
      using namespace std;
      using llt = long long;
      using llf = long double;
      using ull = unsigned long long;
      #define endl '\n'
      #ifdef LOCAL
      FILE *InFile = freopen("in_out/in.in", "r", stdin), *OutFile = freopen("in_out/out.out", "w", stdout);
      #endif
       
      const int N = 20, MOD = 998244353;
      int n; llt _a[N], _b[N], tt[N], m;
      int Fpw(int a, int b){
      	int ans = 1;
      	for(; b; a = 1ll * a * a % MOD, b >>= 1) if(b & 1) ans = 1ll * ans * a % MOD;
      	return ans;
      }
      int Inv(int a){ return Fpw(a, MOD - 2); }
      
      int main(){
      	ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
      	cin >> n >> m;
      	for(int i = 1; i <= n; ++i) cin >> _a[i] >> _b[i], _b[i] += 1, _a[i] %= MOD;
      	for(int i = 1; i <= n; ++i){
      		int s = 1, v = Inv(_a[i]);
      		for(int j = 1; j <= n; ++j) if(j != i) s = s * (1 - v * _a[j] % MOD) % MOD;
      		tt[i] = Inv(s);
      	}
      	llt ans = 0;
      	for(int i = 0, cs = (1 << n) - 1; i <= cs; ++i){
      		llt c = 0, t = 1;
      		for(int j = 1; j <= n; ++j) if((i >> (j - 1)) & 1)
      			c += _b[j], t = 1ll * t * -Fpw(_a[j], _b[j] % (MOD - 1)) % MOD;
      		c = (m - c) % (MOD - 1); if(c < 0) continue;
      		for(int j = 1; j <= n; ++j) ans = (ans + 1ll * t * Fpw(_a[j], c) % MOD * tt[j]) % MOD;
      	}
      	cout << (ans % MOD + MOD) % MOD;
      }
      
      P




      posted @ 2025-06-08 21:48  xrlong  閱讀(44)  評論(2)    收藏  舉報

      Loading

      主站蜘蛛池模板: 无码日韩做暖暖大全免费不卡| 亚洲av色在线观看国产| 欧美嫩交一区二区三区| 日本一二三区视频在线| 免费国产精品黄色一区二区| 亚洲中文字幕综合小综合| 欧美日韩免费专区在线观看| 国产一区二区一卡二卡| 国产极品粉嫩学生一线天| 大地资源中文第二页日本| 久久久久久曰本av免费免费| 草草浮力影院| 九九久久亚洲精品美国国内| 亚洲中文字幕日产无码成人片| 清河县| 国产av综合影院| 国产对白老熟女正在播放| 久久发布国产伦子伦精品| 18禁无遮挡啪啪无码网站| 亚洲一区二区三区水蜜桃| 欧美日韩亚洲国产| 国产SM重味一区二区三区| 精品国产一区二区色老头| 蜜桃网址| 少妇高潮水多太爽了动态图| 日本一级午夜福利免费区| 亚洲欧美高清在线精品一区二区| 亚洲熟妇无码另类久久久| 亚洲 欧美 中文 日韩aⅴ| 久女女热精品视频在线观看| 无码任你躁久久久久久久| 久久国产精品精品国产色婷婷| 免费一区二三区三区蜜桃| 精品人妻系列无码天堂| 亚洲国产亚洲综合在线尤物| 国产伦一区二区三区久久| 亚洲一区二区在线无码| 99亚洲男女激情在线观看| caoporn免费视频公开| 美乳丰满人妻无码视频| 桑植县|