[2021.4.7多校省选模拟33]A,B,C
文章目錄
- 考試復(fù)盤
- A
- B
- C
考試復(fù)盤
今天的題其實(shí)蠻溫柔的
考完試預(yù)估分160160160,好家伙到手的只有14\frac{1}{4}41?
第一題是原題,做過(guò)的,雖然忘記怎么做了。。?!薛洹?/p>
但是因?yàn)楸旧磔^簡(jiǎn)單,考場(chǎng)上也想到了正解
但是因?yàn)槲业?span id="ozvdkddzhkzd" class="katex--inline">SBSBSB二分midmidmid附近點(diǎn)的取舍可能導(dǎo)致了答案位置的偏差111
(╥╯^╰╥)哭暈在廁所了謝謝
第二題感覺(jué)就是DPDPDP轉(zhuǎn)移,聯(lián)想到了正睿的一道劃段合并DPDPDP
但是passpasspass了,發(fā)現(xiàn)自己需要將i,i+1i,i+1i,i+1與i,i?1i,i-1i,i?1是否相鄰的信息拿到
自己的DPDPDP轉(zhuǎn)移不動(dòng)
反正就是沒(méi)想到再拿一個(gè)DPDPDP,兩個(gè)相互轉(zhuǎn)移
打了303030的暴力表,結(jié)果又出現(xiàn)了經(jīng)典菜譜——沒(méi)輸入完就直接return0return\ 0return?0
丟了
第三題撲臉而來(lái)的NTTNTTNTT味道,但是不會(huì)
只能寫最原始的DPDPDP暴力打表,結(jié)果忘記調(diào)用打表函數(shù),結(jié)果全是000
Y(>_<、)Y
總結(jié):
今天考試的題目比較簡(jiǎn)單,所以暴露出了很多細(xì)節(jié)/策略問(wèn)題,如果在省選場(chǎng)上出現(xiàn)那簡(jiǎn)直就是死亡名單
1.對(duì)細(xì)節(jié)偏差非常講究的題目(T1)一定要構(gòu)造數(shù)據(jù)想辦法測(cè)一測(cè)那個(gè)臨界附近
2.如果暴力分能打表,那么最好先寫了后放在后面自己跑,繼續(xù)往下做。今天是最后十幾分鐘打的表,時(shí)間就很緊張。如果一開(kāi)始看到就直接打表肯定是不會(huì)這么慌張的(T3)
3.對(duì)于自己不喜歡寫對(duì)拍這件事,真的要在后面每一場(chǎng)比賽中強(qiáng)迫自己寫——不一定要拍正解,至少保證自己能拿到的分一定不能丟!(T2如果寫對(duì)拍,肯定發(fā)現(xiàn)自己直接結(jié)束程序的bug)
省選在即,把握細(xì)節(jié),優(yōu)化策略,是最后的王牌!
A
考慮枚舉作為中位數(shù)的妹子
然后兩邊一定是選一段等長(zhǎng)的區(qū)間
設(shè)選的中位數(shù)在mid處
二分兩邊選多少數(shù),設(shè)選了x個(gè)
比較這2x+12x+12x+1個(gè)數(shù)的平均數(shù)和第mid?xmid-xmid?x和n?x+1n-x+1n?x+1個(gè)數(shù)的平均數(shù)
如果前者大,說(shuō)明xxx應(yīng)該變小,否則xxx變大
#include <cmath> #include <cstdio> #include <iostream> #include <algorithm> using namespace std; #define maxn 200005 int n; double a[maxn], sum[maxn];double calc( int pos, int x ) {return ( sum[pos] - sum[pos - x - 1] + sum[n] - sum[n - x] ) / ( x << 1 | 1 ); }int main() {scanf( "%d", &n );for( int i = 1;i <= n;i ++ )scanf( "%lf", &a[i] );sort( a + 1, a + n + 1 );for( int i = 1;i <= n;i ++ )sum[i] = sum[i - 1] + a[i];double ret = 0;for( int i = 1;i <= n;i ++ ) {int l = 0, r = min( i - 1, n - i ), x = 0;while( l < r ) {int mid = ( l + r ) >> 1;if( calc( i, mid ) < calc( i, mid + 1 ) )l = mid + 1;elser = mid;}ret = max( ret, calc( i, l ) - a[i] );}printf( "%.4f\n", ret );return 0; }B
設(shè)計(jì)狀態(tài)f[i][j]f[i][j]f[i][j]表示前iii個(gè)數(shù)的排列且包含jjj對(duì)非法相鄰數(shù)對(duì)的方案數(shù),當(dāng)我們放入i+1i+1i+1這個(gè)數(shù)的時(shí)候發(fā)生3種狀況:
-
打破一對(duì)非法相鄰數(shù)對(duì)
-
增加一對(duì)非法相鄰數(shù)對(duì)
-
不變。
增加的情況就發(fā)生在iii和i+1i+1i+1相鄰的情況
所以我們還需要知道iii是否處于某個(gè)非法相鄰數(shù)對(duì)。
因此重新設(shè)計(jì)狀態(tài)
f[i][j]f[i][j]f[i][j]:前i個(gè)數(shù)的排列且包含j對(duì)非法相鄰數(shù)對(duì),iii存在于非法相鄰數(shù)對(duì)中的方案數(shù)
g[i][j]g[i][j]g[i][j]:iii不存在于非法相鄰數(shù)對(duì)中的方案數(shù)
轉(zhuǎn)移就是考慮i+1i+1i+1的放置方法
f[i][j]f[i][j]f[i][j] 的放置方法有4種:
-
與i相鄰且增加一對(duì)非法相鄰數(shù)對(duì),只會(huì)是[i?1,i,i+1],1[i-1,i,i+1], 1[i?1,i,i+1],1種方法
f[i+1][j+1]f[i+1][j+1]f[i+1][j+1]
-
與iii相鄰且不改變非法相鄰數(shù)對(duì)數(shù),只會(huì)是[i?1,i+1,i],1[i-1,i+1,i], 1[i?1,i+1,i],1種方法
f[i+1][j]f[i+1][j]f[i+1][j]
-
不與iii相鄰且減少一對(duì)非法相鄰數(shù)對(duì),有j?1j-1j?1種方法
g[i+1][j?1]g[i+1][j-1]g[i+1][j?1]
-
不與i相鄰且不改變非法相鄰數(shù)對(duì)數(shù),有i?ji-ji?j種方法
g[i+1][j]g[i+1][j]g[i+1][j]
g[i][j]g[i][j]g[i][j]有333種情況
-
與iii相鄰且增加一對(duì)非法相鄰數(shù)對(duì),可能是[i+1,i][i+1,i][i+1,i]或[i,i+1],2[i,i+1],2[i,i+1],2種方法
f[i+1][j+1]f[i+1][j+1]f[i+1][j+1]
-
不與iii相鄰且減少一對(duì)非法相鄰數(shù)對(duì),有jjj種擺放方法
g[i+1][j?1]g[i+1][j-1]g[i+1][j?1]
-
不與iii相鄰且不改變非法相鄰數(shù)對(duì)數(shù),有i?j?1i-j-1i?j?1種擺放方法
g[i+1][j]g[i+1][j]g[i+1][j]
最終的答案就是g[N][0]g[N][0]g[N][0]
#include <cstdio> #include <cstring> #define int long long #define maxn 1005 int n, mod; int f[maxn][maxn], g[maxn][maxn];signed main() {while( ~ scanf( "%lld %lld", &n, &mod ) ) {memset( f, 0, sizeof( f ) );memset( g, 0, sizeof( g ) );g[1][0] = 1;for( int i = 1;i < n;i ++ )for( int j = 0;j < i;j ++ ) {f[i + 1][j + 1] = ( f[i + 1][j + 1] + f[i][j] ) % mod;f[i + 1][j] = ( f[i + 1][j] + f[i][j] ) % mod;if( j ) g[i + 1][j - 1] = ( g[i + 1][j - 1] + f[i][j] * ( j - 1 ) ) % mod;g[i + 1][j] = ( g[i + 1][j] + f[i][j] * ( i - j ) ) % mod;f[i + 1][j + 1] = ( f[i + 1][j + 1] + g[i][j] * 2 ) % mod;if( j ) g[i + 1][j - 1] = ( g[i + 1][j - 1] + g[i][j] * j ) % mod;g[i + 1][j] = ( g[i + 1][j] + g[i][j] * ( i - j - 1 ) ) % mod; }printf( "%lld\n", ( g[n][0] + mod ) % mod );}return 0; }C
令N=3×7×11×47=10857N=3\times 7\times 11\times 47 =10857N=3×7×11×47=10857
fi,jf_{i,j}fi,j?表示前iii位模NNN為jjj的方案數(shù)
利用倍增思想,優(yōu)化暴力轉(zhuǎn)移
f2i,(j?10i+k)%N=∑fi,j?fi,kf_{2i,(j·10^i+k)\%N}=\sum f_{i,j}·f_{i,k}f2i,(j?10i+k)%N?=∑fi,j??fi,k?
固定2i2i2i,令k=10i,gt=∑j?k%N=tfi,j,ht=fi,tk=10^i,g_t=\sum_{j·k\%N=t}f_{i,j},h_t=f_{i,t}k=10i,gt?=∑j?k%N=t?fi,j?,ht?=fi,t?
則轉(zhuǎn)移變?yōu)?span id="ozvdkddzhkzd" class="katex--inline">f2i,j+k=gj?hkf_{2i,j+k}=g_j*h_kf2i,j+k?=gj??hk?
標(biāo)準(zhǔn)可愛(ài)的卷積形式,FFTFFTFFT加速矩陣
#include <cmath> #include <cstdio> #include <iostream> using namespace std; #define ll long long #define maxn 200100 #define mod 9973 #define N 10857struct complex {double x, i;complex(){}complex( double X, double I ) {x = X, i = I;} }A[maxn], B[maxn];double pi = acos( -1.0 );complex operator + ( complex a, complex b ) {return complex( a.x + b.x, a.i + b.i ); }complex operator - ( complex a, complex b ) {return complex( a.x - b.x, a.i - b.i ); }complex operator * ( complex a, complex b ) {return complex( a.x * b.x - a.i * b.i, a.x * b.i + a.i * b.x ); }int len = 32768; int r[maxn];void FFT( complex *v, int opt ) {for( int i = 0;i < len;i ++ )if( i < r[i] ) swap( v[i], v[r[i]] );for( int i = 1;i < len;i <<= 1 ) {complex omega( cos( pi / i ), opt * sin( pi / i ) );for( int j = 0;j < len;j += ( i << 1 ) ) {complex w( 1, 0 );for( int k = 0;k < i;k ++, w = w * omega ) {complex x = v[j + k], y = v[j + k + i] * w;v[j + k] = x + y;v[j + k + i] = x - y;}}} }int qkpow( int x, int y ) {int ans = 1;while( y ) {if( y & 1 ) ans = ans * x % N;x = x * x % N;y >>= 1;}return ans; }bool check( int x ) {return x % 2 && x % 3 && x % 5 && x % 7 && x % 11 && x % 47; }void mul( int *f, int *g, int L ) {int k = qkpow( 10, L );for( int i = 0;i < N;i ++ ) {A[i * k % N] = complex( int( f[i] + A[i * k % N].x ) % mod , 0 );B[i] = complex( g[i], 0 );f[i] = 0;}FFT( A, 1 ), FFT( B, 1 );for( int i = 0;i < len;i ++ ) A[i] = A[i] * B[i];FFT( A, -1 );for( int i = 0;i < len;i ++ ) {f[i % N] = ( f[i % N] + (ll)( A[i].x / len + 0.5 ) ) % mod;A[i] = B[i] = complex( 0, 0 );} }int n; int c[5] = { 1, 2, 3, 5, 7 }; int f[maxn], g[maxn];int main() {scanf( "%d", &n );if( n == 1 ) return ! printf( "1\n" );else n --;int L = f[0] = g[1] = g[2] = g[3] = g[5] = g[7] = 1;int l = log2( len );for( int i = 0;i < len;i ++ )r[i] = ( r[i >> 1] >> 1 ) | ( ( i & 1 ) << ( l - 1 ) );while( n ) {if( n & 1 ) mul( f, g, L );mul( g, g, L );n >>= 1;L <<= 1;}int ans = 0;for( int i = 0;i < N;i ++ )for( int j = 0;j < 5;j ++ )if( check( i * 10 + c[j] ) )ans = ( ans + f[i] ) % mod;else;printf( "%d\n", ans );return 0; }總結(jié)
以上是生活随笔為你收集整理的[2021.4.7多校省选模拟33]A,B,C的全部?jī)?nèi)容,希望文章能夠幫你解決所遇到的問(wèn)題。
- 上一篇: 2021-4-4 省选模拟赛(灯,十字路
- 下一篇: [CF/AT]各大网站网赛 体验部部长第