2053:?Birthday

樣例輸入

3 3 1 2 1 2 1 2

樣例輸出

5

分析：

考慮費用流時把每個part拆成n個點，選擇第i個點的代表為放置i塊蛋糕和(i - 1)塊蛋糕的時間差，這個時間差是遞增的，因此在費用流的過程中必定會從小到大選擇

具體建圖：左邊n個點代表n個蛋糕，右邊m * n個點代表m個part,每個part拆成n個點。源點向每個左邊的點連一條流量1費用0的邊，每個右邊的點向匯點連一條流量1費用0的編。每個蛋糕向可以放的兩個part的所有點連邊，連向第i個點的費用為i^2 - (i - 1)^2，流量為1。這樣求最小費用流既為答案。

#include <cstdio> #include <cstring> #include <vector> #include <queue>using namespace std;const int INF = 0x3f3f3f3f; const int MAXN = 5005; struct Edge{int value,flow,to,rev;Edge(){}Edge(int a,int b,int c,int d):to(a),value(b),flow(c),rev(d){} };vector<Edge> E[MAXN];inline void Add(int from,int to,int flow,int value){E[from].push_back(Edge(to,value,flow,E[to].size()));E[to].push_back(Edge(from,-value,0,E[from].size()-1)); }bool book[MAXN];//用于SPFA中標記是否在queue中 int cost[MAXN];//存費用的最短路徑 int pre[MAXN];//存前節點 int pree[MAXN];//存在前節點的vector中的下標 bool Spfa(int from,int to){memset(book,false,sizeof book);memset(cost,INF,sizeof cost);book[from] = true;cost[from] = 0;queue<int> Q;Q.push(from);while(!Q.empty()){int t = Q.front();book[t] = false;Q.pop();for(int i=0 ; i<E[t].size() ; ++i){Edge& e = E[t][i];if(e.flow > 0 && cost[e.to] > cost[t] + e.value){cost[e.to] = cost[t] + e.value;pre[e.to] = t;pree[e.to] = i;if(book[e.to] == false){Q.push(e.to);book[e.to] = true;}}}}return cost[to] != INF; }int Work(int from,int to){int sum = 0;while(Spfa(from,to)){int mflow = INF;//SPFA找到的最短路徑的最小容量 int flag = to;while(flag != from){mflow = min(mflow,E[pre[flag]][pree[flag]].flow);flag = pre[flag];}flag = to;while(flag != from){sum += E[pre[flag]][pree[flag]].value * mflow;E[pre[flag]][pree[flag]].flow -= mflow;E[flag][E[pre[flag]][pree[flag]].rev].flow += mflow;flag = pre[flag];}}return sum; }int main(){int N,M;while(scanf("%d %d",&N,&M) == 2){int a,b;for(int i=1 ; i<=N ; ++i){scanf("%d %d",&a,&b);Add(i,a+N,1,0);Add(i,b+N,1,0);Add(0,i,1,0); }for(int i=1 ; i<=M ; ++i){for(int j=1 ; j<=99 ; j+=2){//99 = 2*50-1;Add(i+N,M+N+1,1,j);}} printf("%d\n",Work(0,M+N+1));for(int i=0 ; i<=M+N+1 ; ++i)E[i].clear(); }return 0; }

2088: 電音之王

樣例輸入

1 1 1 1 1 0 1000000007 10

樣例輸出

904493530

分析：

由于取模的常數為64bit整數，所以考慮使用O(1)快速乘，然后暴力求解即可

#include<bits/stdc++.h> using namespace std; typedef long long ll; ll ksc(ll x,ll y,ll mod) {return (x*y-(ll)((long double)x/mod*y)*mod+mod)%mod; } ll m0,m1,c,M,k; ll a0,a1; int main() {int T;scanf("%d",&T);while(T--){scanf("%lld%lld%lld%lld%lld%lld%lld",&a0,&a1,&m0,&m1,&c,&M,&k);ll ans=ksc(a0,a1,M);ll a2;for(int i=2;i<=k;i++){a2=((ksc(m0,a1,M)+ksc(m1,a0,M))%M+c)%M;ans=ksc(ans,a2,M);a0=a1;a1=a2;}printf("%lld\n",ans);}return 0; }

鏈接：https://www.nowcoder.com/acm/contest/205/B
來源：牛客網
?

電音之王

時間限制：C/C++ 2秒，其他語言4秒
空間限制：C/C++ 262144K，其他語言524288K
64bit IO Format: %lld

題目描述

終于活成了自己討厭的樣子。
聽說多聽電音能加快程序運行的速度。
定義一個數列，告訴你a0,a1,m0,m1,c，定義an=m0an-1+m1an-2+c對所有n≥ 2。
求

輸入描述:

第一行一個整數T(1≤ T≤ 1000)，表示數據組數。 每組數據一行7個整數a0,a1,m0,m1,c,M,k，保證1≤ M≤ 1018,0≤ a0,a1,m0,m1,c< M, 2≤ k≤ 106，保證M為奇數。 保證。

輸出描述:

對于每組數據，輸出一行表示答案。

示例1

輸入

復制

1 1 1 1 1 0 1000000007 10

輸出

復制

904493530 #include <stdio.h> #include<vector> #include<utility> #include <assert.h> using namespace std; #define rep(i,a,n) for (int i=a;i<n;i++) #define per(i,a,n) for (int i=n-1;i>=a;i--) #define pb push_back #define mp make_pair #define all(x) (x).begin(),(x).end() #define fi first #define se second #define SZ(x) ((int)(x).size()) typedef vector<int> VI; typedef long long ll; typedef pair<int,int> PII; const ll mod=1000000007;ll gcd(ll a,ll b) { return b?gcd(b,a%b):a;} // headtypedef unsigned long long u64; typedef __int128_t i128; typedef __uint128_t u128; int _,k; u64 A0,A1,M0,M1,C,M;struct Mod64 {Mod64():n_(0) {}Mod64(u64 n):n_(init(n)) {}static u64 init(u64 w) { return reduce(u128(w) * r2); }static void set_mod(u64 m) {mod=m; assert(mod&1);inv=m; rep(i,0,5) inv*=2-inv*m;r2=-u128(m)%m;}static u64 reduce(u128 x) {u64 y=u64(x>>64)-u64((u128(u64(x)*inv)*mod)>>64);return ll(y)<0?y+mod:y;}Mod64& operator += (Mod64 rhs) { n_+=rhs.n_-mod; if (ll(n_)<0) n_+=mod; return *this; }Mod64 operator + (Mod64 rhs) const { return Mod64(*this)+=rhs; }Mod64& operator -= (Mod64 rhs) { n_-=rhs.n_; if (ll(n_)<0) n_+=mod; return *this; }Mod64 operator - (Mod64 rhs) const { return Mod64(*this)-=rhs; }Mod64& operator *= (Mod64 rhs) { n_=reduce(u128(n_)*rhs.n_); return *this; }Mod64 operator * (Mod64 rhs) const { return Mod64(*this)*=rhs; }u64 get() const { return reduce(n_); }static u64 mod,inv,r2;u64 n_; }; u64 Mod64::mod,Mod64::inv,Mod64::r2;int main() {for (scanf("%d",&_);_;_--) {scanf("%llu%llu%llu%llu%llu%llu%d",&A0,&A1,&M0,&M1,&C,&M,&k);Mod64::set_mod(M);Mod64 a0(A0),a1(A1),m0(M0),m1(M1),c(C),ans(1),a2(0);for (int i=0;i<=k;i++) {ans=ans*a0;a2=m0*a1+m1*a0+c;a0=a1; a1=a2;}printf("%llu\n",ans.get());}

?

總結

以上是生活随笔為你收集整理的Wannafly summer camp的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。